STATISTICAL ANALYSIS OF SEISMIC DATA 

AND 

SEISMIC RISK ANALYSIS OF INDIAN PENINSULA 


A Thesis Submitted 

In Partial Fulfilment of the Requirements 
for the Degree of 

DOCTOR OF PHILOSOPHY 


by 

SUSANTA BASU 

I 


to the 

DEPARTMENT OF CIVIL ENGINEERING 

INDIAN INSTITUTE OF TECHNOLOGY KANPUR 

JUNE 1977 



fir, 

Acc, ssfe,. 



B HAY tm. 



CE- i‘57-7--D'-(fcAS- STA 



CERTIFICATE 


This is to certify that the thesis entitled, 
'^Statistical Analysis of Seismic Data and Seismic Risk 
Analysis cf Indian Peninsula" by Suscnta Baeu, has been 
carried out under my supervision and has mt been mb- 
mitted elsewhere for a degree. 


dune - 1977 



Department of Aeronmtioal Engineering 
Indian Institute of Technology, Kanpur 
imiA 


?osT graouatk orriGE 
Th« has been 
feif the »w«r«l the tiegree «>f 
Doctor nf Pi:iiif..3k>pfsy 
in ar"* ' with the ■ | 

regulatimts ol the ituiijw 
{astiivite of ^ 

Dated; %S ' ^ • \ 97f^ ^ R^__ \ 


iiCKNOTLEDGElIEffiOS 


The author expresses his deep sense of gratitude to Prof. N.C. Nigan. 
for his able guidance and constant encouragement during the course of this 
work . 

The author is indebted to litr. H.S. Chaudhury, The Director 
Seismology Division, Indian Meteriological Department, New Delhi, for 
permitting the author to use Earthquake Catalogue of I.MJ), The author 
thankfully acknowledges the help provided hy Prof. A.S. iirya and 
Dr. P.E. j!\garwal of SRTEE, Roorkee during data collection. 

The author is thankful to llir, B, Ghosh, Librarian, Geological 
Survey of India, Calcutta, for extending the library facilities to the 
author . 

The author expresses his sincere gratitude to I\fcs. Katharine Rao, 
who translated the Spanish references cited in this vrork. 

Part of the financial assistance to the author came from H.B.T.I., 
Kanpur and the author is grateful to the authorities of the H.B.T.I. 

The author is grateful to his colleagues in H.B.T.I. in general 
and to Dr. C .V .S.K. Rao and Mr. H.S. Niranjan in particular for their 
generous help and cooperation. 

The author acknowledges the help provided by Mr. T.M.S. Raghavan, 
Mr. G.C. Baral, Dr. S.K. Ray and Mr. S. Banerjee for helpful discussion 
and correction of manuscript. 



IV 


The author is very muoh thankful to Ifc. R.C. udhikari for 
his help in various ways during the course of the work. 

The author is thankful to the staff of the Computer Centre, 

I. I. T , /Kanpur, for their cooperation. 

Finally, thanks are due to If . G.L. Misra and Me. S.K. Tewari 
for their skillful and excellent typing. 



r 


LIST OF TABLES 


CHAPTER 2 s STATISTICAL ANALYSIS OF SEISMIC DiiTA OP 
BIDIi'iN PENINSULA 


Table 2 .1 s Properties of Earthquake with Magnitude 

above Five 25 

Table 2.2 s Estimated Serial Autocorrelation 

Coefficients 25 

Table 2.3 ° Estimated Correlation Matrix 25 

Table 2 4 ° Mean Squares for Various Block Sizes in 

Grieg-Smith' s Pattern Analysis of Epicentral 
Location 27 

Table 2.5 s Analysis of Variance Table for Linear 

Regression P.n N = a + bM 55 

Table 2 .6 s Hyperexponential Distribution of Inter- 
arrival Time 43 

Table 2.7 • State of Markov Chain, its Stationary 

Probabilities and Parametersof Holding 
Time Distribution 44 

Table 2.8 s Transition Probability Ifetrix of Markov 

Chain 47 


CHAPTER 3 s MODULliR EARTHQUAKE SOURCE ilND SENSITIVITY 

ANiiLYSIS 

Table 3*1 « Ratio of Return Periods for Normal (^) and 

Mxed Lognoraml (mfi.) Distribution to Focal 
Depth, H 69 

CH/iPTER 4 8 SEISfflC RISK ANALYSIS OP BDLiN PENINSULA 

A MODEL HOMOGENEOUS IN TECTONIC FEATURES 

Table 4 *1 ° Analysis of Variance Table for Multiple 

Linear Regression Model - 
C(v) = exp [ aQ+ a^ Hn M + a^M ] 


87 




Ti 


Tatle 4 .2 


Tatle 4 *5 


CHiiPTER 5 


Table 5*1 


CHAPTER 6 


Table 6.1 
Table 6.2 


CHAPTER 7 
Table 7 *1 


Table 7 *2 


APPEHDIX A7.1 
Table A7 .1 


Analysis of Variance Table for Linear 
Regression Model - 

S,n[n { 1+C^(y' )}/(nd1 ) ] =Y ^n(V/V< ) 90 

Estimated Intensity and Coefficient of 
Variation Sq^uare for Earthquake and Peak 
Generali2ed Seismic Intensity at Grid 
Points. Lognormal Focal Depth Distribution, 

HI E (5f9 ] f b = 150 kms. and H^=600 kms* 93-101 

SEISMIC RISK analysis OP INDIilN PENINSHLA- 
A MODEL NOffiOMOGENEOUS IN TECTONIC PEATHEES 

Geometric and Associated Seismic Activity 

of Tectonic Dnits and Faults of Indian 

Peninsula 121-122 

EFFECT OF SCATTER IN THE ATTENUATION LAW 
ON SEISMIC RISK 

Parameter of Envelopes 157 


Values of C^ and for Different 

Probability Levels for Envelopes A and B 158 

SEISMC ZOHNING lA'iPS OF IlffilA 
Comparison of IS s18'93-1975 with Present 


Work 


158 

Peak Ground Acceleration in g for 
Different Cities 

159 

LIST OP EARTHQU/iKEB i^PCEE 

1917 


List of Earthquake Before 

1917 

'250 



CHAPTER 

Pig. 

Pig. 

Pig. 

Pig. 

Pig. 

Pig. 

CHAPTER 

Pig. 

Pig. 

Pig. 

Pig. 

Pig. 


Tii 


LIST OP PIGHEffiS 


Page Ho . 

2 5 STATISTICAL ANALYSIS CP SEISMIC DATA 

OP INDIAN PENINSULA 

2 .1 s Number of EaxtHcu-aAes in each Year 

with I’ifagnitude > 5 

2.2 ; Histogram of Spatial Occurrence of 

Earthquakes of Magnitude, M > 5, from 
January 1917 to December 1975> total 
Number 1575 

2 .3 : Histogram of Pocal Depth with Magnitude, 

M > 5 

2 4 " Frequency-Magnitude Relation 

2 .5 s Seismicity of Indian Peninsulas Energy 
Contours in 10^5 ergs km"^ yr”"' 

2 .6 5 Arrival of Earthquakes 

3 s MCDULAR EARTHQUAKE SOURCE AND 

SENSITIVITY ANALYSIS 

5.1 s Earthquake Source at sites Volume of 

Revolution of Sector about Xj 53 

3.2 t Effect of arc Length, b, on Return 

Period, T , of Peak Ground Acceleration, y 65 

y 

5 .5 s Effect of arc Length, b, on Return 

Period, T , of Peak Ground Velocity, y 66 

y 

5.4 % Effect of arc length, b, on Return Period 

T , of Peak Ground Displacement, y 67 

y 

3.5 * Effect of Maximum Pocal Depth H^ on the 

Probability of the Peak Ground Acceleration 
Y Exceeding yi m s (5,9j 7^ 


21 

23 

30 

56 

37 

39 



Vlll 


Eig. 

5 .6 

Pig. 

3.7 

Pig. 

3.8 

Pig. 

3.9 

Pig. 

3.10 

Pig. 

3.11 

Fig. 

3.12 

Fig. 

5.13 

CHAPTER 

4 

Fig. 

4.1 

Pig. 

4.2 

Fig. 

4.3 

Fig. 

44 

Pig. 

4.5 


Effect of Jlaximim Eocal Depth on the 

Probahility of the Peak Ground Velocity 

Y Exceeding yi m e (5 j93 71 

Effect of Maxiinum Pocal Depth Hq ok. the 
Probability of the Peak Ground Displacement 

Y Exceeding yi m e (5,9] 72 

Effect of Maximum Eocal Depth on the 

Probability of the Peak Ground Acceleration 

Y Exceeding y| m e (m^, “ ) 73 

Effect of Maximum Pocal Depth on the 

Probability of the Peak Ground Velocity Y 
Exceeding y? m e ) 74 

Effect of Maximum Pocal Depth H on the 

0 

Probability of the Peak Ground Displacement 

Y Exceeding yj m e (m^, “) 75 

Effect of Iifegnitude Cutoff on the Probabi- 
lity of Peak Ground Acceleration Y Exceeding 
y 76 

Effect of Magnitude Cutoff on the Probability 
of the Peak Ground Velocity Y Exceeding y 77 


Effect of Magnitude Cutoff on the Probability 
of the Peak Ground Displacement Y Exceeding y 78 

SEISMIC RISK ANALYSIS OP INDIAN PENINSULA 
A MODEL HOMOGENEOUS IN TECTONIC FEATURE 


Grid Numbers 82 

Coefficient of Variation of y and Magnitude 
Relation 88 

Variability of Local Seismicity within Indian 
Peninsula for Magnitude, M > 5* 91 

Peak Acceleration (cn/sec ) Contours Focal 
Depth Lognormal ly Distributed &Me(5>9] 102 


Peak Velocity (cn/sec ) Cent our s-Logncaanal 
Distribution of Focal Depth & M e (5? 9 ] 


105 



ix 


Pig. 

4 .6 


Peak Displacement (cm) Contours - 
Lognormal Distribution of Focal Depth & 

M e (5,9] 

IC4 

Pig. 

4.7 

s 

Variation of VT with Peak Ground 
Acceleration, Y 

107 

Fig. 

4.8 

S 

Variation of pT with Peak Ground 

Velocity, Y 

108 

Pig . 

4.9 

2 

Variation of PT with Peak Ground 
Displacement, Y 

109 

CHAPTER 

5 


SEISMIC RISK ANALYSIS OP 11© IAN PENINSULA 

A MODEL NOIEOMOGENEOUS IN TECTONIC FEATURES 


Fig. 

5.1 


Seismo-Tectonic Hap of India 

116 

Pig. 

5.2 

g 

Idealized Tectonic Features of Indian 
Peninsula 

117 

Fig. 

5.5 

' 

2 

Peak Acceleration (cny^sec ) Contours- 
Lognormal Distribution of Focal Depth, 

M e (5,9] & Nonhomogene ous in Tectonic 
Features 

124 

Pig. 

54 

% 

Peak Velocity (cny^sec) Cent our s-Lognormal 
Distribution of Focal Depth, M e (5,9 ] 

& Nonhomogeneous in Tectonic Features 

125 

Pig . 

5.5 


Peak Displacement (cm) Contours-Lognoraml 
Distribution of Focal Depth, M e ( 5,9 ] & 
Nonhomogeneous in Tectonic Features 

126 

CHAPTER 

6 

O' 

0 

EFFECT OP SCATTER IN THE ATTENUATION LAWS 

ON SEISMIC RISK ANALYSIS 


Fig. 

6 ,1 

0 

Attenuation Law for Peak Ground Acceleration, 
Y 

129 

Pig. 

6,2 

2 

Attenuation Law for Peak Ground Velocity ,Y 

130 

Pig. 

6 .3 

g 

Envelope of Scatter 

134 

Pig. 

6 .4 

2 

Effect of Model-1 on the Exceedance 
Probability of Peak Acceleration, Y 

139 

Fig. 

6.5 

0 

Effect of Model-1 on the Exceedance Probahi- 



lity of Peak Velocity, Y 


140 



X 


Fig. 6.6 s 

Effect of Model-1 on the Exceedance 
Probability of Peak Displacement, Y 

14-1 

Fig . 6.7 s 

Effect of Model-2 and Gaussian Model on 
the Exceedance Probability of Peah 
Acceleration, Y 

144 

Fig, 6.8 8 

Effect of Model-2 and Gaussian Model on the 
Exceedance Probability of Peak Velocity ,Y 

145 

Fig, 6.9 s 

Effect of Model-2 and Gaussian Model on the 
Exceedance Probability of Peak Displacement,Y I 46 

CHAPTER 7 s 

SEISlilC ZONING MAPS OP INDIA 


Fig. 7.1 8 

Peak Acceleration Contours 

150 

Fig . 7 .2 8 

Zoning Map Based on Acceleration 

152 

Fig. 7.3 8 

Peak Velocity (cm/sec.) Contours 

154 

Fig. 74 8 

Zoning- Map Based on Velocity 

155 

Fig. 7.5 • 

Jeak Displacement (cm) Contours 

156 

Fig, 7 £ 8 

Zoning Map Based on Displacement 

157 

APPENDIX A3 .8.1 

PROBABILITY DISTRIBUTION FUNCTION AND 
PROBABILITY DENSITY FUNCTION OF FOCAL 
DISTANCE, R, FOR VCLUME SOURCE AT SITB 


Fig. A 5 . 8 .I.I 

Earthquake Source at Sites Volume of 

Revolution of Sector About X, 

3 

180 

Fig . A3 .8 .1.2 

r c [ 0 ,a sin v ] 

187 

Fig. A 3 .8 . 1.3 

r e [a sin v, 2 a sin (v/ 2 ) ] 

187 

Fig. A3 .8 .14 

r e [ 2 a sin (v/ 2 ), ] 

187 

Fig. A3 . 8 . 1.5 

X c [ Rq ] 

187 

APPENDIX A3 .8 .2 

PROBABILITY OP GENERALIZED PEAK SEISMIC 
INTENSITY 


Pig. A 3 .8 .2 .1 

Domain of Different [Y > y] ; z e (zt;^,<») 

196 

Fig. A 3 . 8.24 

Domain of Different [Y > y] 5 z e (z^jz^] 

197 



xi 


APEENDn; A5.I SECTION OP A SPHERE (EARTH) BY A 
PLANE (fault) 


Fig* » 

A5.I 

Rotation of Axis for Fault Plane X* OY' 

202 

APPEM)IX A5.2 

PRCBiiBILITY DISTRIBUriON FUNCTION AND 
PROBilBILITY DENS ITT FUNCTION OP PCC.iL 
DIST.ilCE, R, FOR AREA SOURCE 


Fig. 

A5 .2 .1 

Area Sovirce due to Fault with Projection 
of Site 0 at Fault Plane 

206 

APPEMDIX A64.I 

THE PROBABILITY DISTRIBUTION FUNCTION 

OF Z 


Pig. 

A6 .4 .1 .1 

Z < 

216 

Pig. 

A6 .4 .1 .2 

2-, 1 2 < ^2 

216 

Pig. 

A6 .4 .1 .5 

Domain of Z 

216 

APPEND IK A6.4.2 

THE PROB/EILITY DENSITY AND DISTRIBUTION 
FUNCTION OP W 


Pig. 

a 6 4.2.1 

V 1 

0 

222 

Pig . 

a 6 4 .2 .2 

Event [ ^ W £ min (w2»w^) ] 

222 

APPENDIX A6.4.5 

PROBABILITY OP PEAK GENER/iLIZED SEISMIC 
INTENSITY, Y 


Pig. 

A6.4 .5.1 

Event [ Y 1 ] 

227 

Pig. 

A6 4.5.2 

Event [ y^ 1 Y 1 72 ] 

227 

Pig. 

A6 .4 .5 .5 

Domain of Z and W 

22 J 



XI 1 


SYNOPSIS 

STATISTICAL ANALYSIS OP SEISMIC DATA 
AND 

SEISMIC RISK ANALYSIS OP INDIAN PENINSULA 
The Thesis Submitted 

In Partial Pulfilment of the Requirenents 
Por the Degree of 
DOCTOR CP PHILC60PHY 
Ly 

SUSANTA BASU 
to the 

Department of Civil Engineering 
Indian Institute of Technology Kanpur 
June 1977 

Earthquakes constitute a major potential cause of natural disaster 
over a large part of the earth-. During an earthquakie the surface 
of the earth and all earth-based structural systems are subjected 
to dynamic loads . Aseismic design and analysis of such systems, 
therefore, requires a prediction of the occurrence of earthquakes 
and accompanying ground motions at a location during the service life 
of the structure. The history of recorded earthquakes and ground 
motions during earthquakes indicate that the occurrence of earthquakes, 
their size, temporal and spatial characteristics, and the ground 
motion at a site are random in nature. The prediction of the design 
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seisaaic loads during the service life of a structure, therefore, involves 
statistical analysis and estrapolation ofavailable seismic data. Ihe 
paucity of seismic data available in most parts of the world, however, 
puts severe limitations on the use of mere statistical results. For 
meaningful results, it becomes necessary to synthesise the statistical 
results with judgement based on geotectooic features and any other 
related information available, even if it is subjective. 

hue to the random nature of the earthquakes, a probabilistic approach 
provides a rational framework for aseismic design. Engineering interest 
in earthquakes encompasses location of a structure at the 'best' site, 
based on relative seismic, risk, and establishing design loads for the 
structural system to provide the optimum balance between cost and accept- 
able risk to life and property. To a designer this information is 
available through codes containing zoning maps and algorithm for computatiOE 
of design loads. 

Ohe current code of practice for earthqijiake resistant design of 
structures in India is contained in IS; 1895-1975. The code divides 
the Indian peninsula in five zones and gives lateral force coefficients 
for each zone. The seismic zoning map is based on the subjective 
judgement of the seismicity of the different parts of the country on : 

the basis of epicenters of the past earthquakes, isoseismals and geotectoni 
features. The code specifies the design acceleration in each zone 
without accounting for the service life of the structure in the evaluation ! 
of the seismic risk. The code also does not provide information 
regarding other design parameters, such as, peak displacaaent and i 

velocity which govern the design of some structures. 



The obj active of the present work is to improve the available ' 
techniques for seismic zoning based on stochastic models of earthquake 
phencmena, and to develop rational seismic zoning maps of Indian 
peninsula, The Indian peninsula is defined by latitude (6°', 40 °] and 
longitude (66°, 98° ] and the seismic data available for a period of 
55 years (January 1917 to December 1972) is used to determine the 
parameters of the stochastic models. An attempt is made to quantity 
the notion of seismic risk based on T~year generalised intensity. 

Statistical analysis of the past seismic activities of the Indian 
peninsula is carried out in the first part of the investigation. 

Histograms of origin-times, latitude, longitude and focal depth are 
obtained, A stochastic model based on origin-times only is made to 
investigate the structure of mainshock and aftershocks for different 
magnitudes. However the exact location of earthquake focus is largely 
ignored in the model. The model parameters are estimated based on the 
fact that an apriori discrimination between mainshocks and aftershocks 
from the earthquake catalogue is indistinguishable. An attempt is 
made to obtain a semi-Markov model of earthquakes occurrence for magnitude 
greater than five. The focal depths are considered to be spatially 
homogeneous and tomporally stationary. The focal depth data is fitted 
with truncated lognormal and a bimodal distribution by the method of 
maximum likelihood. Tests of goodness of fit are performed in each case, 

The seismic risk analysis of the Indian peninsxila is investigated 
for two different cases. The earthquake process is assumed to be Poisson 
in both the cases. In the first case the seismic risk is determined 
assuming that the earthquake source is homogeneous in tectonic features. 
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Sensiti-vity analysis of seismic risk with respect to various dimensional 
parameters of earthquake source, magnitude cutoff and three assumed 
distributions of focal depth is performed. Based on the sensitivity 
analysis model parameters are chosen to evaluate seismic risk of the 
Indian peninsula, !IIhe resifLts are plotted in a contour form for lOO 
year return period. Method of extrapolation to obtain seismic iu.sk 
for any other return period and specified acceptable risk frcm this map 
is indicated. In the second case, the Indian peninsula is assumed to 
be nonhomogeneous in tectonic features based on available geotectonic 
data, Ihe earthquake source now consists of area sources associated 
with faiflts, and volume sources corresponding to major geological zones. 

The seismic risk is obtained and ccmpared with the risk evaluated in the 
first case. 

The available data shows a large scatter around assumed attenuation 
laws, A true evaluation of seismic risk must incorporate the effect of 
this xmcertainity. Two new models are proposed to incorporate the 
effect of the scatter and compared with the deterministic and G-aussian 
model due to Esteva, Steps to incorporate the effect of scatter in 
the seismic risk analysis at a site are indicated. 

Based on the stochastic models discussed in this thesis, the 
available information regarding the geotectonic features of Indian 
peninsula and the current seismic zoning maps of the I.S, Code, 
seismic zoning maps of peak ground acceleration, velocity and displacement 
are prepared for a 100 year return period. The maps reflect the syn- 
thesis of all available information to the best judgement of the author. 
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CHAPTER 1 


INTROJUCTIOH 


1.1 Geneiai 

Soil, air and mter are basic to the creation and sustenance of 
life and civilization on the earth. It is ironical that each one of 
them is also the cause of major nature disaster through earthquakes, 
hurricanes, floods and tsunamis • Of these, earthquakes occur over a 
large part of the earth, without prior warning, and have, therefore, 
caused severe loss to life and property. With increasing population, 
spillover of cities and industiiesiiCo seismically active regions are 
inevitable. The risk to life and economic loss during future earthqu^es 
is, therefore, increasing at an accelerated rate. The engineering 
aspects of earthquake cover: 

i) Minimising damage to existing structures, throu^ strengthening 
and modification, if possible. Demolition of unsafe structures 
to avoid loss of life, 

ii) Ctontingent^r plans to deal with the likely disruption of life-line 
facilities, such as, water-supply, electricity, sewarage after an 
earthquake. 

iii) location, design and construction of new structures to withstand 
future earthquakes. 

An earthquake is caused by Ihe release of energy below the surfece 
of the earth. The energy is propagated in the form of the wave and causes 
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the shaking of the soil deposits, reservoirs and all manmade structtcres 
in contact with the ground* The shaking causes dynamic loads vtiicii 
must be resisted to avoid damage during an earthquake, Aseismic desigp. 
and analysis of any system, therefore, require a prediction of the 
occurrence of earthqviakes and the nature of accompanying ground motions 
at a site during the service life of the system. 

Earthquakes must have occured since the earth came to acquire 
its present form. However, the documentation and scientific study 
of eenrthquakes called seismology, is relatively recent. The history 
of recorded earthquakes indicates that the occurrence of earthquakes, 
their size, temporal and spatial characteristics and the ground motion 
at a site axe random in nature. The prediction of the design stismic 
loads during the service life of a structure, therefore, involves 
construction of stochastic models, based on available seismic data, and 
its extrapolation. The paucity of seismic data available in most 
parts of the world, however, puts severe limitations on the reliability 
of mere statistical predictions, Por meaningftiL results, it becomes 
necessary to i^rnthesise the statistical results with judgment based on 
geotectonic features and any other related information. 

Due to the random nature of the earthquake phenomena, a probabilistic! 
approadr provides a rational framework for aseismic design. Engineering 
interest in earthquakes encompasses location of a structure at the 'best* 
site and adopting desigi loads to provide optimum balance between cost I 
and accept^le risk to life and property. The codes containiag s»ning 
maps and algorithm for the camputation of design loads are expected to 



3 


provide this information. There are several ways by which the zoning maps 
may be prepared, Houaaer and Jenings^^^ have classified the maps as 
(a) seismicity maps, (b) fault maps, seismotectonic maps and seianic 
p3X)b ability maps, and (c) engineering maps. Seismicity maps provide 
subjective information regarding shaking in terms of Modified Mercalli (mm) 
Intensity Scale or other similar scales. Maps showing active fatilts are 
called fault maps, Seismotectonic maps axe augmented fault maps with 
geological information (e.g, local geology, tectonic process etc. X Seismic 
probability or seismic risk maps assign probability of occurrence of 
earthquakes of different magnitudes at each fault and provide information 
regarding areal distribution of intensity. Each of these maps is 
largely subjective in predicting the gromd shaking. Engineering maps 
provide quantitative information regarding ground motion which can be 
used to fLx the seismic load on a structure. Preparation of accurate 
engineering maps is, therefore, of utmost importance, 

1,2 Engineering Aspects of Earthquakes 

1,2.1 Causes and Mechanism. 

She earthquakes of engineering interest are primarily associated 
with tectonic activity in the earth’s crust. The process leading to an 
earthquake is essentially one of the transformations of the vast store 
of energy in the earth's interior. According to Belousov, as reported 

(2) 

by Medvedev'’ \ the vertical oscillatory movement of^tfae jearth’s cr ust 
is the basic factor in the tectonic process of earthquake occurrence. 

The earthquake meohanisa. consists of disruption of continuity in earth's 
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crust composed of more rigid blocks adjoining less rigid blocks. Relative 

movements of these blocks develope high shearing strain at the junction 

of the blocks causing earthq.uakes. These junction zones are knom as 

fatats. The medaanical rupture at a sin^e point within a fault leads 

to a redistribution of shearing strain over the entire fault vihich causes 

possibility of ocouirence of earthquakes in the fault or an adjoining 

fault, A second theory due to Evison^'^’ describes the earthq\iake 

process as the effect of phase change in rocks accompanied by a volume 

change in relatively small volumes of the crust. The validity of the 

above theories cannot be established conclusively due to lack of seismic 

data. The generally accepted theory of earthqT;©ke occ\mence is the 

slip along geological feults. However, these faults cannot extend 

beyond a certain depth owing to hi^ confining pressure and temperature 

in the mantle. The evidence of the occurrence of earthquake at a depth 

(5 ) 

between 600 to 800 kilometers, as catalogued by Gutenberg and Richter , 
stands against the fault slip mechanism. On the other hand, the study 
of Southern California earthquakes by Allen et, al^^^ strongly supports 
the view that most of these motions originate as slips along geological 
faults • It is now recognised that more than one mechanism may be i 

(7) 

possible for the generation of tectonic earthquakes , 

Two entirely different approaches, namely, physical and stochastic 
are available to understand 1he mechanism of the origin of earthquakes, 

(s) 

A physical model, proposed by Bunidge and Knopoff^ , consisted of 
masses resting on a lev^ iy.iface with covCLomb friction and 
oonnected by linear springs and dashpots. This model yielded insists 
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into the eaithq.ugfce mechanisias regarding foreshocks, aftersohocks and 

fault creep, A second approach due to Yere-Jones^^* treated 

earthqviake mechanism as a stochastic process analogous to a leaking 

reservoir. An essentially CJompound Poisson Process model was proposed 

by Shiein and iDoks'oz'' ^ where mainshocks were assmed to trigger the 

( ^ a\ 

earthquake mechanism, later Bavkes' * * introduced a class of 

statioaa3:y marked point process in which the value taken by marks 
(e,g, the energy released by earthqxiakes) influenced the future intensity 
of the process. Althou^ mechanical and stochastic models are based on 
apparently different analogies, they provide a realistic mathematical 
treatment of the behaviour of earthquake mechanisms. 

The in^jortance of the study of earthquake mechanism, was evident 

from the study of Eascon and Cornell'' ^ showing that strong motion 

characteristics near the focus of earthqxiakes were dependent on the 

( l6 ) 

conditions leading to earthquakes, Housner'' ' with some assumptions 

regarding bedrock mechanical properties concluded that the peak ground 

acceleration was 50fo of gravity, aLmilar conclusion could be made 

C 1T 

regarding peak ground velocity^ * , 

1,2.2 Earthquake Parameters and Attenuation law 

An earthquake is characterized by its location (latitude, 
longitude, depth of focus), its time of origin and the amount of 
^ergy released by it. The latter Is a quantita'tlve measure of the 
size of the eairthquskes. 
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( 19 20 ) 

Bichter'a'* * ' magaitu(te, M, is intended to rate an earthquake 

independent of the place of occurrence, The magnitude i/vas originally 
defined as the common logarithm of maximum amplitude of surface waves 
in microns, measured at a distance of 100 Kms. on a Wood Anderson 
Seisomograph, !I3ae empirical relation between energy, E, released and 
magnitude, M, 

log E = i,5M + 11,4 (l* 1 ) 

vdiere E is in ergs. The magnitude, M, characterizes the relative 
strength of earthquake at focus. 

The intensity, I, of earthqu^se is a subjective measure of 
shaking at a point on the earth's surface. This is usually reported 
on a non instrumental ordinal scale. Most commonly used scale is 
modified Mexcalli Scale. Empirical relation between magnitude, M, and 
intensity, I^, at epicentre for different focal depth, h (in kms, ), 

^g(5, 21,22) 

1.5M - 3.5 log h + 3 ( 1 . 2 ) 

Eelationship between the magnitude, M, the intensity, I, and the 

/ n’Z 

focal distance, R (in kms,), has been proposed for different regions'^ 

(24) 

The relation given by Bsteva and Rosenblueth is 

I = 8.16 + 1.45M - 2,46 log E (l.3) 

To quantify the historical data vhich is available in the form 
of the des<sription of damage during past earthquakes, an attempt has 
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been made to correlate the intensity with peak ground acceleration or 
velocity* !Bhe correlation proposed by Esteva and Bosenblueth^^^^, 

(25) 

and Eosenblueth between MM intensity, I, and peak gromd velocity, 
V, is 

I = log (l4-v)/log 2 (1.4) 

(23 ) 

Several authors have proposed correlation between peak ground 
motion (displacement, velocity, acceleration), Y, with magnitude, M, and 
focal distance, R, !Chese relations are also called, attenuation laws, 
and can be expressed in the following general Ibrm 

Y = esp [ CgM - tn E ] (1.5) 

where Y is in c, &s, units, E is in kms. C^, 0^ and are empirical 
constants* 

1.3 Seismic lata 

The historical data regarding past earthq,tiakes are available in 
the following forms* 

(i) Narration of tiie human reaction and effect on structures, 

(ii) Isoseismal maps or intensity infoimation on MM or similar scale, 

(iii) Instxximental data in the form of time of occurrence, magnitude, 
latitude, lon^tude and depth of focus. 

(iv) Ittstaninental data giving gromd motion records and attenuation 
of earthquake parameters, such as, ground acceleration, velocity 
or displacement with distance. 
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!Dhe -iJistrumejaial .data. have become available recently. However, this 
data is not complete for many events, Eor a quantitative treatment of 
risk due to earthquakes, data under serial (i) and (ii) can be converted 
into quantitative measures through relations given in the previous 
sections. One of the major problems in the quantitative treatment of 
earthquake rids is the paucity of data at the regional level, 

1.4 Geology 

A realistic evaluation of ^ismic risk requires information regarding 
geological features in particular the structure of geological formation. 
Links between seismicity and geological formations can be estimated by 
compaiirLag the geological evidence with seismometric materials. However, 

(2) 

these links vary widely \ An important general feature is the fact 
that the stronger the differences in the rate of vertical movement within 
a certain zone, the stronger is the seismic activity. However, zones 
of different seismic inten^ties cannot be demarcated by seismometric 
data owing to its validity for short period in ccmparison with geological 
times, 

!Che geological information is given in the map of tectonic features 
giving various folds, faults, thioist etc, and the lithologic group map 
showing the various type of rock groups in diffei?eat parts of the country. 

1.5 Seismic Zoning 

Seismic zoning maps divide a coimtry or a region in zones with a 
Specified intensity of ground shaking due to earthquakes. The intensity 
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of ^akiag may be specified in teims of MM intensity. Peak ground 
displacement, velocity or acceleration or similar site-based parameters. 
Seismic zoning maps represent a synthesis of hystorical data, geological 
features, attenuation laws and an acceptable balance between cost and 
risk of damage due to earthquakes. Depending upon the treatment of 
uncertainty associated with the occurrence of earthquakes and attenuation 
laws, the available zoning methods may be classified in two groups 

1 * Deterministric 
2, Probabilistic 

1.5.1 Deterministic Approach 

luring the last foirr decades considerable efforts have been made 
to prepare and update seismic zoning maps in all countries lying in the 
seismically active regions. Although the ultimate purpose of seismic 
zoning maps is the same, there exist differences in the approach used 
by different c 0 unt 3 l.es. Early maps were prepared on a deterministic 
approach. In Soviet Union, the first zoning map was prepared in 1937 
by the seismological Institute (Academy of Sciences of the U, S. S. R. ). 

The procedure adopted to prepare the maps consisted of two steps. The 
first step was the identification of various zones on the basis of 
seismometric (magnitude or energy) and geological data. The second step 

( 2 ) 

consisted of forecasting intensity and joining them by isoseismals . 

The procedure adopted in C3aina was understood to be same as in U. S. S.R, 

( 26 ) 

The basis of preparing zoning maps in Japan was seismicity map^ \ . The 
development of seismic zoning in U, S.A. was closely related to bxjilding codes 
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Each coimty /state has its ov?n btiilding provision. Uniform Building 

(27) 

Code of 1967 ’ provided a seianic probability map of United States 

of imerica. In 1959 » Bichter^ ^ produced a zoning map using Soviet 

(2^ 

concept of zoning. According to Medvedev^ ' the map was an over- 
estimation* The woife on the seismic zoning map of U. S. A. based on 
seismic risk an alysis is currently being done by the Applied Technology 
Council, 

The first seismic zoning map of India indicating three earthquake 
zones, demarcating on qualitative basis, as heavy, moderate and li^t 
was produced by Geological Survey of India in. 1934. Jai Krishna^ in 
1959 published a four zme map in terms of likely peak ground acceleration 
in different zones. In 1962, Indian Standard Institution introduced a 
seismic zone map in the code, namely, Indian Standard Recommendation for 

(31^ 

Earthquake Resistant Design of Structure IS: 1893-1962^ The map was 

produced by plotting all known epicenters with magnitude five and above, 

drawing assumed circular isoseismal using the average intensity-aagnitude- 

distance relaticnship end modifying marginally with regional tectonic 

features. Seven seismic zones were demarcated in this map according 

to MM intensity, A revision of previous map was made in 1966 (IS; 1893- 

1966)^^^^, The zone I of IS: 1893-1962 was extended to cover only the 

’taarginal depression” in accordance with Tectonic Map of Geological Survey 

of India, Ihe zones at the fopt of Himalaya were made parallel to its 

teotonio trend, Itevision of shape of location and grade of isoseismal 

around Delhi, Eashtnir, Assam, Uagaland, Andaman, Manipur, Tripura and 

( 33 ) 

Satpiara were made. Seismic zone map IS: 1893-1970^ was produced in 
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accordance with Oldham's hypothesis of the faifLted coast line along 

West Ooast owing to the occurrence of Koyna Earthquake of 1967. 

Different seismic zones were demarcated on the probable distribution of 

the general int^sity of future shocks as guided by the specific 

tectonic features and past seismic history. In the second revision 

of IS: 1893 seismic zones were reduced to five by combining zone 0 

and zone I as zone I and by combining zone V and VI as zone V. Eour 

separate islands were added in Kashmir, H.P. , U.P, and near Nepal, 

Andaman was entirely xinder zone V, A zone IV in the Koyna region 

surrounded by a zone III all along West coast was introduced. Boundary of zone 

II and zone III were redemarcated raising the seismicity in parts 

(34') 

of Gangetic plane. Latest version of IS; 1893-1975 contained the 
seismic zone map of IS: 1893-1970. The earthquake likely to occur in 
the five zones were of magnitude less than, 5, 5 to 6, 6 to 6,5, 6.5 to 7 
and greater than 7 respectively in Bicter's scale, Kaila et, al, * ' 

produced a seismicity map of India which agreed quite well with other 
maps prepared by different methods, 

1,5.2 Probabilistic Approach 

A large measure of uncertainty associated with the occurrence 
of eartiiquafces, their sizes, locations and attenuation laws makes it 
necessary to use the probabilistic approach to develop a rational basis 

(7) 

for earthquake resistant desiga. To. quote Nevanark and Bosenblueth 

"In the past the orthodox viewpoint maintained that the objective 
of design was to prevent failure; it idealized variables as deterministic. 
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fEhis simple approadi is still fruitful when applied to design under 
mild uncertainty, and in situations in which the possibility of failure 
may be contemplated at such a distant future as to be almost irrelevant; 
but vshen confronted with the effects of earthquakes, this orthodox view- 
point seems so naive as to be sterile". 

While the need for a probabilistic approach is generally recognised 
but the difdK-culty lies in the scarcity of data. During the last decade 
intensive efforts have been made all over the world in the application 
•f probabilistic method. Assessment of seismic risk forms the basis of 

(37 ) 

probabilistic approach , A seismologist defines a seismic risk as 
the probability of the occurrence of an earthquake of a specified size 
or more in a region during a specified interval of time, lo an engineer 
it is more meaningful to define seianic risk as the probability of 
occurrence of generalized intensity (e.g. peak ground acceleration, 
velocity, displacement etc, ) greater than a specified value at a site 
in a specified interval of time (service life). To distinguish the two 
definitions, henceforth the former is termed as geological The 

seismic risk at a site not only calls for the assessment of geological 
liak of an earthquake source located at a distance from a site but also 
accounts for its effect by the attenuation of ground motion Eq, (l,5). 

It is generally assumed in both types of analyses that the 
earthquake arrival process is Toisson and the magnitudes are exponen- 
tially distributed. These two assumptiorslead to Gumbel's''^ ^ 

O^e I extreme value distribution for maximum earthquake magnitudes and 
generalized intensities. 



13 


Choosing a critical size of earthqtiakes, Lomnitz^^^^ produced 
a geological ri^ map of Chile in 1964, for a time period less than 
the time period of the availahle seismic data. He subdivided the 
entire region by grid of m x n points, vhsce m is the past epicenters 
of earthquake in the region, and assuming statistical independence of 
origin of critical size earthquake at each grid point obtained geological 
ri^ at each node, 

Iiomnitz and Epstein^^^^ used extreme value theory to obtain 
geological risk of California, The results were in return periods 
of extreme magnitude of earthquakes. However, California region was 
too big to be considered a single unit. 

A major contribution to the seismic rd^ analysis was made 
f 4 l) (42,43) 

independently by Cornell^ ' and Esteva^ ^ , In their work, 

earthquake source geometry was assumed to be either a point, or a line, 
or an area located at constant depth. An attenuation law incorporating 
tte effect of attenuation during wave travel from source to site was 
used to evaluate the peak generalized seismic intensity at the site for 
a Specified retiim period, The unobservable Gaussian error, obtained 
from the linear regression analysis in fitting the attenuation law to 
data, was also incorporated in the analysis of Esteva, He estimated 
the intensity of so\irce by Bayesian statistics to account for regional 
variation in the parameter of magnitude-frequency law. Cornell on the 
other hand opted for classical statistics to obtain the intensity of 
occurrence of earthquake at a source. 




14 


later, Ctorsaell and Vanmarcke^^^ ^ performed sensitivity analysis 
of the parameters involved in the seismic risk model^^^^, IChey arrived 
at the conclusion that major contribution to the risk comes from the 
more frequent smaller events at close projdmity to the site, 

(45 ) 

Goto and Kameda used discrete analogue of Poisson process, 
the Bemouli random variables, to model the arrival process of earth- 
quakes of different intensities. Phis result was a multinomial 
distribution of earthquakes of different intensities during a specified 
time interval, Kie return period of exceedance intensity v?as obtained 
from the above distribution. They further assumed that the gromd 
motion record was a Gaussian process and level crossing was a Poisson 
process in obtaining the probability distribution of maximum grovind 
motion in a single earthquake. Intensity parameters of the model were 
matched with the empirical formula. Several samples of accelerogram were 
simulated for different correlation time^^^^ to obtain probability 
distribution of maximum ground acceleration. Combining asstmied model 
of earthquake process and probability distribution of maximum ground 
acceleration, the probability distribution of maximum ground acceleration 
at a site in a 75 year period was obtained and the expected maximum 
ground- acceleration was computed. Same procedure was used for tte 
ground velocity. The results were presented in contour forms showing 
expected wa yi-nBim ground acceleration and velocity to occur in 75 years 
in Japan, 

Benjamin^^"^ ^ assumed a miiltincmial model for earthquake of 
different intensitie, Bayesian forecasting technique was adopted 
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,(48) 


with diffuse prior^^''^ to optimize desi^ decision at a site. 

id obi 
,(39) 


(49 ) 

lomnitz' incorporated the historical data and obtained a geological 


risk map of Ghilie using his previous technique' 

Mine and Ihvenpoirt^^*^^ obtained a least squS-re solution with 
number of earthquakes at a site that produced an acceleration equal 
to or greater than a specified value at a site per annim. A risk 
map in terms of the return period of an acceleration of 1C^o gravity for 
eastern part of Canada was prepared. They also compared their result 
with that obtained from extreme value theory. 


(5l) 

Algeimissen obtained a geological ri^ map of United States 

of America using distribution of M intensity and strain release data 
and association of strain release data with geological features. 


Relation between fault length and magnitude, the spatial distribution 
of intensity of ground shaking and magnitude, the felt area and magnitude, 
and tbe frequency of occurrence of earthquakes and magnitude were given by 

(52) 

Housner , Sbr an average site in California effective upper bound 
of magnitude was given as a function of seismicity, 

A comparison of geological risk for various models of earthquake 

(53) 

arrival process at a site in Alaska was presented by Chou et, al, , 
Inter-arrivals of earthquakes were considered to be exponential distri- 
bution and ^*'eibull distribution. Both classical and Bayesian estimation 

procedure were used for comparison of gedcgLcal risk for various 
distributional model. They also obtained geological risk Of Alaska 
site using Benjamin approach. 
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(54) 

In 1971, Cornell incorporated the Gaussian errar term in the 
attenuation law in his seianic risk model discussed earlier. He 
also introduced a random aspect in this paper for the structiaral response. 
Assuming damage measure a monotonic increasing function of intensity of 
ground motion the probability distribution of cummulative damage of 
Palmgren -Miner type was obtained, ELnally the probability distribution of 

c cmmulative damage of structure was obtained approximately using approjcuuate 
probability distributions, ejdireme value distribution of Type I, of 
peak structxiral response and assuming coefficient of variation of peak 
response given that the site intensity is small, 

(55) 

Algermissen and Perkins following somewhat similar technique 
(41 ) 

proposed by Cornell obtained a seismic risk map for Utah and 
Arizona, Hiey prepared maps for the region on the basis of peak ground 
acceleration of 50— year return, period and 50-year maximum acceleration 
with 10^ exceedance probability, 

Merz and Comell^^^^ developed a seismic risk model using a 

(57 ) 

quadratic-magnitude-frequency law based on observed data to provide 

a better fit for large magnitude, later, Merz and Cornell^^^^ used 

a trigger model in seismic ri^ analysis to account for aftershocks. 

Spatial distribution of aftershocks were assumed according to a 

distribution obtained by Utsu^^^*^^^ and temporal characteristics in 

(61} 

accordance with Omari’s' ^ law. Based on this investigation it was 
concluded that equivalent event model (no discrimination between shocks) 
provided a conservative estimate of the seismic risk. 
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( 62 ) 

In 1975, Liu and m.gel'' provided geological rii^ of United 
States using available methods, !l?he ri^ was presented in terms of 
magnitude and return period. 

Ihe essential feat^ire of the probabilistic approach to seismic 
zoning is the synthesis of the stochastic model of the earthquake 
phenomena incorporating the temporal and spatial uncertainties, a 
distribution function for magnitude, attenuation laws and geotectonic 
features. In view of the paucity of instrument data, Bayes* theorau 
is used for parameter estimation to make the best use of historical 
and future data. In the invest! ^tion known to author the random 

(7 ) 

nature of focal depth has not been incorporated. There is also 
scope for improvement in the available methods to accomt for the 
scatter in attenuation laws and extrapolation of specified intensities 
at a site for different return periods, cutoff magnitude and acceptable 
risk. The goodness-of-fit of assumed distributions is rarely reported. 
The approach adopted to prepare the c\irrent seismic zoning map for the 

(54) 

Indian peninsula e2j>licitly states, "The Sectional Ctommittee has 
appreciated that there cannot be any statistical approach to the 
problem of earthquake intensity and an entirely scientific basis for 
zoning is also not possible in view of the scanty data available", 

Mae map has been prepared for ground acceleratioi only and does not 
account for the return period. 

1,6 Present Woik 

The present woik is an attonpt to improve the available technique 
for seismic zoning based on stochastic models of earthqviake phenomena, 
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and to develop rational seismic zoning map of Indian peninsula, (Hie 
Indian peninsula is defined by latitude (6°, 40°] and longitude 
( 66 % 98 ° ] . 

Chapter 2 presents an analysis of the seismic data of the Indian 
peninsula for the period of 55 year (January 1917 ~ December 1972). 
ifi.stogram for spatial location of earthquake is obtained first, A 
contour map showing energy released by earthquakes of magnitude above 
five is presented, Ihequency-magnitude relation for Indian peninsula 
is obtained, A stochastic model ignoring exact location of earthquake 
focus and based only on origintimes is constructed to investigate the 
structure of mainshodcs and aftershocks for earthquakes above different 
magnitude threshold, a Semi-Markov model incorporating spatial 
(epicentral) and temporal location of earthquakes is fitted to the 
past seismic data of the region, Ihe focal depth cteta is asstmed to 
be spatially and temporally homogeneous and stationary, She focal 
depth data is fitted to the truncated lognormal and bimodal distributions, 
!Ihe goodness-of-fit is discussed in each case. 

Seismic risk formulation for a modular volume source centered at a 
site is presented in chapter 5. Sensitivity analysis with respect to various 
geometric parameters of the earthquake source, magnitude cutoff and three 
assumed focal depth distributions is performed, 

Ihe Indian peninsula is divided into a 2° x 2° grid and seismio 
risk is evaluated assuming homogeneous in tectonic features at each 
grid point for a lOO year return period based on modular volume source. 
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Hiese results are presented in Chapter 4. The geometric and other 
parameters are diosen on the basis of the sensitivity analysis in 
Chapter 3. Regional variation of magnitude-frequency law is taken 
into account by Bayesian procedure. The resiXLts are plotted in 
contour forms, Eictrapolation and intea?polation procedure for different 
retiOTi periods as well as procediore to obtain T— year maximum 
generalized intensity with exceedance probability a is discussed. 

Chapter 5 gives a brief description of the structural geology 
of Indian peninsula and their idealisation for the study. Seismic 
risk analysis of Indian peninsula is carried out assuming it to be 
nonhomo geneous in tectonic features. The earthquake source consists 
of area source due to known faults and volume sources based on 
geological character. Contour maps obtained from this analysis are 
compared with that of Idae homogeneous analysis prepared in the previous 
chapter. 

Three different methods of incorporating the scatter in the 
attenuation laws are discussed in Chapter 6, Examples are given to 
show the method of accounting the scatter in the contour maps. 

Based on the results of the present investigation, available 
geotectonic information and the current I. S, Code, seismic zoning 
maps are prepared for a 100-year return period in terms of peak ground 
aoceleication, velocity and displacement. These are presented in 
Chapter 7« Thesis is concluded with a summary of results and a 
discussion, of the scope of future research in Chapter 8. 



CHAPTER 2 


STATISTICAL AWALYSIS OP SEISIIIC DATA OP IHDIAH PEHIHSIILA 


2.1 Introduction 

The earthq_ua3s:e process is essentially a stochastic field process. 
However, at present it is not possible to carry out complete and successful 
investigation assuming it as a stochastic field process, primarily, 
because of the insufficient nature of available seismic data. Never- 
theless the statistical analysis of the available data has its in^ortance 
towards proper understanding of the mechanism of earthq.uake in a region. 

It is with this aim that the statistical analysis of seismic data of 
Indian peninsula (latitude ( 6 *,40° ] and longit-ude (66®, 98° ] ) has been 
carried out and presented in this chapter. 

The data analysis is primarily aimed at to investigate the earth- 
quake characteristics with respect to space, time and its magnitude. 

To this erd distribution of focal depth and magnitude frequency relation 
are established in this chapter. Further, some stochastic models on 
interarrival times are developed and parameters of these models alPet 
estimated. The goodness of fit criteria is determined in each case and 
is discussed. 

2.2 Data and General Observation 

The seismic data are obtained from Indian Meteriological Department, 
New Delhi. The catalogue of earthquakes contains data from the year 895 • 
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But upto the year 1917 "the data are very sparse and cannot meaningfully 
he used in the statistical analysis. It is for this reason that 
earthquake data from 1917 onwards have been used in this study. The 
early part of the catalogue contains earthquake data with magnitude 
above 5 only. However with the availability of sophisticated instruments 
during the last decade the catalogue contains earthquake data with 
magnitude less than five from 1965 onwards. In the present analysis 
except for' magnitude-frequency relation the data of earthquake with 
magnitude above five are used . The focal depth data of all the earthquakes 
are not mentioned in the catalogue . Further the catalogue contains a 
large number of earthquakes with focal depth exactly equal to 35 kms. 

It is suspected that 53*0 kms. is attributed to all those cases T/here 
focal depth has. not been recorded properly or some doubt exists on the 
value of focal depth. Hence, all data with focal depth of 35 *0 kms have 
been rejected. It also appears that magnitude data of some of the early 
earthquakes of the catalogue might have been derived from intensities. 
However, this indication is available only for a few cases. In some 
magnitude data only range is specified. The mid point of the specified 
range is taken as the actual data. 

Pig . (2 .1 ) shows the yearly number of earthquakes with magnitude 
above five in the region for 55 years. The mean rate is 25 earthquakes/yr 
A linear regression L = a+bt is performed. The parameter L represents 
earthquakes/yr and t is in years. The results of the analysis are 

a =11 .57 earthquakes/ yr 

b = 0.488 earthquakes/yr ,/yr 
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Multiple correlation coefficient square = O.16O5 

The null hypothesis h = 0 is tested by F-test , The regression 

coefficient b obtained is significant because F =18.05 > F 

O 0 •yyj 

(1,53) ^'12.2, The result implies that on an average about 26 more 
earthquaikes/yr of magnitude above five were recorded at the end of 
period than at the begining. 

Fig., (2 .2 ) gives the histrogram of latitude and longitude data 
at 2° interval. It can be seen that maximum earthquakes during the 
period have occurred between longitude of (70°, 72° ] and between 
latitude of (36°, 38 °] . 

The Table 2.1 provides length of the period; the number of sample, 
the sample mean dispersion and the coefficient of variation, for inter- 
arrival time, magnitude, longitude, latitude and focal depth data. 

Table 2.2 gives the sample serial correlation coefficient of interarrival 
time, magnitude, longitude and latitude of lag upto 10. None of the 
series is completely random since values of the serial correlation 
coefficient of lag one is far above the standard error^ , a ^0.03. 

The correlation matrix of longitudes, latitudes, magnitudes and 
interarrival times are given in Table 2.3. The correlation between 
latitude and longitude is significant and it confirms the presence of 
Himalayan belt and indicates that on an average epicenters of earthquake 
have high latitude and small longitude or viceversa. However, other 
correlations are vary small. 
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Table 2. 1 

Properties* of Skrthquakes with Magnitude above live 



Interarrival 

times 

(days) 

Magnitude 

(log(erg)) 

Latitude 

(degrees) 

Longitude 

(degrees) 

Eocal 
depth, 
(bas, ) 

length of mta 

20064. 95** 





Uos, of Samples 

1374 

1374 

1374 

1374 

624 

Mean 

14.603 

5.743 

29.597 

82.635 

98.861 

ELspersion 

455.979 

0.234 

70.772 

120.954 

C:z__ 

5829.710 

Coefficient of 
variation 

1.463 

0,082 

0.284 

1.333 

0,773 




Origin is fixed at the first event after 1917 

” (1917-1972) lass time to the 

The number is high due to widely separated data. 


lag 


Table 2,2 

Estimated Serial Autocorrelation Coefficients 
1 2 


3 4 5 6 7 8 

int^va. 0.234 0.18 , 0.185 0.170 0.281 0.26 , 0.209 o.aoo 0.^79 o’m 

2^tuha 0.191 0.^1 0.198 0.214 0.224 0.2,3 0.202 0.202 0.219 0.190 

0.147 0.183 0.134 0.066 0.073 0.078 0.060 0.098 0.048 0.097 

0.235 0.224 0.184 0.184 0.162 0.108 0,110 0,114 0.1410.122 


longitude 


Interariival 

times 

Magnitude 

latitude 

Longitude 


Table 2.5 

Estimated Correlation Matrix 

1.0 


0.0396 

0.0015 

-0.0265 


1.0 

-0.0429 

0.1038 


1.0 

-0.6192 


1.0 
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2 .3 Spatial Pattern 

In this section spatial pattern of earthq.uakes is studied • Only 
position of epicenter is considered. Focal position is not considered 
owing to the nonavailability of focal depth data at all origin times. 

In the first part location of epicenter is tested for two different 
time periods (from 1917 to 1947 and from 1917 to 1972) to investigate 
stationarity of spatial pattern. In the later part the data are tested 
for space-time correlation. 

The spatial pattern is analysed by the method suggested by Greig- 
Smith^^^^, The method uses contiguous quadrats (small areas ) covering 
the entire region. The number in each quadrats is counted * M jacent 

pair of quadrats are then combined to give quadrats of double area. 

This procedure is followed until the quadrat itself becomes the region. 

The total sum of squares about the mean of the smallest quadrat or single 
unit block is apportioned as in the analysis of variance. The results 
of the analysis provide mean squares for different block size and are 
reported in Table 2.4. However, choice of period (1917-1947) is arbitrary. 
According to Greig-Smith the mean square value remains small for block 
size being small in relation to the mosaic patches of the pattern. 

However, the mean square increases with increase in block size till area 
of block becomes equal to the area of patch. The mean square value 
remains in this high level with further increase in block size. This 

indicates that the patches themselves are random or aggregated. The 

fall of mean square value from this high level ensures regular arrangement 


of patches 



Table 2. 4 


Mean squares for various block sizes in G-reig-Saith 
pattern analysis of epicentral location 


Period 

1917-1947 

1917-1972 

Block size 

Mean Squares 


1 

0.711 

3.095 

2 

0.641 

3.083 

4 

0.798 

2.728 

8 

0.429 

3.172 

16 

1.392 

9.581 

52 

1.031 

5.743 

64 

2.221 

22.683 

128 

4.073 

30.477 

256 

5.956 

51.090 

512 

6.686 

62.940 

1024 

2.337 

11.124 

2048 

29.907 

242.678 


Number of quadrats is 4096 

Number of earthquakes for the period 1917~1947 is 515 
Number of earthqiiakes for the period 1917-1972 is 1375 
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The result for the period (1917-1947) shews three different mean 
size of patch area while that of period (1917 to 1972) shows two 
different mean size patch area. This indicates that the pattern in 
different time periods are different i.e. the pattern is not stationary. 
For the period (1917 to 1972) the spatial pattern is composed of 
patches of mean size of area "between 256-1 024 units within which patches 
of mean size area "between 8-32 units exist and these patches have 
regular pattern. A regular pattern of patches of three different mean 
size of area "between 256-IO24 units, 8-32 units and 2-8 units constitutes 
the epicentral pattern for the period (l 9 l 7 - 1947 )* 

A second investigation has heen carried out to ascertain virhether 
earthquakes occuring at small epicentral distances from each other have 
also occurred close in time. Existance of space-time clustering will 
imply the S 5 raipathetic behaviour of earthquakes' occurrence. The method^' 
requires a priori specification of closeness of event both in space and 
time, A 2 X 2 contingency table is made by taking every possible pair 
of earthquakes and classifying them in accordance to a priori choice of 
closeness. The a priori closeness specification used in this study are 
spatial distance of 15 O kms. and time interval of 180 days. The 2 2 

contiaency table obtained from data is given "below 


Time in 
days 

Distances 

< 150 

in Kms . 

^150 

Total 

<180 

X = 3329 

23010 

26339 

>180 

38659 

879627 

918286 

Total 

41988 

9082637 

944625 
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The null hypothesis of no spaxse-time correlation, sxiggested 
hy Knox^^^ ^ and later proved by Barton and David by graph 
theoretic approach, is to test that the statistics X is Poisson 
with mean 

X = (26559 x 41988)7944625 = 1170.75 

The null hypothesis is rejected by approximating the Poisson by Normal 

-6 

random variable at P value less than 5^10 . This test confirms 

that earthquake data are composed of main shock and aftershock data. 
However, the test performed has the drawback of a priori choice of 
critical distance and time. The end result of the test is strongly 
dependent on this choice, 

2 ,4 Focal Depth 

In this section focal depth distribution is chosen and estimated. 

The focal depth data contain missing data which are supposed to be that 
of 'normal' depth earthquakes. The variation of focal depth of earth- 
quakes with respect to time cannot be investigated because of these 
gaps in data. Hence focal depth is assumed to be statior^y parameter 

of earthquakes. It is also assumed to be homogeneous in Indian 

^ ' ' - 

peninsula. The maximum focal depth of earthquake in the data is 583 kms. 

Histogram of focal depth is shown in Fig. (2.5). The histogram clearly 
shows a major peak at 5^ kms. and a minor peak at 250 kms. for earth- 
quakes with magnitude above five. The sample properties of focal 
depth are reported in Table 2.1. Two distributions, truncated lognormal 
which is a unimodal distribution and a mixed truncated logncrmal 




31 


distribution which incorporates both the peaks observed in histogram, 
are tried to fit the focal depth data. The truncation level is 
chosen as 600 kms. The subjective judgement in the cutoff is the 
assertion that earthquakes with focal depth greater than 600 kms are 
not likely to cause damage to the structure with foundation in ground. 

The focal depth data is first fitted to a lognormal distribution 
with maximum focal depth, H^= 600 kms. The focal depth data is a 
sample from 

fg.(x) = exp [-0.5(iln x-e^ )^/e2] / [/Oire^ x ^{(iln H^-0^ )/t^ }] 

|xe(0,H^] (2.1) 

where 8^ , Sg are parameters of the distribution and $( .) is probability 
density of N(0 , l). The details of the method of maximum-likelihood 
estimation of the parameters 0^ and are given in Appendix A2.4«1 • 

The estimated values are t 

0^ = 4*252 and 6^ = 1 *0761 

and the large sample covariances of the estimated parameter ares 

2.02 X 10“^ 

02 7*38 X 10"^ 5.55x10"^ 

The goodness of fit is tested by Kolmogorcv-Smirncv test. The 
test statistics for two-sided test is D=0.0866 ^^24 0 05~ ’^•^544 • 

Hence the null hypothesis that the sample is from lognormal distribution 
is rejected at significance level, a =0*05. 



52 


A bimodal distribution for focal depth is used to fit the data. 
The distribution is mixed lognormal with probability density function 


fjjCx) = f p^exp [-0.5(jln X- 

[/2^r0^ X f {(AnH^-e^^)/ ;xe(0,H^] (2.2) 


where 


I Pi - 1 

H = 600 kms 
0 


and 6 ^^, 821 ’ are parameters of the distribution and ^ (•) is 

the probability density of 1T(0, 1 ). The details of the method of 
maximum-likelihood estimation are given in Appendix A2 . 4 .1 • The values 
of the estimated parameters are; 

p = 0.8258 

1 

= 5.966 
8^2 = 5.373 

02 ^ = 0.8475 
8^2 = 0.00466 


The large sample covariances of the estimate are; 

p^ 5*227 X 10"4 

8 ^^ 1.678 X lO""^ 2.016 X 10"^ 
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«12 

1 .561 X 10"^ 

4.653 X 10“^ 

3.525 ^ 10"^ 


®21 

1 .772 X 1 0"^ 

2 JoSa X10”^ 

1 461 X 10"^ 

6.838 X 10"^ 

®22 

-5.996 X 10*^ 

-7.023 X10"^ 

-6.107 x-io"^ 

-1.811 X 10 
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The goodness of fit is tested both by Chi-square and Kolmogorov-Smirnov 
test. The 'test statistics for two sided test is given by D = 0,0614 
^624 0 005 ^ *^525 • Hence the null hypothesis can be accepted at 
significance level a = O.OO5. For Chi-square test the number of class, 
k, is chosen according to Mann and Wald^ . The class limits are 
chosen such that the number of theoretical frequency in each class is 
equal to N/k, where N is the number of samples . The number of class 
is 49 for 624 samples and significance level, “ =0,05. The test 
statistics is T=155»87 59-034 = Xq (43) * Hence the null hypothesis 

■that the focal depth distribution is mixed lognormal is rejected at 

a= 0,05. 

Basic weaknesses of Chi-square test are the requirement for large 
sample, heavy dependence upon the choice of number and position of 
intervals which are generally subjective in nature, and tbe possibility 
of very high type II error for some feasible alternative distributions. 

The parameters of focal depth distribution for different truncation 
are estimated for the distribution given byEq, (2.2) and is reported 
in the followings 


H 

0 

A 

®1 

02 

400 

4 .5808 

1 .289 

500 

4.2900 

1429 
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However, none of these satisfies goodness of fit criteria. Graphically 
mixed lognormal distribution (Fig. 2.3) gives the best fit of the data. 

2,5 Magnitude and Energy 

In this section frequency-magnitude relation^^’^ indicating 
correlation between the number of earthquakes with magnitude above a 
level, M, and. the magnitude, M, is established. The data of earthquakes 
below magnitude five are used in this section. In the latter part a 
seismicity map giving contours of average ^arly energy released by the 
earthquake is obtained . 

A linear regression analysis 

Jtn H = a + bM (2.3) 

where 

M = magnitude of earthquakes 

m = /^ earthquakes /day with magnitude above M is performed 
in this section to obtain "magnitude-frequency" relation. The 
analysis of variance table, confidence intervals of the estimate are 
given in Table 2.5. The test of hypothesis indicates that the 
coefficient b is highly significant. The regression line with 9 ^fo 
confidence intervals is shown in Fig. (2.4). 

A contour map of early energy released is shown in Fig. (2.5). 

The energy released, E, by an earthquake is obtained from the empirical 

4.- (5,19,20) 

relation^ ’ ’ 
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Table 2.5 

Analysis of Variance Table for linear Regression Model 
Hn R = a + bM 

lumber of sample = 23 

Source of Variation sum of squares Degree of Root Mean 

freedom square 


Due Regression 

146.609 

1 


12.11 

Due Residual 

1.957 

21 


0.093 

Total 


148.566 

22 




Value 

Standard deviation 


T-value 


a 

6.997 

0.304 


23.00 


b 

-1.903 

0.048 


39.66 



The covariances of the estimate are 
a 0.09256 

b -0,01428 0.002303 

The multiple coeffidLent correlation square = 0.9868 
The null hypothesis b = 0 is rejected at P << O.OOl 

since R = 1572.87 » F cna^tj^l) = 14.6 

The 95^ confidence interval for a and b are 

a = (6.365, 7.630) 
b = (-2.003, -1.803). 


P-Ratio 

1572.87 



mag n it 
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log E = 1 ,5 M +114 (24) 

where M is the magnitude of an earthquake . The energy released 
hy earthquake at 1 ° by 1 ° grid point is obtained . Smoothing of the 

result is done prior contouring over a 5° by 5° window and using 'a cimre 
which gives a weight O.5 at the central square and weight in other 
square in propoortion to the inverse square of their distance from 
the centre. The contour map shows that maximum energy is released 
in Hindukush and Andaman-Nicobar region. Entire Himalayan, Assam, 
Delhi, Kutch and Koyna region also release considerable amount of 
energy. 

2 .6 Interarrival Times 

The purpose of this section is to investigate properties of 
interarrival times of earthquakes. In the labter part of the section 
stochastic models based on origin time are developed to describe 
earthquake occurrences. The parameters of the models are estimated 
and goodness of fit criteria is discussed in each case . 

Sample statistics of interarrival times are given in section 2.2. 

Fig, (2,6) shows average of successive group of 20 interarrival times 

of earthquake. It is evident from the Fig. (2,6) that a quadratic 

( 68 ) 

trend exist. Test for homogenity of variance^^ for the entire 
period with 69 groups yields 

over all average = 14*603 day 

2 

estimated variance =348.857 (day) 



■ ibter arrival times (days) 




40 


and corrected chi-sqTiare test statistics, T = 1452 53 • ilie test 
statistics is greater than '.'5^ point of chi-sq.uare distribution 
with 68 degrees of freedom. This indicates that there is less than 
0,05 chance that these variance estimates are samples from populations 
with the same variance. This test implies either intensity of the 
earthquake arrival process is different for different time intervals 
or there exist gaps in the data. 

Because of the simplicity of Poisson Process and its success in 
describing many natural phenomena, identification of a point process 
data is first tested for Poisson process. A Poisson process is a 
stationary point process with independent increment. The interarrival 
time of Poisson process is exponentially distributed. A rough indication 
of departure of the earthquake data from the Poisson process is indicated 
in Tables (2.1 and 2.2). The coefficient of variation of earthquake 
data is greater than one and the existance of significantly nonzero 
positive serial correlation coefficimt of lag one are the most striking 
features of the results presented in those tables. However, estimated 
mean intensity of earthquake process with magnitude above five is 
l/y= 14«603 day/earthquake , The goodness of fit criteria is tested 
by chi-square test at significance level, ct=0.05. For earthquake 
sample of 1374» the number of classes is obtained as 67 • The intervals 
are chosen so that the theoretical frequency in each interval is equal » 

The test statistics is, T = 873*47 > 84*821 = (65). Hence the 

null hypothesis of exponentially distributed interarrival time is rejected 
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at significance leTel, a = O.O5. The Kolmogorov -Smirnov Test also 
rejects the null hypothesis because the statistics for two sided test 

is D . 0.1515 > i>i5-,4,.05 - 

A two-state SemioJilarkav process^^^"'^*^^ with two types of interval 
is considered next to fit interarrival data. The probability density 
of the interarrival time is hyperexponent ially distributed for this 
process and is given by 
2 

f^(t) = I ' Pi exp (- y^t) (2.5) 

i=1 

where 

IPi = 1 (2.6) 

Eq . (2.5) states that the earthquake process is conposed of 
two types of shock, namely, mainshock ani aftershock but a priori 
discrimination from the catalogue is not possible. Both main shock 
and after shock posses exponential distribution with parameter 
and respectively. The aftershocks are likely to occur within 

a short interval in comparison to mainshocks. This is ensured by the 
condition 

U2 > P-i (2.7) 

The paran^ters of the distribution Eq . (2.5) are estimated by the 
method of maximum-likelihood* The details of the method of estimation 

are given in Appendix A2 .6 .1 . The estimation and test for goodness of 

fit are carried out for earthquakes vfith magnitude above various levels 



42 


arid are reported in Table 2.6, The i^rpothesis that the interarrival 
time is hyperexponentially distributed is rejected for earthquake with 
magnitude above 5*0) 5*5 a'fid 6.0. A peculiarity for earthquakes with 
magnitude above 5»5 is that it contains less mainshocks than the main- 
shocks in earthquakes with magnitude above 5» The mean value of main- 
shock distribution reduces with increase in magnitude. The mean value 
of aftershock distribution reduces with increase in magnitude excluding 
the earthquakes with magnitude above six. The results (Table 2,6) 
indicate that the catalogue of earthquake contain very few aftershocks 
(less than 2C^) . 

Lastly in view of the analysis of epicentral pattern result, the 

Indian peninsula is arbitrarily divided in 16 equal areas and the 

boundaries of these areas are given in Table 2.7. k test for Markov 

(71 ) 

Chain proposed by Kullback et al is made for the transition of 
earthquakes in the areas defined by Table 2.7. The results of the 
test are 

Statistics Degree of freedom 

Two way independence 503,1974 225 

Iferkovity 881 ,4558 5^00 

Two way/one way independence 1384 •655 3825 

These statistics are to be compared with chi-square distribution 
with given degree of freedom. The two way independence test implying 
P£E^jE^]= P , where E^ is the state of the Markov chain, is 

rejected at significance level, a =0,05. However, both Markovity 



i^yperexponential Distribution of Interarrival time 
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Table 2.7 


State of Markov Chain, its Stationary Probabilities and 
Parameters of Holding-Time Distribution 


State 

latitude 

longitude 

X 10^ 

X^x 10 

1 

(6^.,.13.83] 

(66,74 ] 

1.39 

6.209 

2 

(6.00,13.83] 

(74,82] 

0.71 

3.798 

3 

■ (6,00,13.83] 

(82,90 ] 

0.76 

1.619 

4 

(6^00,13.83] 

(90,98] 

87.23 

5.706 

5 

(13.83, 21.94] 

(66,74 ] 

13.02 

20.891 

6 

(13.83, 21.94] 

(74.82] 

2.87 

6.093 

7 

(13.83, 21.94] 

(82,90 ] 

3.03 

3.021 

8 

(13.83, 21.94] 

(90,98 ] 

27.68 

4.372 

9 

(21.94, 30.55] 

(66,74 ] 

63.32 

8.962 

10 

(21.94, 30.55] 

(74,82 ] 

32,62 

6.871 

11 

(21.94, 30.55 ] 

(82,90 ] 

33.67 

5.473 

12 

(21.94, 30.55] 

(90,98 ] 

258.05 

8.377 

13 

(30.55, 40.00] 

(66,74 ] 

309.73 

7.117 

14 

. (30,55, 40.00] 

(74,82 ] 

84,57 

5.057 

15 

(30.55, 40.00] 

(82.90 ] 

30.63 

7.467 

16 

(30.55, 40.00] 

(90,98 ] 

50.72 

6.082 
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and "two way/one way independence are accepted. Hie two way/one way 
independence test says whether rr x r tables are homogeneous. The 
acceptance of Markovity implies that transition between these regions 
can be considered as Markov chain of order one i. e, P [Ej^ j = 

P £Ej^lE^ j , In view of this result a Semi -Markov process is defined 
in which the states are the areas given in Table 2.7 and holding time 
in the states are exponentially distributed, Formally, let p. . 

Xj 

represent the probability of the next transition to state j given 

that it just made a transition to state i and the time between these 

transitions has an exponential distribution which depends only on the 

state it just entered and is given by H^(t). During a tine interval 

t, let lir(t) represents the total nianber of transition, IT.(t) 

represents the number of transition to state j and N.-(t) represents 

J 

the next transition to state j given that it has just entered state 

i»i.e* 

m mm 

ir(t) = I ir.(t) = I I N (t) (2.8) 

j=i i=1 3=1 

where m is the number of states of the embed ed Markov-chain. 
lyke^"^^^ has shown that for a fixed time interval t and when the 
transition probability matrix, [ P ] , and holding time distribution 
function, ^ = 1,...,m are not functionally related, the 

probloa of maximimi-likelihood estimation reduces to three separate 
maximum— likelihood estimation problons. These are (i) maximum- 
likelihood estimation of p. . the transition matrix of Markov-chain; 
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(ii) maxamum-likelihood estimation, based on non random sample size, 
of H^(.) for i=1,,,.,m and i Jj where J is the last state in 
vshioh the p3?ocess has entered before t, and (iii) maxinivan likelihood 
estimation of I^(») of the renewal process. 

The solution to the problem of maxLnium-likelihood estimtion 

(72) 

of the transition matrix [ T ] of the Markov chain is given' ^ by 

p. . = .(t)/N.(t) (2,9) 

The transition probability matrix for the data is estimated using 
Bg., (2*9) and is reported in Table 2,8. 

" let X. . denotes the holding time of jth visit to state i, 

03a e maximum likelihood estimate of the distribution i=1,...,m 

and i 7^ J is 

\(t) 

= N (t)/ ^ . X.. (2.10) 

where is the parameter of the distribution function, Il(.), 

A 

The maximum~likelihood estimation for is 

Bj(t) 

V (2.11) 

Ki— 1 

where 

m 

u-t- I Z 

i=1 3=1 

The parameters of the distribution of holding time are 
estimated fcoa the data and are reported in Table 2,7. 



Transition Probability Matrix of Markov Chain 
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0,002 0,000 0.000 0.096 0,007 0.002 0,000 O.O33 0,038 0.028 O.C45 0.186 O4O8 0.080 0.026 0,(49 
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0.000 0.014 0.000 0.057 0,000 0.000 0.000 0.014 0.029 0.014 0.000 0*286 0.314 0.072 0.014 0.186 
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The stationary distribution^^^^ associated with [P] of 
the Markov-Ghain is determined by 


<II> = <II> [ P ] 


(2.13) 


where <ii> is a row 
Table 2,7. 


vector of i='|,..^,m. and is reported in 


then simply 


The probability density fmction of the interarrival time is 

( 63 ) 


fj(t) = i n^h^(t); t > 

i=1 


0 


(2.14) 


■vdiere h^(.) is the probability density fimction corresponding to 

The goodness of fit criteria is tested first by chi-square test 
at significance level, a = 0.05, with 1374 samples. The number of 
class is 67 and the test statistics is T=80,092 > 48.602 = ^i-(34), 

0,05 

Hence the null hypothesis that interarrival of earthquakes is distributed 
as given by Bi. (2.14) is rejected. Hie Kolmogorov-anirnov test 
yields statistics for both-sided test, I>=0,01^68 > 0.03669 = 

Hence the null hypothesis is rejected. 


1374,. 05* 


Among the models tried to fit the interarrival time data the 
best fit is obtained by B3.. (2. 14). However all these models donot 
satisfy any of the goodness of fit test for earthquakes above magnitude 
live. The saai-Markov model can be used for seismic risk analysis if 
suf£loient data in a region are available. 
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2,7 Sumnary of the Results 

The results of the statistical analysis can be sunmerized as: 

(i) the spatial pattern is not stationary , 

(ii) the earthquake data is space time correlated, 

(iii) the intensity of arjrival process might be different in different 
periods, 

(ir) a long cycle trend mi^t exist in the data and 
(t) among all the models tested the Semi-Markov model provides the 
best fit. However, these conclusions are arrived through first order 
and graphical analysis only. 



CHAPTER 3 


MODIILAE EARTHQUAKE SOURCE ARE SERSITIVITI ARAIYSIS 

3. 1 latroduction 

To construct the seismic zoning maps of the Indian peninsvCLa, 
on a macroscopic scale, the seismic risk at a location is expressed 
in this thesis in two ways; 

1, Generalised intensity for a specified return period. 

2. Generalised intensity corresponding to a specified exceedance 
probability and return period. 

The peninsula is divided in 2 ° x 2° grid points and the 
generalised intensity is determined at each point. The zoning maps 
are prepared by drawing the contours of equal intensity. The require- 
ments for evaluation of generalised seismic intensity at a location are: 

1, Identification of the sources of earthquakes which may be 
idealised into point, area or volume soxirces. 

2, Construction of stochastic model of the earthquake activity in 
the sources incorporating the temporal, spatial and magnitude 
characteristics based on seismic data, 

3, Assumptions regarding the atteniiation laws, 

Tbe identification of earrthquake sources near a location can be 
done on the basis of geotectonic infomation and/or epicentral ^ta of 
past earthquakes. lue to tectonic natxire of earthquakes, epicentral 



data, if available over a long period of time, should esdiibit a strong 
correlation with the geotectonic features. To construct the seismic 
zoning map of the peninsula on this basis, earthquake sources can be 
identified at each grid point and the seismic risk can be evaluated, 

&jch an approach explicitly accounts for the non-homogene ous nat^lre of 
the geotectonic features in a large peninsula like India and is adopted 
in Chapter 5 . 

While explicit incorporation of nonhomogeneotis nature of geotectonic 
feature is desireahle in calculating the seismic risk for preparing the 
zoning maps, there are two aspects connected with it which should be 
considered: 

1 , In many locations, particularly in alluvial areas, it is difficult 
to identify distinct earthquake sources and estimate parameters associated 
with the stochastic model at each source. 


2, Use of zoning maps for design of engineering systems may require 
extrapolation of specified intensities for different return periods 
(service life), exceedance probability, cutoff magnitudes and effect 
of scatter on attenuation laws, Itor maps based on non-homo geneous 
tectonic features such extrapolations will be different for each grid 
point and, therefore, too cumbersome to be incorporated in the codes of 
practice for earthquake resistant desigp,. 


In view of the above conments the viability of assum ing the 
In^an peninsula to be homogeneous in tectonic features is examined, 
A modular volume source centered at each grid ^ ^ 
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standard earthquake source, Eie earthquake occurrance at the soirrce is 
assumed to be spatially homogeneous and temporally stationary, Uri^cr-mj 
truncated lognormal and mixed truncated lognormal distributions are 
assumed for the focal depth. !Ehe mean arri'/ol rate is, howe-ver, 
calo\3lated at each grid point based on regional seismic data which, as 
pointed earlier, is expected to show strong correlation with tectonic 
features. Contour maps based cn homogeneous and non-fiomogeneous 
tectonic feature assumptions are compared in Chapter 5, 

Hie geometry of the modular volume so’jrce, distribution functions 
of magnitude, focal depth,- temporal model of earthquake occurrence and 
assumed attenuation laws are described in tl'is chapter. The influence 
of eadh of these parameters is synthesised to determine the seismic 
risk at the grid point for the modular source. Effects of ^ometric 
parameters of the source and magnitude cutoff are sttdied to establish 
the parameters for the modular source and are used in the next chapter 
for preparing seismic zoning map, 

3.2 Modular Earthquake Soxirce-Geometry 

The modxilar earthquake source, to determine the peak intensity 
at a site, is a 'volume elaaent formed by a sector of revolution about 
the site with an arc, length, b, and maximum focal depth, in which 
future earthqr&kes can potentially occur as shown in Tig. (3. l)« 33ae 
volume element is assxmed to be homogeneous in tectonic features and 
pTOsence of faults, if ary, within it is disregarded, IiOcation of an 
earthquake within the earthquake source is temporally stationary and 
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Spatially homogeneous. Occurrence of earthqtyce is eqt^lly likely 
in the azimuth direction, a , and in the latitude direction, 0 , 

Sig. (3.1). 

3»3 Magnitude-Erequency law 

Hypothesis of magnitude stability, i,e, in any given realization 
of the earthquake process the magnitude distribution is (a) stationaxy 
and (b) indepejadent of rate of occurrence of earthqi^kes^^^**^^^, is 
accepted in this study. This assumption implies that the magnitude of 
earthqviakes is a regional variable and the magnitude distribution is 
represented by a 'space-series’ indepoident of time. Consequently the 
magnitude distribution is assumed to posses a continuous probability 
density function denoted by fj^(, ), . 

In early forties Japanese and American scientist independently 
discovered that the number of earthquakes in a region decrease 
er^nentially with their magnitudes. This relation^ip is known as 

(79) 

"ma@aitude-frequency equation" due to Ishimoto and lida and 

(5) 

Gutenberg and Eidater and expressed as 

log = a- bm (3. 1) 

vdaere is the number of earthqmkes greater than a specified 

threshold magnitude m and a and b are empirical parameters, 
normalization of Bq, (3. 1) and subsequ^t differentiation with respect 
to m yields the probability density ftmction of magnitude, M, as 

fj^(m) = e [ - am ] J (3.2) 


vtiere g = b Jin 10 
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where 

lA = 1 - exp [6(m^-mp ] (3,5) 

3.4 Ibcal Depth 

In the investigation of seismic ri^ analysis known to the author, 
focal depth of earthqv&kes has been treated as a deteiministic parameter. 

It is well known that focal depth has a measvire of uncertainty associated 

(7) 

witii it , The focal depth, H, is assumed to be temporally stationary 
and posses a continuous probability density function fg(, ) with a 
maximum focal depth, 

If past seiSEiic data of focal depth is not available, it is logical 
to assime that earthquake is equally likely to occ\ir upto a depth, 

The probability density of focal depth, H, cctrre spending to this assumption 
is 'uniform* and is given by 

1/H she (0,H^ ] 

fjj(h) = (3.6) 

0 I otherwise 

An analysis of available focal depth data indicates that the 
shallow focais earthquakes predominate over intermediate and deep focus 
earthiuakes in Indian peninsula, Ihe histogram in KLg, (2.3) shows a 
distinctive minor peak at 250 kms. To incorporate the above features, 
a truncated lognormal and mixed trmcated lognormal distribution are 
proposed. The probability dmisity function of these distributions are 
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1* Lognomal distribution 

fg(h) = exp [ -0.5 Un /2^ h ${ (to } ], 

; h e (O, ] (3.7) 

2. Mixed lognormal distribution 

2 

fjj(ii) = I Pj_ ®^P [-0«5(£n h-e^^) /e2i]/[ / 27102 h 

X 1 

HUn Hq-6^^)/ ^ } ] ; h e (0, ( 3 . 8 ) 

wfaere 

2 

I Pi = 1 (3.9) 

i=1 

and *(.) is the probability distribution function of Iif(o, 1 ), 

and { 0 ^^ } i,j = 1 , 2 } are the parameters of the distributions, 

3.5 SDemporal Model of Biirthq\;eike Sburce 


An earthquake pix>oess is essentially a clustered marked point 
process with memory owing to the well accepted hypothesis of elastic 
rebound theory, presence of aftershocks and often presence of swarms 
of earthquakes. However, quite extensi-ve amount of past seismic data 
is required to model earthquake process as a clustered marked point 
process with memory. In this analysis, the earthquake p3?ocess is 
assumed to be Poisson process with an intensity, y , for magnitude 


greater than a threshold mgaitude, m^, Ihis assumption is physically 


inconsistent but is considered adequate for engineering purpose 


(82,83) 

. # 


Statistioal test on past seianic data after removal of swaims and 
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aftershocks from the catalogue do not reject the null hypothesis 


that the earthquake process is Poisson decisively 


( 84 ^ 8 ) 


3.6 Attenuation law 

To include the effect of wave p3?opagation within liie modular 
earthquake source, deterministic attenuation law which estimate the 
peak seismic intensity (e,g, acceleration) at a site as a function of 
earthquake magoitude, M, and focal distance, E, from the site to •the 
location of the hypocenter of earthquake is assumedi The generalized 
attenuation law proposed by Kanai^^^ and Esteva and Eosenbleeth^^^ ^ 
is used and is given by 


Y = exp [ C^M - in (R + C^)] (3»10) 

where Y is the required peak ground intensity, C^j Cg) and 
are empirical cmstants. The unit of E is in kilometers and the 
•values of (c^, C^, C^, C^) are taken as (2000i0, 0*8, 2.0f 25.0) , 
(l6,0, 1,0, 1.7j 25w0) and (7.0, 1.2, 1.6, 25. o) for peak ground 
acceleration (cm/sec^), peek ground velocity (ca^/sec) and peak gromd 
displacement (cm) respecti'vely, 

3,7 General Boimulation of Seimnic Bisk at a Site 

3»7. 1 Probability distribution function of generalized pe^ seismic 
intensity 

In the Eq. (3. 10) let 
Z = G^ exp (c^) 


(5.11) 
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and W = exp [ -0^ M(e + C^) ] 


(3.12) 


The joint p3?obability density function of Z and W, assuming 
statistical independence of focal distance, E, and magnitude, M, of 
earthq.uafce is gi-ven by 


= f^(z) f^(w) 


iff\ 


where 


(C +1)/C_^ -1/C_ 

f^(w) = l/[ C„ w fj^(w ^ - C^) 


} w e[w^, Wg) 

fgCz) = l/(zC2) {(m z - in C^VCg} 


with 


Wi = (E^ + C^) 

“^3 

^2 = 


-C, 


and E = maximum focal distance 
o 


= /a^+ (a - Hq)^- 2a (a - H^) cos v 


(3.13) 


(3.14) 

(3.15) 


(3. 16 ) 

(3.17) 


(3. 18 ) 


The parameters a, and v are the r^ius of earth, maximum focal 
depth and the angle subtended at the cent3:^e of earth by an exfci^e 
point in the modiXLar earthquake source (SLg. 3.l) respectively. The 
probability density ftoction, fg(,), of focal distance, E, is derived 
on the assumed probability density functions of the location variables 
of future earthquakes in the source. The exceedance probability of 



generalized seismic peak intensity, Y, is obtained from Bj. (3. 13) and 
is given by 

Sy(y) = rCY>y]= [ j f^(zjw) dz dw (5.19) 

{(z,w0:2w > y} 

The above result yields the probability that peak seismic intensity, 

Y, at the site will esxceed a certain value, y, given that an earthquake 
with a magnitude in the domain of interest occvirs somewhere in the 
modxjlar earthquake source. 

3.7.2 T— yeeir generalized peak seismic intensity 

To consider the random number of occunrenoes in a given time 
period, it is assumed that the arrival of earthquakes is a Poisson 
process with intensity, y , in the entire earthquake source. The 
number of earthquakes, {lT(t), t ^ 0 } , which occur in the earthqtxake 
source during a time interval, t, is given by 

Pjj(n,t) =? [l(t) = n] = (vtf/nl exp ^yt); n=0, 1,2,3,. .. (3.20) 

Mong all these events, the earthqu^es yjaich produce peak 
generalized seismic intensity at site ^ove a threshold value, y, is 
given by Eq, (3.19), Since the events, [ Y > y ] , arrive independent 
of each other aid the earthquake process as such is a Poisson arrival 
process, it can be shoiiaa^^*^^ that the number of earthqiiakes, {L(t),t ^ O) 
wincii produce peak generalized seismic intensity at a ate greater tiian 
a value, y, during a time interval, t, is given by 

l^(jl,t) = P [L(t) =i]= (Vt Sy(y)>%i exp [ -Vt Sy(y) ] 
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!I3ae probability of no occurrence of the event, [ Y > y j, 
during a time interval, t, is obtained by substituting il = 0 into 
Eg., ( 5 . 21 )* Ihis gives the probability distribution function of 
masdmum peak generalized seismic intensity in a time interval, t, and 
is formally written as 


ly (y;t) = P 1 y;t J = Pj.( 0 ,t) = exp [ -pt S^Cy)] (3*22) 


max 


!Sie annual maximum peak generalized seismic intensity is 
obtained by putting t = 1 in the Eq, ( 3 . 22 ) and is gi-yen by 


Ey (y;i) 

max 


If Y < 
*• max — 


y;i ] 


-1 - vSy(y) 


Pj^(0,l) = exp [ -p SyCy) ] 
(*/ hSy(y) « 1 ) 


Therefore 


]= VS^(y) 


(3.23) 


The average retixm period, T , of peak generalized seisoaic 

V 

intensity exceeding y is defined as 


T =i/p[Y > yji] =i/[vs„(y)] (3.24) 

y max a 

QSie T-year peak generalized seismic intensity at a site is 
obtained by solving Eq, (3.24) for y given a service life, T^. 

The ^neral formulation of the problem of seiaaic risk analysis 
at a site for the modular earthquake source involves two probability 
density functions, namely, the probability density function, %(•)» 
of the magnitude, M, and that of the focal depth, H, fg(, ). These 
must be established on the basis of local seismic data. 
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3»8 Probability IDLstribution. of Peak Generalized 3Bisinic Intensity 
for Assumed listributions of Bbcal Depth, H, and Ma^aitude, M, 
of Earthquakes 

3.8.1 Magnitude, M, cutoff at lower level i, e, me (m^, ") 

The probability of peak ^neralized seiaaie intensity, Y = ZW, 
exceeding a specified value, y, is obtained from Eq, (3,19) for 
magnitude, M, cutoff at lower level, m^, and is given by 

Sy(y) = P [y > y ] = 1 ; y £y^= exp [ Sn(R^+C^)] 

E (3.25) 

o 

Sy(y) = / s(r) fpCr) dr + F^(xp ; y^ i. Y 1 y2 (3.26) 

^1 

R 

o 

Sy(y) = / s(r) fg(r) dr j 7-^ 7^ C^exp: [ 0^] (3,2?) 

THhere 

S(r) = eap [ ${^ 0 + C^- sn(r + C^) - m y J/Cg } } (3.28) 

x^ = eap C An Jin y)/Cj ] - 0^ (5.29) 

!Ehe fjj(. ) 2 ind T'g(. ) are the probability density and prcbab^ity 
distribut icn function of focal distance, R, for different aaaamed 
distribution of focal depth Ejs, (3.6-3. 8), The detailed derivatim 
of these functions are given in Appendix A3.8, i. The details ^ 

Iqs, ( 3 . 25 - 5 . 27 ) are given in Appendix A3, 8, 2. 

3.8.2 Magnitude, M, cutoff at both ends i.e, m e (m^, m^ j 

The probability of the event, [ Y > y] is obtained from 
Eq. ( 3 . 19 ) when magnitude, M, of earthquakes is being cutoff at both 
ends and is given by 
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Sy(y) = P [y > y ] = 1 j y < y^= exp [ C 2 m^-C^)m(E^+ c^) 
Sy(y) = 1 - [ 1 -Pj^(x ^)3 K + K J S(r) fjj(r)dr 


] 

( 3 . 30 ) 


5 1 y ( 3 . 3 l) 

Sy(y) = (i-K) ^’^(xg) +K PjjCx^) + / ^ S(r) fjj(r)dr 

; yj < y ±7^= exp [ c^m^- in c^] (3.32) 

E 

O 

Sy(y) = 1-K-« / s(r) fjj(r)dr ; y ^<7 <. y^ = 


exp [Cgm^- Me^+C^) ] ( 3 . 53 ) 

^2 

Sy(y) = ( 1 -K) ^^(xg) +K / S(r) fjj(r)dr; y>max(y 2 ,y ) ( 3 . 34 ) 


inhere 

Xg = exp [ (C 2 m^+ £n G^- to y)/c^ ] - ( 3 . 35 ) 

One of the Bq, ( 3 . 32 ) and Bq. ( 3 . 33 ) is to he used according 
to the numerical -values of the rel e-want parameters. Details of the 
aho-veare given in Appendix A 3 . 8 , 2 . 

The integrand in Eqs, ( 3 . 25 - 3 . 27 ) and Eqs, ( 3 . 31 - 3 . 34 ) cannot 
be obtained analytically, Numerical integratim by Idle method of • 
Gaussian quadrature is used to obtain the results. 

3.9 Eesults aM Discusstion 

The paarsHBeters of the probabili-ty density function of focal 
depth) H, for -various assumptions and -ttiat of magni-fcude, M, are estimated 
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in Chapter 2. The intensity of the arrival process, y , is estimated 
in this chapter by the ratio of number of past earthquakes in the modular 
souirce at a site and the time period of the seismic data. The radius 
of earth, a, is assumed to be 6371 kms. 

The effect of variation in modular source geometry parameters, 
namely, arc length and maximum focal depth, on seismic risk is inves- 
tigated. 

The sensitivity of seismic risk due to variation in arc length, 
b, is investigated assuming the intensity, y , of the earthquake arrival 
process at a source is directly proposrtional to the volume of the source. 
The return periods for different arc length, b, are calculated at a site 
(latitude 31 longitude 95°) for various level of peak ground 
displacement, velocity and acceleration. The results are daown in 
5igs, (3, 2-3.4) for uniform and lognormal focal depth distribution with 
maxLnm focal depth, H^= 600 kms, and magnitude, m e (5,9] . iigs. 

(3, 2-3,4) indicate that the smaller the arc length the larger is the 
return period (lesser seiaaic risk) for all levels of generalized peak 
seianic intensity. However, the increase in return period due to the 
reduction of arc length, b, reduces as the peak generalized seiaaic 
intensity increases, Erom Eigs,(3. 2-3,4) it is also observed that the 
return period remains constant beyor^ arc length, b = 150 kms, for all 
practical purposes. The seismic risk is less for uniform distribution 
of focal depth in comparison with lognormal focal depth distribution 
with the assxmed maximua focal depth, H^= 600 kms., and a given level of 
peak ground intensity. 




= 50 cm/ 
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The effect of Tariation of maxim-uiri focal depth, on the 
seismic risk is investigated firstly for uniform distribution of 
focal depth assuming constant intensity parameter, u , of the earthquake 
arrival process for various geometric configuration of sources. The 
results for magnitude, m e (5,9] are shown in Tigs. (3. 5-3. 7) and 
that for magnitude, m e (5, ”) are shown in Tigs. (3.8-3,10), It is 
concluded from these figures that for a constant return period peak 
grotuad intensity increases with decrease in maximum focal depth, 

Seismic risk changes marginally for various peak seismic intensities 
(Table 3. l) when intensity, y , of arrival process for various geometry 
of sources having lognormal or mixed lognormal distribution of focal 
depth is constant. 

1 difference in behaviour Table 3.1 and Tigs, (3. 5-3. 7) is observed 
between the cases of uniform and lognormal distribution of focal depth 
when maximum focal depth, is cutoff at different level. Since 
maximum focal depth of past earthquake data of Indian peninsula is 
383 bms. the minimum cutoff of is considered to be 400 Ms. for 
the lognormal distribution of focal depth. The estimation of parameter 
involved in uniform distribution of focal depth does not require the past 
earthquake data and hence does not properly reflect the earthquake 
activity in the source. Moreover this assumption increases the number 
of deep focus earthquakes in comparison with lognormal distribution of 
focal depth and consequently lesser risk is predicted for equal maximum 
focal depth cutoff. The risk for various cutoff of maximim focal 
depth, remain almost constant in the lognormal distribution of 
focal depth. The insignificant variation of risk is essentially the 
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Table 3. 1 

Eatio of return periods for lognormal («,) and mixPd 
lognormal (mfc) distribution of focal depth, H, 

b = 150 kms, 

m e ( 5 , 9 ] 


Eatio of return period is defined as liae ratio of return period 
of any assumed focal depth distribution and lognormal focal depth 
distribution with maximum focal depth, 600 kms. 


Distributioa of focal depth 
maximum focal depth, 

mil 

600 loasi 

1 

500 kms. 

£ 

400 kms. 

peak gromd acceleration 
in cm/sec^ 

Ratio 

of return periods 


25 

0.983 

1*007 

1.065 

50 

0.973 

1*002 

1.078 

100 

0.970 

0.987 

1.084 

200 

0.962 

0.994 

1.097 

500 

0.999 

.0v984 . . - 

1.115 


peak ground velocity 
in cm/sec 


5 

0.996 

1.005 

1.053 

10 

0.999 

1.005 

1*059 

25 

0.975 

1.004 

1.076 - 

5.0 

0.958 

0.998 

1-095 

90 

0.956 

0.998 

1.112 

peak ground displacement 
in cm 

10 

1.005 

1.0C6 

1.038 

20 

1.005 

1.005 

1.041 

30 

1.006 

1.006 

1.044 

50 

1.009 

1.006 

1.051 

100 

0.982 

1.006 

1-069 



y (cm/sec ) 








Uniform focal depth distributi! 
MU ( 5,9] 


EFFECT OF MAXIMUM FOCAL DEPTH, Ho 
ON THE PROBABILITY OF THE PEAK 
GROUND VELOCITY, Y.EXCEEDING y 



Uniform focal depth distribution 













y ( c m) 
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FIG/ 3-13 EFFECT OF MAGNITUDE CUTOFF ON THE PROBABlLiTY OF THE 
' PEAK GROUND DISPLACEMENT , Y, EXCE E DING y 
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reflection of iiie estimated parameter values involved in that distri~ 
bution. Contrary to the sin^jle conjecture that the seismic risk 
for lognormal distribution of focal depth is upper bound, for uniform 
distribution is lower bomd and for mixed lognormal distribution of 
focal depth is inbetween two extreme bounds, it is observed that the 
seismic rislc is maximum for mixed lognormal distribution, This 
observation suggests that in the Indian peninsiSLa a better fit of 
past seismic data in the intermediate focus region increases the 
frequency of occurrence of earthqu^e in the shallow focus region 
also and hence the increment of risk. The modular source with 
uniform focal depth distribution and having niaximum focal depth, 

®o~ kms,, gives compairable seismic risk at a site with a modular 
soirrce of maximum focal depth, 600 kms,, and having a lognormal 
distribution of focal depth. The arc length of Ihe both earthquake 
sources is 150 kms. 

The sensitivity of seismic risk due to magnitude cutoff is 
shorn in jBigs, (3.11-3.13). Jhr constant return period the increase 
in Upper magnitude, m^, cutoff increases the peak ground generalized 
intend ty. The dashed line in Hgs. (5.11-3.13) represents the 

range of exceedance probability of Y in Indian peninsula. The 
effect of magaitude cutoff at upper level on seismic ride is 
insignificant for larger exceedance probability of Y. ihom these 
figures it is seen ttiat the distribution assuaaption of focal depth 
is more pronouuaced for peak ground acceleration ani peak gromd 
velocity than for peak ground displacement. 
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jPram. the results presented in this section it is concluded that 

modular source with arc length, b = 150 kms, and maxMoum focal depth, 

H = 600 kms is to be used for evaluation of seismic risk of Indian 
o 

peninsula, The distribution of focal depth is lognormal and magnitude, 
m e (5,9 j . 



CHAEPER 4 


SEiailC BISK AMIYSIS OP INDIAN PMINSUIA-A MODEL HOMOGENEOUS 

IN TECTONIC PEATURIS 


4. 1 Introduction 

The modular volme source of earthquakes presented in the 
previous chapter is used for the seismic risk analysis of the Indian 
peninsula in this chapter. The peninsula is divided in2° ^ 2° grid 
as shovm in Pig. (4. l). Peak generalised seismic intensity in terms 
of peak gromd acceleration, velocity and displacement are determined 
at each gifLd point. The earthquake source at grid points are assumed 
to be homogeneous in tectonic feattxres i. e. precise information re^rding 
location of faults within the modular earthquake source at grid point 
is either not available or ignored. Every region within Indian 
peninsiola is assumed to be prone to earthquakes. Nonavailability of 
evidence of occurrence of earthquakes in the available data of past 
earthquakes does not imply that the region is free from earthquakes. 

The acceptance of the hypothesis that any region is free from seismic 
activity is synonymous of a dying planet. 

The scarcity of data leads to a situation where classical 
statistics fail to provide meaningfifL results. A Bayesisin approacii 
iS| therefore, used to estimate the intensity parameter, p , of the 
earthqiiake arrival process at a source, !Hie local seismicity i, e. 
the probabilities of occurrence of earthquakes in the Indian peninsiila 
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is obtained in Section 4,2. The regional seianicity in different 
regions vdthin the Indian peninsula is estimted next. The peak 
generalised seismic intensities for lOO-year return period are 
estimated at each grid point and contours of equal intensity are 
drawn. Curves to extrapolate 100-year intensity at a site for a 
different return period, or for other specified exceedance probability 
are presented. 

4.2 local Seismicity 

!Hae local seismicity of Indian peninsxola possessing a volme, 

Yf of earth is assumed to be constant. The intensity parameter, v , 
of the earthquake arrival process for magnitude greater than a threshold 
value, M, is assumed to be a function of magnitude, M, and volume, V, 
and is given by 

S-n p = ln{aY) - BM (4. l) 

vdaere a and B are empirical constants. The relation Bi. (4. l) 
is a representation of magnitude-frequency law. The parameter, a, 
of the magnitude-frequency law (Eq. 3.l) is a function of volume. 

Since the assumed arrLv^ of earthquakes is Poisson process , the 
interarrival of earthquakes is a sample from exponential distribution. 

The probability density of interarrival time is given by 

f (x) = y exp (- yx) ; X > 0 (4,2) 

X 

The parameter, y , will be different for various regions of Indian 
peninsula as stated in the Introdxiction and is estimated by Bayes* theorem. 
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4.2.1 Bayes* Theoraa 

Bayes' Hieorem states that 'posterior probabilities* are equal 
to the 'prior probabilities' and the coirespondlag likelihood -values 
disregarding the normalising constant. IThis formula can, of course, 
be stated in logarithmic form, in which Bayes' Theorem is nlLe of 
addition of information. In the language of Bayesian statistics, 
the model density is 

m^(x) = y exp (- y x) ; x > 0 (4.3) 

The likelihood with respect to y of the observations (interarrival 
times) is 

(y) = exp (m fi-n y - y ^ z^) (4.4) 

The natural conjugate Bayes density may be taken to be 

b(y) = lA(m) exp[m £n( z^) + (m-l) £n y -y ] (4.5) 

vhexe r (,) is a Gamma function. 

The natural conjugate Bayes density with m = z = 0 reduces 
to the density l/ y i, e. uniform on 6 =: £n y , which is a valid 
indifference density but pro-vides neither a stochastic upper bound nor 
a stochastic lower bound on y , 

let the first arri-val occurs at x > 0. The posterior 
density is given by 

bjj^(v) = l/r(m+l) exp [(m-fl) £n(x-4- ^ z^) + m £n y 

- y(x -*-12;^)] (4.6) 
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With m = z = 0 the posterior is proper as long as x > 0, aad 
is given by, x exp (-yx). As x 0, b(y) approaches a unifoim density. 

Let n events have occurred with interarri-val, t The 
posterior density is 

b(y) = l/r(n4m) exp [ (n-4i]i) ^n ( I z. + I t.) + (n+m-l) iln v 

^ t5 

- pil ; y ^0 (4.7) 

The Bayes estimate of y under quadratic loss function is 
y and is given by 

y = E [y ] = (n-Hii)/( I z^+ (4.8) 

and coefficient of variation of y is 

c(y) = l//n^m" (4.9) 

4.2.2 Estimation of Local Seismicity 

In the estimation of local seismicity the historic evidence 
is neglected. .T’faus,m = z = 0, in Bj. (4.8) and 

V - E [ y ] = h/t^ (4. 10) 

G(y)= 1/’^ (4.11) 

where = the period of record 

= I (4.12) 

and “ refers estimated value and E [ . ] is expectation. 
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Neglect of historic evidence will not chan^ U sufficiently 
but will increase disperssion. The above method of estimation leads 
to a point estimation of y for different magnitude ranges. To 
obtain a continuous variation of y the values obtained from Bii. (4.10) 
are fitted to the Bj., (4, 1). It is to be noted that this fit will be 
identical with the magnitude-frequency law given in Chapter 2. 

The coefficient of variation of p, c(vi), Bq. (4. II) is 
fitted with a Horel function 

c(y) = ejip S-n M + a^M ] (4.13) 

The analysis of variance table is given in Table 4.1. The 
fitted plot is shown in lig. (4.2). The 95^ confidence interval of 
coefficients a^, a^ and a^ are (5.826 _+ 2.605), (-17.097 + 3,345) 
and (3.653 + 0.558) respectively. 

The local seismicity of a smaller volume, say V’ , is obtained 
assuming y’/ y and y are statistically independent and y’ refers 
to the intensity of arrival for volume V’. Bie spatial variation of 
seismicity within Indian peninsula is accounted for by y’/ y . It 
can be shown from Bg,. (4.l) that 

yt = v'/v y ( 4 . 14 ) 

o^( y') = y'/y ) + C^(y»/y ) G^(p) + o\v) (4.15) 

where C^(,) is coefficient of variation square, Esteva^^^) jjas 
shown that by restriction of condition V T=:->y' y to Bq. (4.15) 
yields the solution 
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Table 4. 1 

Analysis of Variance Table for Multiple linear Regression 
Model - g(p) = exp [ a^+ to M + a^M ] 


Number of sample = 

23 




Source of Variation Sum of square 

Degree of 

Root Mean 

P-ratio 



freedom 

square 


hie regression 

30. 1023 

2 

3.8796 

530.2 

Due residual 

0.5677 

20 

0.1685 


Total 

30.6700 

22 



Multiple correlation coefficient 

square = 0.981 


i Coeff. a^ 

a t -value 

a. 

1 

max(x^) min (x^) 


0 5.826 

1.249 4.66 




1 -17.097 

1.605 10.65 

2.128 1. 

,386 


2 3.653 

0.268 13.65 

8.400 4. 

,000 


Covariance of coefficients 





1.561 -1.987 

— 

0.325 



Cot [ ^ ] = 

-1.987 2.575 

-0.427 




0.325 -0.427 

0.072 ■ 




The null hypothesis: a^= a 2 = 0 is rejected at P < 0.001 since 

P = 530.2 > ^ ^ = 9*5 and the null hypothesis re^jding 

o .999 ; 2,20 

a^=0, i = 0, 1,2 is rejected at P--value much less than 0,001. 
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C^( v»/v ) = (Y/r)'^ -1 (4.-16) 

Use of Bg.s. (4.15, 4. 16 and 4,1l) yields 

(V/V)^ = n [l + C^( y«)] /(n+l) (4.1?) 

Bor each volume ratio 250 random volme samples of V’ are 
generated within Indiaa peninsula to estimate G ( y’). !Ehis estimated 
value is used to obtain Y and shown in iig. (4.3). fhe results of 
this analysis are given in fable 4.2. The 95^ confidence interval 
*of Y is (0.408 + 0.02l). 

The relations developed are used to estimate the intensity of 
arrival process and its coefficient of variation at each grid point, 

4,3 Estimation of Hegional Seisoiicity at Grid Points 

To estimate intensity of arrival process at a grid point, the 
number, r, of earthquakes over a given magnitude threshold originated 
at earthquake source, centred at a grid point, in the time period, T, 
of available seianic data is counted. The intensity of arrival y^, 
and coefficient of variation, c( are obtained from Tigs, (2.4 

and 4,2). The intensity, y’, for earthquake source of volume, Y’ , 
is obtained by substitution in Bj. (4.14). The parameter, m, 

and, I z^, of the natural conjugate Bayes density S3., (4,5) are 

m=l/C^(y.|) (4.18) 

I z^= 1/ [ C^( y.j) ] = m/ i); 

where C^( P') is obtained from Bi. (4.15) and Bj. (4.16) as 


(4.19) 
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Table 4.2 

Analysis of Variance Table for Linear Regression Model 
- £n[n {1+C^(vi')> /(n+l) ]= Yto (v/V) 

Humber of samples = 20 


Source of variation 

Sum of Squares 

Degree of 
Preedom 

Root Mean 

sqi;iare 

P~Ratio 

Lue regression 

37.233 

1 

6.102 

400.76 


Due residual 

1-756 

19 

Total 

38.989 

20 


Multiple Correlation Coefficient Square = 0.955 
y = 0.408 

Oy = 0.012 

t-Talue = 59.93 

Since P = 400.76 > P. ^ 15.4, the null hypothesis 

0 0.999; 1,18 

Y = 0 is rejected at P < 0.001. P~value = 2P [t(l8) > 39-95 ] 

< 2x 0.0005 = 0.001, hence Y = 0 is rejected at P <<0.001. 

A 95^ confidence interval for y is (o.408 + 0.02l) 
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v/v 



F!G.4'3VARIABIL!TY OF LOCAL* SEISMICITY WITHIN 
INDIAN PENINSULA FOR MAGNITUDE, M>5 
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G^( y!,) = {(v/v)^ -1} {1 + C^( V^)} + c^( y^) (4.20) 

The Bayes estimate of under q\:^dratic loss function when 
r earthquakes occured in the earthquake source during a time interval^ 

T, is {i^ and is given by 

y' = (r4m)/( T) (4.21) 

and coefficient of -variation of y’ is from Bq,. (4.9) as 

C( y') = 1/ J (r+m) (4.22) 

The regional seismicity at eadi grid point is estimated from 
Bqs. (4.21 and 4.22) and is reported in Table 4.3. 

4*4 Results and liscussion 

The resiJlts dn terms of peak ground accel8ration,velocity and 
displacement are obtained ntmerically from Bq. (3.24) for a specified 
return period, T = 100 years. The generalized intensity is obtained 

y 

in an iterative procedure wi-fch tolerance such that 

100 years T ^ 10 1 years. 

y 

The genea?alized intensities are obtained (Table 4.3) at each 2° x 2° 
grid points and are plotted in contour forms in Tigs. (4. 4-4.6). The 
general ^ape of the contoiir in the three maps is similar. However, 
spread in contours for displacement and velocity are more than that 
for acceleration. Significant differences are present in the central 
and north-eastern region. Although the tectonic features are not 
taken into accomt eiqjlicitly, a clear indication of Himalayan belt is 
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Table 4.3 


Estimated intensity and coefficient of variation square for earthquake 
arrival process and peak generalized seismic intensity at girid points. 

lo^ormal focal depth distribution, m e (5,9j , b = 150 Kms. and 


H = 600 Kms. 
o 


A S 

arid point y xiO 

{ ^ /<iay) 

C^(y)x • 
10 ^ 

100 years ^ 

return period 

101 years 

Acceleration 

(Cm/Sec^) 

Velocity X 10^ 

(Gta/ Sec) 

Displacement 

xIO^ 

(On) 

1 

8. 287 

849.6 

0.2651 

4.009 

9.101 

2 

8. 287 

649.6 

0. 265 1 

4.009 

9.101 

3 

8.287 

849.6 

0. 265 1 

4.009 

9.101 

4 

78.70 

89.47 

8.290 

89.00 

193.0 

5 

8. 287 

849.6 

0. 2651 

4.009 

9. 101 

6 

8. 287 

849.6 

0.2651 

4.009 

9.101 

7 

78.70 

89.47 

8.290 

89.00 

193.0 

8 

8.287 

849.6 

0. 2651 

4.009 ; 

9.101 

9' 

8. 287 

849.6 

0.2651 

4.009 

9.101 

10 

8.287 

849.6 

0. 2651 

4.009 

9.101 

11 

8.287 

849.6 

0.^51 

4.009 

9.101 

12 

8, 287 

849.6 

0.2651 

4.009 

9.101 

13 

78.70 

89.47 

8.290 

89.00 

193.0 

14 

1346 

5.2307 

43.77 

487.3 

128.8 

15 

994.0 

7.083 

37.86 

415.9 

106.9 

16 

78.70 

89.47 

8.290 

89.00 

193.0 

17 

8.309 

851.5 

0.2651 

4.009 

9-101 

18 

8.309 

851.5 

0.2651 

4.009 

9. 101 

19 

79.06 

esA9 

8.320 

89.31 

193.7 

20 

8.309 

851.5 

0.2651 

4.009 

9-101 

21 

8.309 

851.5 

0.2651 

4.0C9 

9.101 

22 

8.309 

851.5 

0.2651 

4.009 

9. 101 
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Table 4,5 contd. 


23 

8.309 

851.5 

0.2651 

4.009 

9. 101 

24 

8.309 

851.5 

0.2651 

4.009 

9.101 

25 

8.309 

851.5 

0.2651 

4.009 

9.101 

26 

8.309 

851.5 

0.2651 

4.009 

9-101 

27 

8.309 

851.5 

0.2651 

4.009 

9.101 

28 

8.309 

851.5 

0.2651 

4.009 

9-101 

29 

8.309 

851.5 

0.2651 

4.009 

9.101 

30 

1636 

4.326 

47.90 

538.9 

145.1 

31 

291.3 

24.29 

19.72 

210.8 

491.2 

32 

8.309 

851.5 

0.2651 

4.009 

9.101 

33 

8.337 

853.9 

0.2651 

4.009 

9. 101 

34 

8.337 

853.9 

0.2651 

4.009 

9.101 

35 

8.337 

853.9 

0. 265 1 

4.009 

9.101 

36 

8.537 

853.9 

0.2651 

4.009 

9.101 

37 

8.337 

853.9 

0.2651 

4.009 

9.101 

38 

8.337 

853.9 

0.2651 

4.009 

9. 101 

39 

8.337 

853.9 

0. 265 1 

4.009 

9.101 

40 

8.337 

853.9 

0.2651 

4.009 

9.101 

41 

8.337 

853.9 

0.2651 

4.009 

9. 101 

42- 

8.337 

853.9 

0.2651 

4.009 

9- lOi 

43 

8.337 

853.9 

0.2651 

4.009 

9.101 

44 

8.537 

853.9 

0.2651 

4.009 

9.101 

45 

293.1 

24.29 

19.78 

211.6 

493.2 

46 

1076 

6.615 

39.34 

433.6 

1122 

47 

862.6 

8.253 

35.30 

385.8 

978. 2 

48 

8.337 

853.9 

0. 2651 

4. 009 

9.101 

49 

8.370 

856.8 

0.2651 

4.009 

9. 101 

50 

8.370 

856.8 

0.2651 

4.009 

9.101 

51 

8.370 

®6.8 

0.2651 

4.009 

9. 101 

52 

8.370 

856.8 

0.2651 

4.009 

9.101 

53 

8.570 

856.8 

0.2651 

4.009 

9.101 

54 

8.370 

856.8 

0.2651 

4.009 

9. 101 

55 

8.370 

856.8 

0.2651 

4.009 

9.101 
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56 

8.370 

856.8 

0.2651 

4.009 

9.101 

57 

8.370 

856.8 

0.2651 

4.009 - 

9.101 

58 

8.370 

856.8 

0.2651 

4.009 

9.101 

59 

8.370 

856.8 

0.2651 

4.009 

9.001 

60 

80.08 

89.55 

8.402 

90.19 

195.7 

61 

8.370 

856.8 

0.2651 

4.009 

9.101 

62 

1443 

4.971 

45.21 

505.1 

1344 

63 

510.4 

14.05 

26.94 

290.2 

704.1 

64 

80.08 

89.55 

8.402 

90.19 

195.7 

65 

8.409 

860.2 

LA 

O 

4.009 

9-101 

66 

8.409 

860.2 

0.2651 

4.009 

9-101 

67 

8.409 

860.2 

0.^51 

4.009 

9.101 

68 

8.409 

860.2 

0.2651 

4.009 

9.101 

69 

8.409 

860.2 

0.2651 

4.009 

9-101 

70 

8.409 

860.2 

0.2651 

4.009 

9-101 

71 

80.74 

89.59 

8.455 

90.75 

197.0 

72 

8.409 

860.2 

0.2651 

4.009 

9- 101 

73 

80.74 

89.59 

8.455 

90.75 

197-0 

74 

8.409 

860.2 

0.2651 

4.009 

9.101 

75 

8.409 

860.2 

0.2651 

4.009 

9. 101 

76 

80.74 

89.59 

8.455 

90.75 

197.0 

77 

8.409 

860.2 

0.2651 

4.009 

9.101 

78 

514.7 

14.05 

27.07 

291-6 

707.9 

79 

153. 1 

47.25 

13.27 

141.8 

317.9 

80 

225.4 

32.09 

16.93 

180.8 

414.5 

81 

8.454 

864.1 

0.2651 

4.009 

9.101 

82 

8.454 

864.1 

0.2651 

4.009 

9-101 

83 

8.454 

864.1 

0.^51 

4.009 

9.101 

84 

1031 

7.084 

38.54 

424.0 

109.3 

85 

8.454 

864.1 

0.2651 

4.009 

9 -10 1 

86 

8.454 

864.1 

0.2651 

4.009 

9-101 

87 

8,454 

864.1 

0.2651 

4.009 

9-101 

88 

154.6 

47.27 

13.35 

14.27 

320.1 
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89 

8.454 

864. 1 

90 

8.454 

864. 1 

91 

8.454 

864.1 

92 

8.454 

864.1 

93 

8.454 

864.1 

94 

154.6 

47-27 

95 

81.51 

89.63 

96 

373-7 

19.55 

97 

8.505 

868.6 

98 

8.505 

868.6 

99 

82.38 

89.68 

100 

82.38 

89.68 

101 

8.505 

868.6 

102 

8.505 

868.6 

103 

8.505 

868.6 

104 

8.505 

868.6 

105 

8.505 

868.6 

106 

8.505 

868.6 

107 

8.505 

868.6 

108 

8.505 

868.6 

109 

8.505 

868.6 

110 

8.505 

868.6 

111 

451.8 

16.35 

112 

156.3 

47.28 

113 

83.38 

89.73 

114 

8.563 

873.7 

115 

8.563 

873.7 

116 

83.38 

89.73 

117 

83.38 

89.73 

118 

8.563 

873.7 

119 

8.563 

873.7 

120 

8.563 

873.7 

121 

8.563 

873.7 


0.2651 

4.009 

9.101 

0.2651 

4.009 

9.101 

0.2651 

4.009 

9.101 

0.2651 

4.009 

9-101 

0.2651 

4.009 

9.101 

13.35 

142.7 

320.1 

8.517 

91.40 

198.5 

22.73 

243.5 

577-3 

0.2651 

4.009 

9.101 

0.2651 

4,009 

9.101 

8.587 

92. 14 

200.2 

8.587 

92. 14 

200,2 

0.2651 

4.009 

9.101 

0.2651 

4.009 

9.101 

0.2651 

4.009 

9. 10 1 

0.2651 

4.009 

9.101 

0.2651 

4.009 

9-101 

0.2651 

4.009 

9.101 

0.2651 

4.009 

9.101 

0.2651 

4.009 

9.101 

0.2651 

4.009 

9.101 

0.2651 

4.009 

9.101 

25.23 

271.1 

651.7 

13.45 

143.7 

322.5 

8.665 

92.97 

202.2 

0.2651 

4.009 

9,101 

0.2651 

4.009 

9.101 

8,665 

92.97 

202-2 

8.665 

92.97 

202.2 

0.2651 

4.009 

9.101 

0.2651 

4.009 

9.101 

•.2651 

4.009 

9.101 

0.2651 

4.009 

9.101 
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122 

8.563 

873.7 

0.2651 

4.009 

9.101 

123 

158.2 

47.29 

13.55 

144.8 

325.3 

124 

83.38 

89.73 

8.665 

92.97 

202.2 

125 

8.563 

873.7 

0.2651 

4.009 

9.101 

126 

158.2 

47.29 

13.55 

144.8 

325.3 

127 

532.3 

14.06 

27.55 

297.0 

723.0 

128 

158.2 

47.29 

13.55 

144.8 

325.3 

129 

84.49 

89.79 

8.753 

93.90 

204.3 

130 

84.49 

89.79 

8.753 

93.90 

204.3 

131 

84.49 

89.79 

8.753 

93.90 

204.3 

132 

8.628 

879.3 

0.2651 

4.009 

9.101 

133 

8.628 

879.3 

0.2651 

4.009 

9.101 

134 

84.49 

89.79 

8.753 

93.90 

204.3 

135 

8.628 

879.3 

0.2651 

4.009 

9.101 

136 

8.628 

879.3 

0.2651 

4.009 

9.101 

137 

84.49 

89.79 

8.753 

93.90 

204.3 

CD 

84.49 

89.79 

8.753 

93.90 

204.3 

139 

84.49 

89.79 

8.753 

93.90 

204.3 

140 

8.628 

879,3 

0. 2651 

4.009 

9.101 

14 1 

312.1 

24.31 

20.52 

219.5 

513.8 

142 

919.0 

8.255 

36.43 

399.0 

1017 

143 

2436 

3.114 

57.27 

660.6 

1830 

144 

615.6 

12.32 

29.74 

321.7 

792.3 

145 

239.8 

32.12 

17.57 

187.7 

432.1 

146 

239.8 

32. 12 

17.57 

187.7 

432.1 

147 

85.74 

89.85 

8.851 

94.93 

206.7 

148 

85.74 

89.85 

8.851 

94.93 

206.7 

149 

8.700 

885.5 

0.2651 

4.009 

9.101 

150 

239.8 

32.12 

17.57 

187.7 

432.1 

151 

8.700 

885.5 

0.2651 

4.009 

9.101 

152 

8.700 

885.5 

0.2651 

4.009 

9.101 

153 

8.700 

885.5 

0.^51 

4.009 

9.101 

154 

8.700 

885.5 

0-2651 

4.009 

9. 10 1 
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155 

8.700 

885.5 

156 

259.8 

32.12 

157 

1550 

4.972 

158 

1087 

7.086 

159 

2859 

2.694 

160 

2782 

2.769 

16 1 

522.2 

24.52 

162 

165.5 

47.54 

165 

8.780 

892.4 

164 

8.780 

892.4 

165 

8.780 

892.4 

166 

8.780 

892.4 

167 

87.15 

89.92 

168 

87.15 

89.92 

169 

245.8 

32.13 

170 

400.5 

19.56 

171 

557.2 

14.06 

172 

245.8 

32. 13 

175 

1889 

4.147 

174 

1541 

5,843 

175 

949.0 

8.256 

176 

1968 

3.982 

177 

1206 

6.6 18 

178 

2164 

3.689 

179 

727.1 

10.98 

180 

88.68 

90.00 

181 

88.68 

90.00 

182 

168.5 

47.37 

185 

168.5 

47.37 

184 

2004 

3.982 

185 

567.5 

14.06 

186 

88.68 

90.00 

187 

647.5 

12.33 


0.2651 

4.009 

9.101 

17.57 

187.7 

432.1 

46.72 

524.1 

1404 

39.54 

436.0 

1130 

61.39 

715.7 

2037 

60.67 

706.0 

2004 

20,89 

223.6 

524.6 

13.95 

149.1 

535.7 

0.2651 

4.009 

9.101 

0.2651 

4.009 

9. 101 

0.2651 

4.009 

9. 101 

0.2651 

4.009 

9.101 

8.958 

96.07 

209.4 

8.958 

96.07 

209.4 

17.75 

189-6 

436.8 

23.62 

255.5 

603.5 

28.23 

304.6 

744.2 

17.75 

189.6 

436.8 

51.15 

580. 1 

1585 

45.69 

486.5 

1284 

37.01 

405.8 

1038 

52.09 

592.3 

1625 

41.55 

460.2 

1203 

54.55 

621.6 

1721 

32.40 

352.1 

879.4 

9.076 

97.52 

212.3 

0.265 1 

4.009 

9. 101 

14.12 

150.8 

339.9 

14.12 

150.8 

339. 9 

52.52 

597.8 

1643 

28.50 

307.7 

752.8 

9.076 

97.32 

212.2 

30.53 

330.6 

817.7 
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188 

248.3 

32.14 

189 

647.3 

12.33 

190 

248.3 

32. 14 

191 

4159 

1.919 

192 

2403 

3.321 

193 

497.4 

16.37 

194 

■>067 

7.628 

195 

578.8 

14.06 

196 

90.38 

90.08 

197 

171.8 

47.39 

198 

171.8 

47.39 

199 

904.5 

9.001 

2D0 

660.3 

12.33 

201 

334.6 

24.35 

202 

660.3 

12.33 

203 

497.4 

16.37 

204 

497.4 

16.37 

205 

1800 

4.522 

206 

171.8 

47.39 

207 

741.7 

10. SB 

208 

578.8 

14. #6 

209 

175-4 

47.42 

210 

591.4 

14.07 

211 

175.4 

47.42 

212 

341.8 

24.34 

213 

591.4 

14.07 

214 

175.4 

47.42 

215 

425.0 

19.57 

216 

92.26 

90.17 

217 

258.6 

32- 16 

213 

175.4 

47.42 

219 

92.26 

90.17 

220 

9.070 

917.2 


17.94 

191.7 

442.1 

30.53 

330.6 

817.7 

17.94 

191.7 

442. 1 

71.91 

861.8 

2545 

56.92 

656.0 

1834 

26.58 

286. 1 

692.7 

39.19 

431.7 

1117 

26.80 

311.1 

762.3 

9.206 

98.69 

215.4 

14.29 

152.7 

344.5 

14.29 

152.7 

344.5 

36.14 

395.6 

1007 

30.84 

334.2 

827.9 

21.35 

228.6 

537.6 

30.84 

534.2 

827.9 

26.58 

286.. 1 

692.7 

26.58 

266.1 

692.7 

50.04 

566.0 

1559 

14.29 

152.7 

344.5 

32.73 

355.9 

890.3 

28.80 

311.1 

762.3 

14.49 

154.7 

349.6 

29.13 

314.7 

772.6 

14.49 

154.7 

349.6 

21.61 

251.3 

545.0 

29. 13 

314.7 

772.6 

14.49 

154.7 

349.6 

24.40 

262.0 

626.8 

9.347 

100.2 

218.9 

18.38 

196.5 

454.3 

14.49 

154.7 

349.6 

9.347 

■K0O.2 

218.9 

0.2651 

4.009 

9 . 1< 
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221 

341-8 

24.34 

21.61 

231.3 

545.0 

222 

924. 1 

9.002 

36.52 

400.2 

1021 

223 

341.8 

24.34 

21.61 

231.3 

545.0 

224 

175.4 

47.42 

14.49 

154.7 

349.6 

225 

349.8 

24.34 

21.89 

234.4 

553.2 

226 

1116 

7.629 

40.04 

442.0 

1148 

2 ZI 

1457 

5.845 

45.41 

507.6 

1352 

228 

690.4 

12.33 

31.55 

342.4 

851.3 

229 

264.6 

32.18 

18.63 

199.2 

461.2 

230 

520.1 

16.38 

27.21 

293.3 

712.5 

231 

264.6 

32.18 

18.63 

199.2 

461.2 

232 

690.4 

12.33 

31.55 

342.3 

851.3 

233 

94.33 

90.26 

9.501 

101.8 

222.7 

234 

9.186 

927.0 

0.2651 

4.009 

9.101 

235 

94.33 

90.26 

9.501 

101.8 

222.7 

236 

179.4 

47.44 

14.70 

157.0 

355.1 

237 

775.5 

10,98 

33.48 

364.5 

915.5 

238 

179.5 

47.44 

14.70 

157.0 

355.1 

239 

9.186 

927.0 

0.2651 

4.009 

9.101 

240 

264.6 

32.18 

18.63 

199.2 

461.2 

241 

620.5 

14.07 

29.87 

323.1 

796.3 

242 

2192 

3.983 

54.66 

626.5 

1735 

243 

17210 

0.5073 

125.8 

1663 

5554 

244 

1406 

6.209 

44.67 

498.5 

1323 

245 

358.6 

24.35 

22.20 

237.8 

562.1 

246 

533.2 

16.38 

27.58 

297.3 

723.8 

247 

96.62 

90.36 

9.669 

103.6 

226 .9 

248 

96.62 

90.36 

9-669 

103.6 

226.9 

249 

9.312 

937.7 

0. 2651 

4.009 

9.101 

250 

533.2 

16.38 

27.58 

297.3 

723.8 

251 

96.62 

90.36 

9.669 

103.6 

226.9 

252 

183. 9 

47.47 

14.93 

159.4 

361.2 
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253 

96.62 

90.36 

9.669 

103.6 

226.9 

254 

96.62 

90.36 

9.667 

103.6 

226.9 

255 

96.62 

'90.36 

9.669 

103.6 

226.9 

256 

96.62 

90.36 

9.669 

103-6 

226.9 

257 

188.8 

47.50 

15.18 

162, 1 

367.7 

258 

996.1 

9.005 

37.90 

4 l 6.3 

1070 

259 

6288 

1.426 

85.06 

1053 

3234 

260 

2880 

3.115 

61.58 

718.2 

2046 

26 1 

2611 

3.436 

59.02 

683.9 

1929 

262 

1534 

5.846 

46.52 

521.4 

1395 

263 

727-0 

12.34 

32.40 

352.1 

879.3 

264 

188.8 

47.50 

15. 18 

162. 1 

367.7 

265 

9.449 

949.3 

0.2651 

4.009 

9-101 

266 

9.449 

949.3 

0.2651 

4.009 

9.101 

267 

99.15 

90.47 

9.851 

105.5 

231.4 

268 

9.449 

949.3 

0.2651 

4-009 

9.101 

269 

278.5 

32.20 

19.21 

205.3 

477.0 

270 

9.449 

949-3 

0.2651 

4.009 

9. 101 

271 

278.5 

32.20 

19.21 

205.3 

477.0 

272 

368.2 

24.36 

22.54 

241.5 

571.8 






FIG. 4*6 PEAK 


DISTRIBUTION OF FOCAL DEPTH 



reflected by the contour maps. It is expected that a substantial 
diange in contour pattern is likely to occur in the proximaty of 
geological faxAlts ■when faults are taken into account in the seismic 
risk analysis. The focal depth distribution is identical at each 
grid point in this analysis. This distribution reflects shallow 
focus earthquakes. However, from past focal depth data it is observed 
that grid point 243 has a dominant intermediate focus earthquakes. 

The analysis with identical focal depth distribution will overestimate 
generalised intensity at this grid point. 

4.5 Extrapolation of Generalized Intensity at a location 

The contour maps prepared in the previous section provide 
estimates of peak ground acceleration, velocity and displacement at 
any location in the Indian peninsula for lOO-year return period. The 
return period is estimated from Bj. (3.24) which can also be- written 
as 

hT P[ Y > y ] = 1 (4.23) 

Substituting t = T from Bq. (4.23) into Bq. (3.22) gives 

y 

By (yjT )= exp (-1) = 0.37 (4.24) 

max 

Thus a generalized intensity based on T -year return period has 
a probability of 0.63 being exceeded. 

While -using the contour maps to establish the desiga loads 

for a structure at a location, it is quite often necessary to 

exrtrapolate the generalized intensity for a 100-year return 
period given in the map to some other return period 
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(service life). It may also become necessary to determine the 
generalized intensity for a specified exceedance probability over a 
given period (service life). 

Since the maps presented in the previous section have been 
drawn on the basis of seismic risk determined from modtilur earthg,uake 
source at each grid point, these extrapolation can be done in a very 
simple manner as explained below. 

Bigs. ( 4.7 - 4 . 9 ) give the plot of generalized intensity 
against J.n( u T ) for 0.63, 0.25, 0.10 and 0.05 values of exceedance 
probabiliti^ at any point in the Indian peninsula. 

4.5.1 Extrapolation for Return Period 

let the design intensity required for a retircti period be 
I = CT. Generalized intensity is determined by solving for y from 

y 

the equation 

An( vl) + an 0 = - an p[r>y] (4.25) 

which is obtained from Eq. ( 4 . 23 ). 

Extrapolatiaa can be done in the following steps as shown in 
Bigs. ( 4.7 - 4,9): 

1, Determine the generalized intensity for a T = 100-year 
return period frcm the contour map. 

2. locate point B on the curve with p = O .63 Add an C to the 
ordinate and locate point C on the curve with p = 0,65. 





y (cm/sec) 


FIG.4-8 VARIATION OFajT WITH PEAK GROUND VELOCITY , Y 




■iwiiiiii*^*^^*^fciiiiaiiiiiii®i*^ 
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3. 


Project vertically from C to the abscissa to get the 

extrapolated generalized intensity for 1 year return 

If 

period. 


Example : 


let T =200 years at a location where 100-year peak acceleration 

O' 

2 

is 52.5 cm/sec , Following the above procedure the intensity is obtained 
2 

as 72 cVsec . 

4.5.2 Extrapolation for specified exceedance probability 


Let the generalized intensity at a location be required for a 
specified exceedance probability p, over a 1 year return period. From 

Bl, (3.22) 

P=1-Ey (y;T) = 1-Pj^(0,!r) = l-exp { - yPP [ Y >y]} (4.26) 

max 

Graphically the extrapolation can be done in the following steps: 

1, Determine the generalized intensity for a T-year return 
period from the contour map, 

2, locate the point B on the curve marked p = O.63. Project 
horizontally to the curve corresponding to required p and 
locate the point D on it. 

3, Project vertically from D to the abscissa to get the extra- 
polated generalized intensity for the specified exceedance 
probability, p, and return period, T. 
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Example 

2 

Let peak acceleration is 52.5 cm/sec at a location for 100-year 

return period. Bellowing the above procedure the intensity is 

2 

172.5 cm/sec if design intensity for 100-year return period with 
exceedance probability, p = 0.05, is required. The exceedance 

(55) 

probability, p = 0.05, is the specified level for nuclear facilities , 

4.5.3 Extrapolation for specified exceedance probability and 
specified return period. 

If the period T is different from the return period in the 
contour map, adopt upto step 2 of the procedure for extrapolation for 
return period. Then 

3 . Project horizontally from the point C to the specified p-cinve 
and locate point E on it. 

4. Proj ect vertically from E to the abscissa to obtain the 
required generalized intensity. 

Example 

Let p = 5/^ over a 200-year period at a site vdiere 100-year 

2 

return period peak ground acceleration = 52.5 cm/sec . The extra- 

2 

polation as ex 5 )lained above yields 220 cm/sec . 

4.5.4 Extrapolation fer magnitude cutoff 

Let the generalized intensity at a location is required vshen 
magnitude, m e (Sjm^), for a T-year return period. Graphically this 
extrapolation is achieved from ligs. (3.11-5.13) by following procedure: 
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1# Determine the generalized intensity for a T—year return 

period from contour map. 

2# Locate the point A on the curve marked, m^= 9# Project 

horizontally to the m^-^curve to mark point B on it. 

3* Project vertically from B to the abscissa to obtain 

extrapolated value. 

Example 

Let CO at a location with peak acceleration of 125 
2 

can/ sec for 100-year return period. By above procedure 

2 

the extrapolated value for 100-year is 135 cm/sec . 

Any additional extrapolation may be carried out by choosing 

the proper procedure (4.5, 1-4, 5. 3) as per reqiiirements. 

As explained above, the generalised intensity obtained frcm the 
100 year return period contour maps can be easily extrapolated for 
different retinrn periods and exceedance probability using a set of 
curves. Such a sin^jle procedure vdiich can be easily incorporated in the 
codes, has been possible due to the assumption that the peninsula is 
homogeneous in tectonic features and use of a modular earthqmke source, 
Hiis is the most significant feature of the contour maps presented in 
this section. 



CHAITER 5 


SEISfllC EE SIC AMIYSIS OP INDIM PMIISULA-A MODEL NOTHCivIOGMEDUS H 

TECTONIC FEATURES 


5.1 Introduction 

The contour maps of peak generalised seismic intensities prepared 
in the previoiis chapter are based on homogeneous source model vdiich 
ignores the geotectonic information. In estimating the seiaaic risk 
the nat-ural conjugate Bayes density at each grid point is assumed to 
hare identical values of the parameter. In this chapter, the earthquake 
sources are identified based on available geotectonic information. The 
natural conjugate Bayes density at different grid points is determined 
based on its location in the Indian peninsula. Por this purpose the 
tectonic features are idealized. CJontour maps of generalised intensities 
are prepared and c<mpared with the maps based on homogeneotis tectonic 
features, 

5.2 Tectonic Peatures of Indian Peninsula 

The tectonogenesis of the Indian peninsula can be pieced together 
on the basis of the continental drift theory. Starting with Pangae in 
the late paleozoic or early Mesozoic period , the formation of Indian 
peninsula can be traced to its present state. The Indian peninsula is 
a part of Eurasia composed of mountains of Aitaid and Alpide systems. 

The geological structure varies considerably within both the systems. 

Hie Indian peninsula can be classified broadly in three different regLons 
as south of Vindhyas, the regions of Alpide systems and the regions of 
Aitaid system. 
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Hie south of Viridhyas is a precambrian shield and coasist of 
Deccan and Sri lanka. The Deccan is a peneplain composed of ancient 
metamorphic rocks and partly covered by Late Cretaceous basalts in the 
western region. The south of Tindhjas has never been submerged totally 
■under the sea. Glie region is a shield of great rigidity unaffected by 
any folding. The structural disturbances in this region are'block 
movement* i.e. displacement in the vertical direction of large land 
masses between radial fa-ults within the region. The mountains in the 
Deccan are of 'relict type’. 

The regions of Alpide system consist of Indo-Gangetic plain 
and Himalya and have evidence, in the form of marine deposits starting 
from Cambrian period, to indicate that the land masses have remained 
beneath the sea for the greater part of its existence. The Himalaya 
j a part of Alpide system , is formed during Tertiary age and is composed 
of comparatively weak and yielding type of rocks. Himalayan tectonics 
have been interpreted as the compression bet'^.veen the Asian main land 
and shield of India. Oliese regions are associated with large-scale 
concentric folding and thrusting. The characteristics of earth-movements 
are of flexible nature. The wrinkles and folds formed by the tangential 
thrust on more or less linear zones of earth's surface are the momtain- 
chains of this region. (Hie indo-gangetic plain is originally a deep 
depression between north of Vindhyas and Himalaya and is covered up 
by allu'vial deposits of river-clays and silts. The Geological Survey 
of India considers this great depression as a ’ibre-deep’, fronting 
the Himalayan earth waves, and a sagging part of the Deccan arresting 
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the southward advance of moimtain-^vaves, A second theoiy considers 
the depression as a ' Rift-Valley' 

A sequence of parallel ranges and high plateaus in the north 
of Himalaya and a part of Altai region at the edge of Siberia shield 
form the region of Altaid systems, The Altaids are mostly Paleozoic 
ranges with intense and wide spread metamorphism. Ghe plateau of 
Tibet is a stable high platform surrounded by seismically active ranges. 

The major tectonic features of India and epicenters of the 
earthquakes in Indian peninsula are diown in lig, (5.l). The tectonic 
features have been idealiaed in this investi^tion following closely 
the suggestion of Srivastava^^^^ and are shown in Tig. (5.2). The 
idealized map consists of twentyone faults and ei^t tectonic units. 
Among ^1 the tectonic units the 'erogenic unit' and the 'foredeep and 
marginal depression unit' have undergone maximum tectonic disturbances 
in comparison to remaining units, Hie geometric property and the 
associated seismic activity of each fault and tectonic unit are given 
in Table 5.1. The significant tectonic features of the idealized map 
are; 

(i) Orogenic unit containing Shillong massif and major thrust 
of Himalaya 

(ii) itoredeep and Marginal depression unit of Himalaya containing 
Patna, Morad^ad, -lucknow faults, 

(iii) West Coast and Earmada-Tapti unit containing the tectonic rift 
that crosses India frem the Gulf of Cambay to Bihar. 
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(It) Gondowana Rift unit - a part of the southern shield with 
Mesozoic or early cenozoic sedimentary cover. 

(v) Shield unit containing non active ancient faults. 

(vi) Bay of Bengal unit containing seismically active Andaman and 
Nicobar island of the Burma arc. 

(vii) Hie Arabian sea unit containing Chagos-Iaccadive ridge of 
ballastic magma. 

(viii) Sri Banka unit containing part of Sri Lanlca which is a part 
of Indian shield 

5.3 Earthquake Sources 

Hie earthquake sources in the seismic risk analysis, based, on 
aon-homogeneous tectonic features, consist of faults, viiiich are area 
sources, and volume sources. (The risk analysis for volume sources has 
been treated in the previous chapter, The probability of exceedance 
of generalised intensity (Sq. 3.19) for faults, vdiidi is idealized as 
segment of great circle, is obtained by passing a plane throu^ the 
fault and transforming the plane to the equatorial position. The 
transformation enabling to rotate a fault plane to the position of the 
equator is given in Appendix A5.1. The formulae (Bqs. 3.25-3. and 
3.30-3.34) for probability of exceedance of peak generalised seianic 
intensity for assumed distribution of focal depth, H, and magnitude, M, 
of earthquakes with suitable changes in the probability density, fp^(*)» 
and probability distribution, Pjj(,), of focal distance, E, remain valid 
even in case of area source. The detailed derivation of these functions 
are given in Appendix A5.2, 
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5.4 !r~year Gteaeralised peak Seismic Intensity 

The faiiLts within a volume source at a site are assumed to be 
mutually statistically independent and also statistically independent 
with respect to the volume source, let the number of faults at a 
particular site be n , The probability distribution function of 

maximum peak generalized seismic intensity can be fbrmally written as 


n+l 

Ey = n P [ Y < y; t ] 

Vx i=i 


(5.1) 


in vhidi Y ^ ^ is 1he maximum peak generalised intensity for the 
source, i. The source n+i represent the volume source at the site 
and the sources i=1,.,,,n represent n number of area sources associated 
with the feults. If {L^^Ct); t ^0, i=l,...,n+l } are the number of 
earthquakes owing to sources which produce peak generalised seismic 
intensity at a site greater than a value, y, during a time inteival, t, 
the Bi. (5.1) can be farther sin^lified as 

n+1 n+1 

Ey (yj't) = n (o,t) = exp [ - I Uj^t Sy_(y) ] ( 5 . 2 ) 

max i=1 i i=1 i 


in vhich is the intensity of earthquiake arrival in the ith source 

and Sy (y) = P [Y^ > y ] is the probability of exceedance of peak 
^i 

generalised seismic intensity, Y^, at a site owing to ith source. The 
probability distribution of annual maximum peak generalized seismic 
intensity is obtained by substitution of t=1 in the Bj. (5.2): 


max 


= exp [- ^ 3 

1 1 - I Sy (y) I Py Sy (y) < < 1 3 

i ’ i 


(5.3) 
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which yields 

1 

Olie average return period, T , of peak generalised seismic 

V 

intoasity etxceeding a value, y, is given by 

ly = l/l > y;1 ] = 1/ { I Sy (y)} (5.4) 

The T-year generalised seismic intensity at a site for a given 
time period T is obtained by solving Bi. (5.4) for y. 

t/ 

5*5 Estimation of parameters of earthquake arrival process for 
various seismo-tectonic provinces and faults* 

The intensity and the coefficient of variation of intensity of 
arrival process in the required magnitude range are estimated from 
Bis. (4. 1 and 4. 1 2) for the Indian peninsijla. The volume of each 
seismo-tectonic province is calculated. Using Bis, (4.13-14) and 
variability of seismicity in Indian peninsula (iig. 4.5) > the intensity 
and coefficiaat of variation of intensity of earthquake process for 
each seianotectonic province is obtained. These results are reported 
in Table 5.1. 

Assuming each favat is of unit thickness, the intensity and 
coefficient of vacation of intensity of arrival process are estimated 
following the procedure of estimation for seianotectonic province. Hie 
results are diown in Table 5.1. 
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Table 5.1 


Geometric and associated Seismic activity of Tectonic mits and Tlaiolts 

of Indian PeninsxXLa. 

Tectonic Units 


Hame 

Y* 

p X lo'^ 
(!5!/day) 

C^(v) X 

Orogenic unit 

110.01 

295.0 

5.085 

Pore deep and Marginal 
depression imit 

20.29 

54.4 

20.060 

West coast and Harmada- 
Tapti unit 

16.12 

43.2 

23.021 

Gondowana rift mit 

8.73 

25.4 

32.347 

Shield mit 

64.76 

173.7 

8.723 

Bay of Bengal xmit 

43.63 

117.0 

11.997 

Arabian Sea unit 

54.96 

93.7 

14.078 

Sf*i Banka mit 

2.12 

5.7 

■ 65.582 


Actual voltxme of unit = V* {JaCs-Hq) ^ ®o 

Eaults 


fault Number 

Length 

(Hus.) 

y 

(#/day) 

c"(p) 

1 

1209-7 

8.38 

41.19 

2 

1512.3 

10.48 

37.51 

3 

2293.0 

5.88 

31.50 

4 

848.5 

5.88 

47.75 

5 

258.0 

1.79 

78.22 
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lable 5.1 contd. 


6 

260.4 

1.80 

77.93 

7 

1846.9 

12.79 

34.50 

8 

388.7 

2.69 

66.03 

9 

279.3 

15.88 

31.50 

10 

342.6 

2.37 

69.57 

11 

464.6 

3.22 

61.33 

12 

205.7 

1,43 

85.90 

13 

265.6 

1.84 

77.30 

14 

233.9 

1.62 

81.46 

15 

279.2 

1.93 

75.71 

16 

316.2 

2. 19 

71.91 

17 

3275.5 

22.69 

27. 10 

18 

2575 

17.84 

29.99 

19 

2079.9 

14.41 

52.82 

20 

218.8 

1.52 

83.74 

21 

302.1 

2.09 

7 3.29 
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5,6 Estimation of f&rameters of Arrival Process at Osrid. Points. 

Ihe procedure of Section 4, 3 is used to estimate the 
parameters of volume source at the site. However, the number, r, 
of earthqmkes originated in the volme source at site exclude 
earthquakes originated in the fetats. All earthquakes \iiich have 
epicenter within 50 kms. from a faiXLt are assumed to be caused by the 
fault. She effective length of a fault at the site is assumed to be 
the laagth of the feult intersected by a circle with centre at the 
site and having a ladius of 150 loas, at the surface of earth. Si irrii a-p 
procedure for volume source (Section 4. 3 ) is adopted to estimate the 

parameters of the arrival process for a fault associated with the site. 

5»7 Eesults and Discussion 

iPhe results of the seismic risk analysis for 100-year return 
period at each grid point is obtained numerically adopting iterative 
procedure suid are plotted in contour forms in H.gs, (5. 3-5. 5). Hie 
tolerance of one year in the hi^er side is used for return period. 

Hie regional effects of active faults are reflected in the contours. 
Comparison of Pigs. (4.4-4. 6 and 5. 3-5. 5) shows that the shape of 
contours for differ aat regions have changed from homogeneous mod^ 
analysis. Hie reason is the acceptance of different Bayes natural 
conju^te density for various seiaaotectonic provinces and that of 
faults. However, values of generalised intensities are quite ccmpairable 
in most of the parts of Indian peninsula for homogeneous and nonhomogeneous 
analysis. In the Koyna region (grid point 84) contours have oilar^d 
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because of the fault (number l) , Slirther the accoxmting of faults 
(number 11-14) causes considerable changes in the shape of conto^ars 
in the central part of India. The contours in the Gondowana rift 
unit (grid point 88) become parallel to the tectonic province. In 
general the assumption of nonhomogenity of tectonic feature builds up 
generalised seismic intensity in the proximity of faults and makes 
the contours parallel to the seism.otectonic province. 

In view of the above, it is concluded that the seismic zoning 
map should be drawn on the basis of seismic risk analysis with sources 
identified on the basis of tectonic features. However, extrapolation 
of the contour maps for the nonhomogeneous model is not possible because 
of the presence of effects of various sources at a ate in the denominator 
of the Bj. ( 5 . 4 ) on vhida computation of T-year intensity is based, A 
fresh computation of seismic risk analysis is necessaiy for information 
other than that presoited in the contour maps in ELgs. (5. 3-5. 5). 



CHAPTER 6 


EPEECT OP SCATTER IH THE ATTMUATIOH LAWS 01 SEISMIC RISK 
6.1 Introduction 


In the preceeding chapters the attenuation laws have been assumed 

to be deterministic. The generalized attenuation law proposed by Kanai^®^^ 

and Esteya and Eosenblueth^^"*^^ (Bg[. 3, 10) is used in this thesis. The 

constants involved in the attenuation laws are taken according to the 

values suggested by Esteva and Eosenblueth. An excellent review on similar 

(23) 

laws is given in Ambraseys paper the parameters of the law are 
determined from regression analysis. The data used by Esteva and Bosenblueth 
shows large scatter around the best fit equation as shown in Pigs. (6.1 and 
6.2). This is because the geological conditions of the data are signifi- 
cantly different. The results of Esteva^"^^^ and show that 

the effect of scatter on the seismic risk is significant. In this chapter 
two models are proposed to incorporate the effect of scatter on seismic 

(42) 

risk besides the Gaussian model due to Esteva ' , Steps to incorporate 
the effect of scatter in the zoning maps are indicated. 


6,2 A Gaussian Model 

(24) 

Esteva and Eosenblueth obtained the parameters of attenuation 
laws by fitting equation to the data. The theory of least sqaere fit 
assumes the existance of an unobservable error. The attenuation law with 
the unobservable error is given by 

Y = C^ exp [Cg M - C^ £n (r + C^) + e] (6.1) 

2 

vshere e is a Gaussian random variable with mean, % , and variance, o • 

C^, Gg and 0^ are empirical constants. The inclxision of the parameter e 



Y exp (-0-8N 



R (Kms) 


FIG.6'1 ATTENUATION. LAW FOR PEAK GROUND 
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changes the exceedance probability of Y in Bi. (3. 19) and therefore the 
seismic risk, 

6.2.1 Probability of Generalised intensity 

Ihe probability of generalised seismic intensity is obtained 
using the fact that the random variables M, R and e are mutually independent. 
The probability of generalised seismic intensity is given by 

P[Y>y] =//P[Y>yl z,w] dz dw (6.2) 

where 

Z = exp (C^ M) (6.3) 

and W = (r + C^) ^ (6.4) 

and is the joint probability density function of Z and W, Using 

the property of statistical indepaidence, the conditional prdaability in 
Eq,. (6,2) can be evaluated as 

P [ Y > y lz,w] = P [esp (e) > y/(zw) lz,w] 

= P [exp (e) > y/(zw)] 

= 1 - (in y-jnz-Jlnw - ^/a} (6.5) 

vUieve $ (•) is the probability distribution function of H(o,i), The 
probability (Sq. 6.2) can be evaluated by substituting Bq. (6.5) into 
it and using the distribution law of Z and W. 

6.2. 1.1 fte probability of generalised intensity for magnitude cutoff 
at lower level i. e. m e (m^,“) 

The probability density function of Z is givm by 



132 


= 6 e3cp[ B{m^+(iln - jjn z)/G^}]/(c^z) ; z e (c^ exp (c^ m^)^) 

... ( 6 . 6 ) 

Substitution of Bj. (6.5) into Eq. (6,2) and use of statistical 
independence yields 

W 00 

r[Y>y] = 1 ^ J J ${(Aa y - £n z - jin w - 5 )/g} f^(z) f^(w) dz dw 


Wi z^ 


where 


eap (C^ m^) 


^2 == ^o^ 


■•°5 


... (6.7) 

( 6 . 8 ) 

(6.9) 

( 6 . 10 ) 


and is the mazinium focal distance. 

Integration with respect to z and subsequent siabstitution of 


w == (r + C^) results in 


R 

Sy(y) = I [Y > y] = 1 - J ° [ ${(a - Hn z^)/a} + 


+ exp [ 6 {m^ + l/Cg [)in - a + $7(202) ]}] 

Q 

X ${(jln z^ + a a/Cg -a )/a } 3 fjj(r) dr 
R 

+ / ° exp [ ^{mp+l/G^ [ HQ + 0 $7(202) -a]} ]fg(r) dr; y > 0 

o 

... (6.11) 

where 

a = jin y + 0 ^ jin(r + 0 ^) -C (6.12) 
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and fg(') is the probability density of E. letails of fp(*) is given 


E" 


in Appendix A3. 8,1. Substitution of Hj. (6.il) in Eq. (5.24) yields 
average return period, of peak generalised seismic intensity. 

6.2. 1.2 ffiie probability of generalised seismic intensity for magnitude 

cutoff at both ends i. e. m e (m 1 

o 1 

Use of statistical independence and s^ibstitution of Bq. (6.5) 
into Eq. (6,2) results in 

w z 

Sy(y) = P[Y>y] = l- / f ^ 9{(m y - z - jin w )/0 } 


Wi 


f^(z) f^Xw) dz dw 


(6.13) 


viiere 


f^(z) = exp + (za - £n z)/Q^}'\ Yi^/izO^', z e (z^jZ^j (6.14) 


1 /K = 1 - exp [(m^ - m^) g] (6.15) 

Z 2 = exp (Cg m^) (6. l6) 

and w , w„ and z^ are defined in Bq. (6.8 - 6.10). 

12 1 

The following result is obtained by proceeding in similar lines 
as in Bq. (6.ll). 

E 

Sy(y) = P [Y > y] = 1+K J ° [ ^{ (a - kn z^/a) } 

2 ' 

-.${(a - 2,n z^)/a }+ exp [ ^{m^+l/C^ [ £a G^+a 6 /(eC^) -a ] > ] 

X {$[ (jluz + /b /C2-a)A 3 - ® [ (nn z^ + /C2 - « )/a] } ] fg(r) dr 


y > 0 


( 6 . 17 ) 
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where a is defined by Bi. (6.12). The probability density function, 
fj^(*) of focal distance, E, is given in Appendix A. 3.8. i. Substitution 
of Bi. (6.17) in Bq. (3.24) yields average return period, T , of peak 

y 

generalised seismic intensity. 

The Bi. (6. 17) is integrated nmerically for different attenu- 
ation laws and the results are presented in Pigs, (6,7 - 6,9). The 
in the attenuation law (Bj, 6.1) is taken to be Kf(o ,0,4225) for peak 
gromd acceleration and N(o, 0.5625) for peak ground velocity and 
displacement in accordance v/ith Esteva^^^\ 

6.3 Proposed Model for Scatter in Attenuation law - Model - 1 

The effect of scatter can also be accounted for by fitting 
envelope functions of the form 

Y = exp [Cg M - n (E + C^) ] (6. 18) 

(24) 

to the data used by Esteva and Boser^lueth \ 

let R* be the focal distance whesre both the upper envelope OB and the 
lower envelope OA, and the attmuation law v?ith the values of the constants 
as proposed by Esteva and Eosenblueth^^^^ (Section 3.6) meet at 0 as 
shown in iig. (6,3). ^ny line, OG, drawn throu^, 0, the point of 
intersection of the envelopes, OA and CB, is considered to be equally 
likely. Eie ratio of the shaded area (Pig, 6.3) to the total area 
bounded by OA and OB gives the probability that the actual generalised 
intensity will be less than or eqx:®.! to the generalised intensity predicted 
by using OC viiich represent the attenuation law. 
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6.3.1 Ihe Choice of envelopes 

Two sets of envelopes are drawn and are marked A and B in 
Pigs. (6.1 and 6.2). The envelope marked A encompases all the data 
points whereas envelope marked B encompases most of the data points 
within focal distance of 300 Ems. These envelopes are drawn visually. 

The values of C^, and for these envelopes are given in Table 6,1. 
Table 6,2 contains values of parameters and for a few probability 
levels. The envelopes A and B for peak gromd displacement are derived 
from the corresponding envelopes of peak groind acceSration and velocity 
and satisfying approximate relation 

2 

ad/v ^ 1 (6. 19) 

2 

where a, v and d are in cny'sec , cm/sec and centimeters respectively. 

The results in term of.B [ f > y] for a few probability levels 
are shown in Tigs. (6,4 - 6.6) for peak gromd acceleration, velocity 
and displacement and are obtained following procedure discussed in 
Ojtapter 3. 

6,4 Proposed Model for Scatter in attenuation law - Model ~ 2 

The generalised attenuation law (Bi. 3. 10) is given by 

Y = Z/W (6.») 

Z = exp (Cg M) 

W = (R + C^)^^ 


where 


( 6 . 21 ) 

( 6 . 22 ) 


Parameters of Bavelopes 
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Table 6.2 

Values of 0 and C_ for different probability levels for envelopes 
' ^ A and B 


Generalised 

Intaisity 

Bcivelope A 


Envelope B 


probability 0^ 

% 

probability 

S 


0.433 

2000.0 

2.0 

0.471 

2000. f 

2.000 


0.500 

7548.5 

2.^20 

0.500 

2760.3 

2-050 

Acceleration 

0.750 

1099.0 X 10^ 

2.95 

0.750 

42682 

2-475 


0.900 

2182.1 X 10^ 

3.40 

0-900 

2207. 1 X 10^ 

2.730 


0.950 

5908.7 X 10^ 

3.55 

0.950 

3816.7 X 10^ 

2.815 


0.316 

l6.0 

1-70 

0.294 

16.0 

1-700 


0.500 

2241.0 

2.40 

0.500 

161.2 

2.050 

Vel ocity 

0.750 

1833.4 X 10^ 

3-35 

0.750 

2662.9 

2.475 


0.900 

1025.6 X 10^ 

3.92 

0.900 

14329.0 

2.730 


0.950 

3922.2 X 10^ 

4.11 

0.950 

25109.0 

2.815 


0.283 

7-0 

1.60 

0.235 

7.0 

1-600 


0.500 

4614.2 

2.60 

0.500 

127.3 

2.050 

Displacement 

0.750 

8055.2 X 10^ 

3.75 

0,750 

1970.4 

2.475 


0.900 

7100.2 X 10^ 

4.44 

0.900 

10195.0 

2.730 


0.950 

3159.8 X 10 

4.67 

0.950 

17634.0 

2.815 


Note : Eie first row of each generalised intensity corresponds to 

h2) 

' ‘ value of parameters and C^. 


Esteva* s^ 
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lo account for the scatter in the seisBiic data and C_ are 

1 5 

assumed to be random variables in this section. Ihe range of these 

variables is determined by the upper and lower envelopes discussed in 

Section 6.3.1. 3!he constants and are empirical and are taken as 

in Section 3.6, Ihe random variables C, , C^, M and R are assumed to 

1 3 

be mutually independent, 23ae random variable is assumed to be 
log - Simpson i.e. 

X = iln 0^ (6.23) 

is a Sirapson distributed random variable- fhe probability density 
function of X is given by 

2(x - xp/{(x^ - x^) (xg - x^)}; in b^ = x^ _< x ^x^ =1^ b^ 

f^(x)=2(x - 3!^)/{(xj - x^ - x^)}; x^ £x ^x^ = Jin b^ 

0 ; otherwise (6.24) 

Eie random variable is assumed to be Simpson. Ihe 
corresponding probability density function is given by 

2(y - a^)/{(a^ - a^) (a^ - a^) } ; a^ £ y 

(y) = 2(y - a^)/{(aj - a^) (a^ - a^) } 5 a^ £y £a^ (6.25) 

3 

0 . ; otherwise 

fhe values of the parameters b^., b^ and b^ in Bj, (6.24) 
and that of a , a and a in Sj. (6.25) are given in feble 6.1, 

1 4 - J 

The probability distribution function of Z{'Sq_, 6,2l) is 
obtained in Appendix A6. 4. 1. ®ie details of probability density and 
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distribution fmction of w(Bq.. 6.22) are given in Appendix A6.4.2. Using 
the distribution laws of Z and W the probability of generalised intensity 
is obtained in Appendix A6.4.3. Kie results for peak grouad acceleration, 
velocity and displacement are obtained nunerically in terms of P[Y > y] 
for both envelopes and are marked as A and B corresponding to previous 
Section, These results are drown in ilgs, (6.7 - 6,9). 

6.5 Results and liscussion 

The resiiLT;s in Rigs. (6.4 - 6.9) are for magnitude cutoff at 
both level i, e. m e (5,9] and focal depth lognormally distributed with maxi- 
mxna focal depth, = 600 Kias. The earthquake source is a volume 
source with arc length b = 150 Kms. Comparison of Gaussian model and 
Model - 2 (Rigs,6,7 - 6.9) shows that Model - 2 predicts higher generalised 
intensity. The Gaussian Model and Model - 1 (envelope A and 50^ chance 
of actural peak acceleration to be less than or equal to the predicted 
peak acceleration) for same return period give comparable peak acceleration. 
However, for peak velocity and displacement, the Model - i with daanee less 
than 50^ is comparable witti the Gaussian model for same return period. 

The Model - 2 and Gaussian model predict average generalised intensity 
when scatter in the data are accounted for. However, for Model - 1 
Specific probability can be specified for p3?edicted intensity. This 
model mi#it be more useful for important structures (e. g. nuclear 
faciliti^). 

6.6 incorporation of the effect of Scatter on Seismic risk 

The steps to incorporate the effect of scatter graphically cn 
the sM.smic risk analysis is given in the following : 













EFFECT OF MODEL-2 AND GAUSSIAN MODEL ON THE 
EXCEEDANCE PROBABILITY OF PE AK DISPLACEMENT , 
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1. Determine the generalised intensity for a T-year return period from 
the contour map given in Chapter 4. 

2. locate the point G on the curve marked Esteva's in iigs. (6»4 - 6.9). 
Choose a model of scatter and project horizontally to the curve 
corresponding to the curve of the chosoa model to locate point D, 

3. Iroject vertically frcm D to the abscissa to get the generalised 
intensity incorporating scatter in the attenuation law. 

Example 

2 

let at a location peak acceleration is 27.5 cm/sec for iOO~year 

return period. The peak acceleration for lOO-year return period 

2 

is obtained following above steps as 37.5 cm/see for Gaussian model, 

as 125 cn/sec for Model — 1 with envelope B and 95/^ chance of 

probability of actual acceleration to be less than or equal to the 

2 

predicted value and as 805 cnv^sec for Model - 2 with the parameters 
corresponding to the envelope B» 

Ihe extrapolation for return period and for specified exceedance 
prctoability may also be incorporated subsequently following procedure 
reported in Chapter 4. 

As explained above, the generalised intensity incorporating the 
effect of scatter can be easily obtained from contour maps given in 
Chapter 4. This simple procedure is only possible for sources - 
homogeneous in tectonic feature. However, for ncnhomogeneous source 
fresh calculatiaa is necessary to include the effect of scatter on 


seismic risk analysis 



CHAPTER 7 


SEiasiIG ZONIH&'MAPS op INDIA 
7.1 latroduction 

Preparation of the seismic zoning maps of region involves a 
synthesis of the seismic history, geological features and an acceptable 
balance between cost and risk. The previous chapters contain tte 
following iofoimatiOE of relevance to seismic zoning 

1. Seismic data of the Indian penins\Jla for the period 1917 to 1972 
and its analysis of special relevance to seismic zoning is the 
seianicity map, Pig. (2.5 )» giving the energy contours. 

2. Contour maps for 100-year return period of the peak ground accelera- 
tion, velocity and displacement assuming the Indian peninsila to 

be homogeneous and nmhanogeneous. Pigs. (4.4 - 4,6 and 5.3 - 5.5), 
in tectonic features. 

3. Haps of Indian peninsula giving idealized tectonic features. Pig, 
(5.2), and epicentre of earthquakes. Pig. (5.l). 

4. A set of graphs throt^ which the generalised intensities can be 
extrapolated for i) upper limit on magnitude, Pigs. (3. 11 -3.13)> 
ii) different return periods Pigs. (4.7 -4.9), iii) specified 
exceedcince probability for a given service life. Pigs, (4.7 - 4,9), 
and iv) scatter in attenuation laws, Pigs. (6.4 - 6,9). 

(Hie seismic risk analysis perfctrmed in the thesis does not take 
into acccunt earthquakes prior to 1917. list of these eart]::qijak», some 
of which are quite large in size, is placed in Appendix A7, 1. 
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IS : 1895 - 1975} Indian Standard Recomaendation for Earthquake 
Resistant Design of structures^ , contains a zoning map of India. 

The code also contains maps showing epicentres and sizes of past 
earthqijakes , principal tectonic features and lithological groups, 

7.2 Seismic Zoning Maps 

In this chapter an attempt is made to prepare a seismic zoning 
map of India by integrating the results of the investigations presented 
in this thesis and the information available in the code IS : 1895 - 1975. 
Three maps have been prepared giving the values of the peak ground 
acceleration, velocity and displacement for a 100 year return period, 
Stepswise procedure used in prepairing the map is explained below t 

1, Contour map of five levels of peak ground acceleration was 
prepared on the basis of nonhomogeneous tectonic features. The 
acceleration values were modifed to include the effect of scatter using 
the Gaussian Error model and corrected maximum magnitude. The map is 
shown in Tig. (7.l). It divides the country into five zones with desiga 
acceleration for vinlform seismic risk. 

2, The lateral force coefficients given in the seismic zoning 
map of IS ; 1895 - 1975 are divifed by 2.6 to give the peak ground 

(93) 

acceleration for 5% damping'’ . , 

3, The map prepared under Step - 1 is compared with jwning map of 
code. The bomdary of the zones are carefully compared and isherever 
major deviations are noticed, geological and epicentral evidaices 


are examined 




Scale' ' 


miles 


FIG-7i PEAK ACCeLERATION CONTOURS 
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4. integrating all the above informat ion8,the seismic zoning map 

of India is drawn, (Slg. 7.2), which represents the best judgement of 
the author. 

Salient features and the deviations between the code map and 
the map prepared by author are : 

In the Iteccan plateau (below 12° latitiade) evidence of earthquake 
origin exists in the form of past earthquakes, This evidence is 
accepted and this region is upgraded. In the region of grid point 70 
earthqu^es of magnitude greater than 6 have occurred in the past and 
consequently this region is upgraded locally. In the Koyna region 
(grid point 84) the seismically active region is enlarged as compared 
to code map. Upgradation of Gondowana rift vinit (grid point 88 and 123) 
is the result of general increment in gradation in the Deccan plateau, 
!Ehe zone boundaries in this region are modifed in conformity with the 
tectonic features. The presence of inactive seismic fault number 7 in 
the West coast and Narmada-Tapti unit is largely ignored in the 
construction of zoning map. In the central 2 !one a narrow higher zone 
is accepted in the region of fault numbers 11 and 12. In the Himalayan 
Tectonic zone evidence of historic eartlmiuake in Srinagar, Ifendi and 
Babhange (grid points 229, 214 and 171 respectively) is taken into 
aocomt by upgrading these regions locally. In the Assam region, 
adjoining Bangladesii and Burma, a lower zone exists inconformity to the 
contour maps previously obtained. Bie region aromd Eiiuj in Gujrat 
(grid point 130) is also upgraded inconformity with historic data. 
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The seismic zoning maps for peak ground velocity and displacement 
have been prepared following a similar procedure. The zoning maps and 
five intensity level contour maps are shown in Fi^. (7.3 -7.6). 3h 
these cases, however, code does not provide a map, 

7.3 Discussion of the Seismic Zoning Maps 

Hie major deviations in the boundaries of the zones between the 
code map and that prepared by the author are along the West coast, in 
Central India and in the Assam, region. The code map is based largely 
on the largest past earthq.xiakes whereas maps based on seismic risk 
analysis are weighted towards large number of smaller earthquakes close 
to the site^^^^ . 

Table 7. 1 shows a comparison of the peak ground acceleration in 
different zones based on two maps. It is seen from the table that code 
provision in different zones are not uniform frcm the view point of 
acceptable risk and is different by as much as a factor of 4,5. Kie 
intensities indicated in the author’s map are for lOO-year return period 
and take into account the effects of scatter. 

Table 7.2 shows the comparison for ii major cities in different 
parts of the country. 

Based on the seismic rii^ analysis of the present work it is 
recommended that the late3ral force coefficient (a^) in the IS : 1893 - 1975 
for zones III, I¥ and V should be incr'eas©! to correspond to 25 - 30 
year return period, Hiese values are indicated below. 




Scale 


FIG-7-3 PEAK VELOCITY (cm/seel CONTOURS 



Islands.: 







FIG-'7-5: PEAK"'©. 
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Table 7- 1 

Ccanparison of IS : 1893-1975 v/ith Present Work 


IS : 1893 - 1975 Present v/ork 


Zone 

lateral 

force 

coefficient 

Accelera- 
tion in g 

(ay 2-6) 

Typical 
retiirn 
period 
in years 

fo exceedance 
probability 
for 100- 
years 

Accelera- 
tion in g 
for 100- 
year 
return 
period 

5^ exceedance 
acceleration 
in g for 
100-years 

Y 

0,08 

0.0308 

10 

99.99 

0.095 

0.270 

TV 

0.05 

0.0192 

10 

99.99 

0.061 

0. 178 

III 

0.04 

0.0154 

15 

99.90 

0.041 

0. 141 

II 

0.02 

0.0077 

25 

95.00 

0.024 

0.087 

I 

0.01 

0.0038 

45 

80.00 

0.010 

0.046 
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Table 7.2 

Peak gromd Acceleration in g for Different Cities 


City 

IS ; 

1895 - 1975 



Present work 


Zone 

No. 

Acceleration 

Return 

period 

in 

years 

Zone 

No. 

Acceleration 
for 100 -year 
return 
period 

5^ exceedance 
probability 
Acceleration 
for 100-years 

Bombay 

Ill 

0.0154 

60 

II 

0.024 

0.087 

Calcutta 

III 

0.0154 

60 

II 

0.024 

0.087 

Chandigarh 

I? 

0.0192 

22 

III 

0.041 

0.141 

Coimbatore 

Ill 

0.0154 

60 

II 

0.024 

0.087 

Delhi 

lY 

0.0192 

22 

III 

0.041 

0.141 

Gauiiati 

V 

0.0308 

10 

V 

0.095 

0.270 

Hydi rlad 

I 

0.00308 

45 

I 

0.010 

0.046 

Kanpur 

III 

0.0154 

60 

II 

0.024 

0.087 

Ifedras 

II 

0.0077 

25 

II 

0.024 

0.087 

Norora 

lY 

0.0192 

22 

III 

0.041 

0.141 

Srina^r 

Y 

0.0308 

23 

lY 

0.061 

0.178 
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Zone 

Present Code 

Recommended Code 


provision for 

provision for 


a 

a 


0 

o 

Y 

0.08 

0.13 

IV 

0.05 

0.08 

III 

0.04 

0-07 


Qlie generalised intensity in the author's map can be extrapolated 
for different return periods and specified exceedance probability for a 
gi-yen service life. For example, at Norora, the peak ground acceleration 
for an exceedance probability of 5/o over a 200 year period works out to 
be 0.257g. 

The generalised intensity values in the seismic zoning map are 
average values for the zones. These values may be exceeded in the close 
proximity of the epicentre of a large earthquake. The generalised intmaity 
values are for firm soil and shall be modified suitably for local soil 
condition, 

Ito.jor difference between acceleration and velocity or displacement 
maps is the two higher zones in region near grid point 244 as is shown 
in Jigs, (7.2, 7.4 and 7.6). In case of velocity map oiiLy zone IV is 
in the region near grid point 199 and only zone II is present near grid 
point 84. Both acceleiation and displacement map contain two zmes near 
these points. Marginal differences exist between boundaries of various 
zones in all the three maps. 



CHA-ITziR 8 


■SmiriAHr os- RESUITS and CONCIUSIONS /J'D SUGG}53TI0N3 roe FUREIER '''OHK 
8,1 Results and Conclusions 

Ihe statistical amlysis of seismic data and the seismic risk 
analysis of Indian perdnsiala have been presented in the preceding 
chapters. The results and conclusions of the investigation can be 
suamail sed as : 

1. Statistical analysis of seismic data for the period 1917 to 1972 
indicates that i) arrival process of earthquakes is not Poisson; (ii) the 
data are space-time correlated, (iii) the spatial pattern is non- 
stationary (iv) Markov chain of the Sani-Markov model giving transition 
of earthquakes from one area to another fits the data. However, final 
inter-arrival distribution of earthquakes obtained from. Semi-Markov model 
does not fit the data. Semi-Markov model for the earthquake process can 
be used for seismic risk analysis provided sufficient data are available, 

2, Effect of focal depth treated as a random variable on the seismic 
risk is significant. A ndassd lognoimal distribution of maximum focal 
depth, = 600 Kias, , gives the best fit to the data. Dogaormal 
distribution wilii 600 iOas. maximum focal depth gives comparable risk 
whereas same is true for uiiform distribution with 15i Sns. cutoff focal 
depth. 

5. Contour maps of equal generalised intensity based on seiaaic risk 
analysis tising homogeneous and nonhcsaogeneous source model show marginal 
deviations. The tectonic features, such as faults, cause local variation 
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in the sbape of the contours but the value of the generalised intensity 
is not changed significantly. 

In view of the above, it is concluded that the seismic zoning 
map should be drawn on the basis of seismic risk analysis with sources 
identified on the basis of tectonic features. Efowever, extrapolation 
for return period and exceedance probability may be made on the basis 
of a modular earthquake source. Curves of such extrapolation can be 
provided in the code. 

4. ffiie scatter in attenuation law has a significant effect cm seismic 
risk. Model - i proposed in this thesis to account for the scatter agree 

(42) 

closely to the Gaussian model due to Esteva for some probability lev< 
idiereas proposed Model - 2 predicts significantly higher generalised 
intensity in comparison to the Gaussian model. The effect of scatter 
should be included while computing the seismic lisk. 

5. A generally accepted approach to earthquake resistant design 
requires the structure to ranain elastic for 'small’ earthqiiakes, to 
xindergo repairable limited plastic deformations for 'medim' earthqu&kes 
and large plastic deformation, without collapse, for 'large' earthquakes 
The criteria for choosing small, medium and large earthquakes are not 
laid. It is proposed that the size of the earthquakes be linked with 
the return period; small, medixma and large earthquake, e.g. , may 
correspond to return periods of 30, 100 and 200 respectively, for 

a site in zone III the desi^ grotiM acceleration is given below* 
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Type of 

Ret\irn period Design peak 

Exceedance 

earthquake 

years ground 

Probability for 


acceleration 

100-year return 


in 

£T 

period 

Small 

o 

30 0.022 

0.50 

Meiium 

100 0.041 

0.67 

large 

200 0.066 

0.85 

The small. 

medium and large earthquakes may 

also be specified 


in teims of exceedance probability for lOO-years as indicated in the 
last column of the above table. 


6. Ihe lateral force coefficients specified in IS ; 1893 - 1975 do not 
correspond to idaitical acceptable risk in each aone. The ratio of 
risk varies as much as 4.5. It is recommended that lateral force 
coefficients in ajne III, IV and V be increased to 0.07, 0.08 and 0,13 
to correspond to 25 - 30 year retijim period. 

7. Since ground acceleration, velocity and displacement govern the 
design of different types of structures, it is meaningfxil to draw 
separate zoning maps for mch of these. These maps show significant 
differoaces in northern and western region. 

8. Ohe proposed zoning maps provide design generalised intensities 
for eqioal acceptable risk and suggests changes in the bomdaries of 
code map. She desigi intensities can be modified for different return 
pahLod and exceedance probabilities, 

8,2 aiggestione for Further Work 

1, The stochastic model of earthiuake ocexurence us^ in this thesis 
and used other investi^tors ne^ects space-time correlation, magaitude 
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focal depth correlation and hare assimed temporal nature to be Bsisson. 
Efforts should be made to incorporate the space— time and magnitude focal 
depth correlation in future work. A simulation approach incorporating 
the above corr^ations and assuming earthquake process is a stocbastic 
field process may be attempted. 

2. The main difficulty in the seismic risk analysis at this stage is 

the paucity of data which makes the estimation of parameters very difficult- 
While efforts to collect more data continue, attempts should be made to 
improve the mechanistic models of earthquake occurrence. 

3. At present the attenuation laws are incorporated at the level of 
local seismicity. Attempt ^oiid be made to establish these laws for 
major geological acmes so that they can be used at the regional level. 

In United States different attenuation laws are used for Eastern and 
Western regions, 

4. It is well knovaa that seismic data reported by different agencies have 
a measure of uncertainty regarding size and location of earthquakes, fhe 
reliability associated with such meertainties my be established and 
incorporated in the seismic risk analysis. 

5. She seismic zoning maps presented in this thesis may be refined by 
a more careful and informed identification of earthquake sources and 
the seismic activity associated with them, use of focal depth and upper 
limit of magnitude as regiocaal variables and use of 'judgement’ factor 
by a grorp of sei biologists and engineers. After such a map is prepared, 
the preset zoning mp in IS : 1893 - 1975 may be reviewed to incorporate 
some of the features suggested in this thesis. 
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Appendix A2 -4 ’1 


ilAXMmi LIKELIHOOD ESTBIATION FOR FOCAL DEPTH 
a) Truncated Lognormal DistriTxition 

Let a random sample of size n is drawn from truncated 
lognormal density: 

fjj(x) -= exp [-0.5()!,n x-e^ )^/e2 ] / [ /2Tr62 x /Sg }] 

; X e (0,H^3 (A2 4.1.1) 

where 9^ and 6^ are parameters ofthe distribution ¥/ith specified 
and $(,) is the probability distribution function of n(0,i). 

The logarithm of likelihood function excliding a constant 
term is 

_n [ 1/2 { jin 62+ (s - 2 6^ x + 8^)702 } + An $ ] (A244 4) 

where 

X = ( y jin x,)/„ (A2 4.1 .5) 

s = ^An x^)^/n (A2 4.I4) 

and is urderstood to be a function of 9^ and To find the 

location of maximum of L , following is computed 

«= -n [ (e^- x)/a2 + c/i>] (424.1.5) 

— = -n [{ 1-(s - 2e^x + e^/Og) /(SSg) + (An H^-6^)C/(202$)] 

(A2 4 .1 .6 ) 
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whei« 

C = -exp [-0.5(iin H^-6^f/e2] ^ (^24-1 -7) 

Solution for 9^ and should be obtained from these 

derivatives (Eq,A2 4.1.5 & A4 4.1.6) by equating it to zero. These 

- ( 94 ) 

equations are solved by the method of parametric differentiation . 
Assume X is a parameter such that X e [0,l] and at X = 0 the 
values of 0 == 0° . The auxili^iy equations from the above 

derivatives are 


alVse^ = (l-^) 00 I i=1»2 (A2 4 4.8) 

The above equations have characteristics that at X = 0, 0= 0° 

ani at X = 1 is the solution of derivatives (Eq. A2 4*'l*5 
X2 4,1,6) equal to zero. 


Differentiating Eq, (A2 4»1*8) with respect to X 
the chain rule yields 

\ d8./dxj \ shVae I 

; > I 

[ H ] 1 ! = - ^ ; 

l (5i' /“2 ie. e° 

where [ H ] is a mtrix whose elements are 


arsi using 


(i2 4.1 »9) 


h^^= aV/ae^ = -n [i/e^ + {[ (an H^-e^Veg-c/^jc }/$] (124.1.10) 


-1 -(a n 9^ )c/^.] C/(2 Bg )3 } /<S ] 


(124.1 .11 ) 
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ii22= » -ti[{(s-2 0^x + e^Ve^- 1/2} /e^ 

+ {[(iln H^-9^)V02-5-(to H^-e^)cA] Un H^-e^ )C/(4e2) }A J 

(124.1 .12) 

Equation (12 4.1.9) represents a set of two simultaneous 

oi^dinary differential equati ore which can he solved numerically fox 

d0^/dX ; i=1 ,2 . resulting ordinary differential equations 

o 

with initial condition Q= © at X = 0 is numerically integrated 
upto X = 1 to give the solution 0 . The Gill variation of 

( 95 ) 

Runge-Kutta formulae is ussi in the nuBerical integration scheme. 

Step length for integration will he governed by the acciiracy required 
and the speed of calculation. It is to he noticed that the elements 
of [ H ] matrix will have to he con^juted at every step of numerical 

A 

integration process . The maximum likelihood estimator 0 for the 
parameter 0 from sample size n, for large n, is approximately 
distributed as multivariate nomnal distribution with mean 0 and with 
covariance 

« = C-ErH]]"’ =[ -[H] §]"^ (il24.1.15) 

h) Truixsated Mixed Lognormal Distribution 

Let a random sample of size n is drawn from truncated mixed 

lognormal density function 

2 ^ 
fg(x) = I p^ exp[-0.5(£n x- 0 ^^)‘'/e 2 i] 

i=1 

X ${(£n H^-0^.)//9~}]5 X € (0,H^j ( 124 . 1 . 14 ) 
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with 

^ p. = 1 (12 4*1.15) 

e^2 >611 (124.1.16) 

and where 9^^ and ^2i* parameters of the distribution 

with specified ^ (-) is tte probability distribution of N(0,l), 

The logarithm of likelihood function is 

n 

L* = I j»n f (x ) (124.1.17) 

j=1 

The method of maximum likelihood consists of solving for those 

* 

estimate 0 of 0 , which maximizes L satisfying conditions 

givm in Eqs- (12 4.1*15 ~ 124*1.16). By transforming the 
variables as below; 


. 2 

= sin y^ 

(124.1.18) 

P = 1 - p 

2 n 

(124.1.19) 

r. 2 

®ir ^2 

(12 4 .1 .20) 

8 ,2- 8^, + (itn H^-e„) 

(124.1.21) 

®21- h 

(124.1 .22) 

. 2 
®22“ ^5 

(124-1.25) 


and changing objective function L to 
L =* -L* 


(12 4.1 .2 
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the stated maximization problem is changed to unconstrained minimization 
problem. The minimization can be done by any standard scheme to 

/s (%) 

estimate the parameters G . However, Pov/el's method of unconstrained 
minimization has been adopted in this work. The covariances of the 
estimate for large sample may be obtained approximately. 
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Appendix A2,6.1 

IIAXE-IUII-LIKSLIHOOD ESTi; lATIGK OP Hti»ERdXP0NI21TI\L 
DISTPIBI3TI01I 


Let an earthquake process has been observed for time t^, starting 
from an event at t = 0, Let n earthquakes have occurired at 
t^,.....,t . Therefore interarrival times are 

x^ = t^ (A2 .6 .1 .1 ) 

X. = t.- t. . ; i=2,...,n (A2.6.1^) 

and they are samples from hyperexponential distribution v;ith probabi- 


lity density function 


^ I p. exp (-v^t) j t ^ 0 

i^l 

with 


Pi = 1 


^2 ^ ’^l 


(A2 .6.1.3) 


(ii2 ,6.1.4) 

(A2 #6 .1 *5 ) 


The log likelihood is given by 

L = (n+l) An p^+n An P^- P^ t^+ An {l+P 2 /p^ exp[ (p^-P 2 )u] ? 

+ ^ An {I+P 2 P 2 exp[ (u^-V 2 )xi]/(p^V^) } (A2.6 .1.6) 

where 

u = tQ - t^ (12.6.1.7) 

The Eq , (A2 .6 .1 .6 ) is maximized subject to the coalitions 
given by Eqs. (A2 .6 .1 4 and 12 .6 .1.5). Following transformation 
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reduces the constrained maximization problem to an unconstrained 
minimization problem: 


minimize 

L 

= -L* 

with 

P 

1 

= sin^fx. ) 
1 


P2 

= 1 - P^ 



II 


^2 



Here x^; i=1,2,5 are slack variables and L is a function 
of slack variables. The solution of the minimization problem can 
be obtained by any standard method. However, Pavel's method^ 


is used in this work 
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PRCBASILI'TY DISTRBUTIOH PUldlOK IIS PROBASILII’Y 
DENSITY PDICTIOH CP PCCAL’ DISTA!SS,3, FOR YOLDlvE 
SOURCE AT SITE 

a) Focal depth, H,i 8 unifoxnO .7 distiributed in (0,H^ J 


It is assumed in this section that an eartbg,uake is equally 
likely to occur anywhere in the Tolume source Pig. (A3 .8.1 .1) which . 
influence the peak generalized intensity of the site. The 
probability that the focal distance, R, will be less than or equal to r is 
equal to the ratio ef the volume of the earth upto r to the 
volume of the earth upte maxifflum focal distance R , i.e. 

P^(r) ^ P[B<rj= T(r)/7 (A5.8.I.I) 

K — 


where 


a cos V 2 ^22 

Y =• ir / (y tan v) dy + ir / (a -y )dy 




(a-H )oos V 


a cos V 

a 


(j tan v) tly t 


ta^H^ )dos V 


(a^-y^ )dy 


ard 


2 itG 


C »= (I-00S u) {a^-(ar-H )^}/5 


(A3 .8 ,1.2) 


(15 .8. 1.3) 


Par e [ 0, a sinuj, the toIub®, Y(r), upto the focal distance r 
is given by from Pig. (A5.8.I.1) as 


Y(r) = ir / (2 ay - y^ )dy + J 'fr(r^-y^)dy 

0 ■ y« 


(A5.8.14) 
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# 



MAXIMUM FOCAL DEPTH 
RADIUS OF EARTH 
maximum focal JtS FAI'CE 

la* (a - Hof- 2a(Q -Ho) Cos v 
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■where 

« a {1 - cos (r/a)} (A5^,1.5) 

The integrand of Eq. ( A3 .8.1. 4) results in 

V(r) = 2 k {rV? - 

Suhsti'tuting this and Eq. (A5.8.I.2) intoEq. (A3.8.I.I) yields 
probability distribution function of R ani is given by 

Eg(r) = {rV? - y^i^^ - ayQ)/2}/C ; C<_r;<a sin v (A3 .3 .1 .6) 


Differentiating with respect to r, the probability density 
fiinction of R is obtained as 


f^^(r) = {r(r--y^) + 1/2 {2ay^-r^} s3-n (r/a)}/C 


j 0 < r < a sin v (A5.8.I.7) 


The volume, V(r), for r e [ a sin v, 2a sin (v/2)j is 
obtained from Fig. (A5.‘^.1»1) as 



where 


and 


^ 


^2 “ '^2 


u^ = a sin V 


2 2 2 

u^ “ COS v/r - a sin v 


(A5,8.1.8) 


(A5.8.I.9) 


(A3 .8.1 .10) 
(A5.8.I ,11) 
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After integration and simplification it can "be shorn that 

V(r) >=> ir[2r^+ JJq {ay^-r^}- 4u2(r^-a^sin^ '^)]/5 (A5,8»1t12) 

Probability distribution function of d is obtained from 
Eq. (A5,a,l,2) and Eq, (A5.B,1-12) as 

Pg(r) - 2u2(r^- a^sin^v)}/5 + y^Cay^- r^)/2 ] /C 

; a sin V < m ^ 2a sin (v/2) (A5»8»1*15) 

TheDcefore, probability density function of S is 

fj^(r) = [r {r[l - l/2 sin (r/a)3 - Zu^)] 

+ {a sin (r/a) - r} y^ ]/G 5 a sin V ^^^2a sin (’^/z) 

(A3 .8 .1.14) 

The TolunBj '^(r), for r e ( 2a sin (v/z), ] is given by 

Pig. (A3.a.1 ,1 ) as 

y^ ^2 

Y(r) = trr J (2ay - y^)ay + tan^ v / (ar-y) dy 
o yj 

+ / ir(r^-y^)dyj 

^2 

o tr£2r - ay^ (2 + cos v) cos v 

- C2y2+ (r^~a^ sin^ v)]/3 (A 3 , 8 . 1 . 16 ) 

where 

y^ *=■ a(l — cos v) (A3.8.1 .17 j 

Hence, the probability distribution function of R is obtain®! as 
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Fj^(r) =■ [ {ay^ (2 + cos v)cos v}/2 - (7^+ Ti^/2) 

(r^- a^sin^ v)] /(5C) (5'»8*1.18) 

Consequently, the probability density of R reduces to 

- r [ r - y^ ] /C (A5.3,m9) 

Prom Fig, (A3.8,1.1 )j the TOlume, T(r), for r e J 

wlaere 

R^= /a^+ (ar-H^)^^ 2a(a-H^) cos v (A5,8.1^0) 

and y^ = - (r^- I]^)/(2a) (A3.8.a.2l) 

is given by 

^3 ^2 

¥(r) = Tr[ J (2ay - y^ )dy + tan^ v / (a~y) dy 

o 

Ho y^ 

+ f + / [ r^-y^" {(a-H^)^-(a-yf }]cly ] 

Above expression after integration and sii!5)lification reduces to 

y(r) = ■rr[{4a Eo(3r^~I^) + 3(r^- H^)^} /(4a) 

“ ay^ (2+COS v)cosv- (272 + )( 3:^^-a^sin“ v)] /? (A3 *8 .1^2) 

Substituting this and Eq. (A5»8.1*2) in Eq* (A5*8,1 ,i) yields 
the probability distribution of R as 

PR(r) = [ {Ho(5r^- 4) + 5(r^~4)V(^a) 

- ay^ (2 + cos v) cc® v}/2 

- (y2+ u^/2)(r^-a^sin^ v)] /(3C) 5 Ho<_r<So (A5.8.1.25) 
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Differentiating above with respectto r,tte probability density 
function of R is obtained as 

f^(r) - r[H^- (r^- 4)/(2a) - /C (A5.8,1.24) 

b) Probability distribution fuixtion of Focal Depth, H, is 
mixed lognormal in (0,H^ j 

The probability distribution function of focal distance , R, for 
volume earthquake source Fig.(A5»‘^*1 *1 ) is derived based on statistical 
independeme of azimuth, A , latitude, 0 , and focal depth, H. The 

joint probability density function of A , 3 and H is given by 


where 


^HAG 

= ^A^"' 

(A 5 . 8 .I. 25 ) 


; X e (0,u] 

(A 5 . 8 .I ^6) 

fe(e) = V(2Tr) 

» 0 e (0>2 it i 

(A3J3.1 ^7) 

Vh) . 1 p, 
1=1 

g.(h)/l .(H ) 

®i^ 01 o' 

(A5.8.1 ^ 8 ) 

I Pi “ "I 


(A3.8.1^9) 


gj.(h) = exp [-(tn h h) ; i=1 ,2 

(A5.8,1.30) 

I^i(Ho) => exp {(e^.- 1/2 me^i) m) $[(ltn 

S m 0 (A5.8#1 . 31 ) 

and 4(«) is the probability distribution function of H(0, 1 ), 

The probability of the event [R < r J is the ratio of the 
volume of earth upto focal distance, r, to the volume q£ eaaHih upto maximm 



185 


laaxim'UHi focal distance, , where earthquakes can potentially occur^, 
In mathaiaatical notation it can be written as 


P [ S < r ] = P_^(r) = Y(r)A(S^) (A5 .8 .1 .52) 

To eraluate volume upto maxiaum focal distance, R^, a new random 
variable Y = a-H is defined, where, a, is radius of earth. It can he 
easily shovm that probability density function of Y is given by 


f^Cy) - (a-y) ; y e [a^H^,a) (A5.8.a.55) 

The joint probability density function of A , 9 and Y may be 
obtained from Eq. (a. 3»8-1.55) and Bq, (A3.8.1.25) ard is given by 

^Yag ^ “ ^Y^^^ ^A^^^ ^0^®^ (A5.8.1.34) 

The volume upto maxinnim focal distaiace, S , is obtained using 


Bq- (A5.B.1.54) as 


2t! 


^(Rq) = / y / sin X J fy <;.’-,XrS)d8 dx dy 


a^H o 

o 


2it 


l/(2Trv) J y f„(a-y) | sin X J do dX dy 


a-H 


Defining 




- 2a / [v ] (13.6.1.35) 


X’ "arc cos {(a^+ y^ -r^)/(2 ay)} 


(A5.S.1.36) 
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the volume ‘V(r), from Pig, (A5.S«1.2) for re[ 0,a sin v] can be 
written as 

a, - V 2w 

■=. / y / sin X / fy.g (yjX , e)d0dXdy (A5.8.I.57) 
ar-r o o 

Substitution of Bq. (A5,8.1.36), Bq. (A5,8.1.26) and Eq. (A5,8.1.27) 
inEq, (A3*8,1*57) and completion of twofold integral yields 

a 

7(r) = l/(2av) / y [r -(ar-y) ] f^.(y) dy (A5.8.I.58} 

Br-X 

Substituting Eq. (A5,8.1.33), Eq. (A3.8.1.28) and Eq . (AJ.S.I.JO) 
InEq. (il5 •8.1*58) and integrating results in 

y(r) - I ~ I^-Cr)} + {I^.(r) 

~ r-I^^(r)} /a] / fev (H^)} (A3-B.1.59) 

Substituting this and Eq. (A5.8.1,55) into Eq. (A5,8,1.52) 
yields the probability distribution funotion of focal distance, E, 

as 

Fj,(r) - I - Ijifc) + (Ijit) - 

2 

- r I^^(r)} /a] j 0 5^^ 5.a sin v (A3 ,8.1.40) 

wheare 

T* » 2(l^os v) ^ p J (H^)-2a I^i(H^)] 

(A5.8.I4I) 

The probability density function of focal distance, H, is obta- 
ined by differentiating Eq. (A5,8.1 i 40) vrith respect to r and 
is given "bj 
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fj^(r) = 2r^p^[l^^(r) - I^^(r)/a] /i* O^r^asinv (i. 3 ^. 1 . 42 ) 

For r e [a sin vj 2a sin (v/2)] , the -volume, V(r), from 

Fig. ( 5^-1 *3) is 

a X' 2tt 

7(2?) = / y / sin X / f„ (y, X , e )<ie dX dy 

yg o ° 

+ J y J sin X J fyj^g (y» X , e )de dX dy 

^2 ° ° 

+ / y / sin X / fy^Q (y» X, e )de dX dy (A 5 . 8 .I 45 ) 

a-*r o o 


where 



2 2 2 

X' =■ arccos [ (a + y ~ r )/(2ay) ] 

(A5-8.144) 


,22 .2 

y . «=» a cos V / r --a sxn v 

(A5.8.145) 

and 

2 2 2 

y„ «= a cos V + / r -a sin v 

(15.8.146) 


Suhstitution of Eq,* (A3«8#1»55)y 29. (A5*8*1»28 )i Eq, (A3*8#1*50) 

and Eq, (A5.8.I.55) into Eq, (A3.8,1.43) ani completion of integration 
results in 

7(r) « I p^[r^I^^(r) - I^i^^) + { J/a 

+ GiCar-y^) - Gi(a-y-, ) ] / { 2 v } (iJ.S .1 47) 

where 

G.(x) = {r^ - 2a^(l-cos v)} 1 .(x) + {4a(l-cos \))-r^/a} I.j(x) 

X ox IX 


+ {2 cosv- 5 } ^x(^) + I^^(x)/a 


{A 5 . 3 .1 43) 
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Probability distribution function of R is obtained by substi- 
tuting Bq, (A5*8.1*47) and Eq.. (A5.8.1.55) into Eq. (A3.3.1.32) and 
is given by 

Fj,(x) = J pj [r^I„i(r) - + {Ij.t) - r==I,.(r) )/a 


+ G^(a-y2) - '» a sin v<^r-^2a sin(v/2) 

( 13 . 8 . 149 ) 

Differentiating Eq. (A3.8-149) with respect to r yields tl» 
probability density function of R as 

f^{x) . I pj2x a„i(x) - 

- [ + I^^Ca-yg)] /a} 

+ q.^(a-y^) + q^Ca-Fg) ] A’’* I a. sin v_^r^2a sin (v/2) 

(A3 .8.1 .50) 


where 


2 2 2 

q^(x) => s[ 2a (l-cos v) - r + x {r /a - 4a(l-cos v) } 


+ x^(3-2 cos v) - x^/a]g:i(x) 


and 


c/V 


2 2.2 


s « T/ y X - a sxn v 


• (A3 .8 .1.51) 


(A5.8.1.52) 


Prom Pig. (A3 . 8 . 1 . 4 ), the volume, ^(r), for re [2a sin(v/2),ii^2 
can be written as 


a 2 


^(r) = J y J sin X J (y» X » e ) dX dy 

T ., 


y o 0 


1 2 


2v 


+ I y / sin X / (yj X » 9)de dX dy 

a-r 0 o 



Simplification after substitution of Eq. (A5«0.1.35 )j (A5»8*'1*26), 

Eq, (A5,<5 .1 .2«l) and Eq. (A 5 .B.I. 54 ) into this and integration results in 

T(r) p. [ r^I^.(r) - (r) + {I^.(r) - r^I^ .(r) } /a 

- 6i(a-y^ ) ] / {2v ^ (A5.8.1 .53) 

Substituting Eq. (A3.8.1«53) and Eq. (3.8.1-35) into Eq. (A3.8.1.32) 
the probability distribution function of R is obtained as 

J?E(r) . I pj + ( I,.(r) - r^I,i(r) >/a 

- ; 2a sin (Vs) (A 3 . 8 .I. 54 ) 

The probability density function of R is obtained by differentiation 
of Eq. ( 15 . 8 . 1 . 54 ) with respect to r and is given by 

yr) = I p.[ 2r {I^.(r) - "C A> 

+ q^(a-y^)] AJ f 2a sin (v/ 2 ) ( 13 . 6 . 1 . 55 ) 

The Tolume V(r),, for r e ] from Fig. (A5-^S.1,5) is 

given hy 

^ 1 X * 2if 

'^i^) “ / J sin X / f„ (y, x , e)de dx dy 

&-H o o 

o 

a „ V 2ir , . 

4- / / I Bin X f (y» X , 6 m dX dy 

y^ 0 o 

Integrating this after sulratitution of Iki, ( 15 .8.1 . 54 ), Eq. (A 5 . 8 .I. 35 ), 
Eq. (15-8.1,28) and Eq. (A5. 8.1 - 26 ) yields 
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T(r) - y p. [x^I .(H ) + {L.(H (H ) } /a 

^ i. *^1 1- oi 0^ 2i 51^ o-^ lx' ' 


Gi(a-.y^)]/ {2v 


(A3 .8 .1 .56 ) 


Substituting Eq , (A 5 .S.I. 56 ) and Eq. (A3 *8 .1.35) into Eq. (A3.S.1.52) 
results in the probability distribution fuostion of R and is given by 

E„(r) = y p. [r^I .(h )-I„.(H ) + {I,.(H )-r^I. .(H )) /a 

R' A L 01 ' o' 2l' 33- 0 ' 1i' o' ‘ 


G^(a-yp] A* 1 ?-o 


(A 5 . 8 .I. 57 ) 


where 


R„ = / a^+ (a-H f- 2a(a-H ) cos v 

o ^ o ^ o 


(A 5 . 8 . 1 .58) 


differentiating Eq. (A 5 .S.I. 57 ) with respect to r, the 
prolDatility density function of E is obtained as 


+ 1 j^(a-y^ ) ] /''^ * 1 ^ ^ 


(45.8.1,59) 


c) Probability Distribution function of Focal Depth, H, is 
lognormal in (0 »Hq } 

The probability distribution and density function of focal 
distance, R, for volume earthquake source for the lognormal focal 

depth H can be obtained from the previous section by follosring 
substitution 

■ » 1 and P 2 = 0 

®1i ®1 ®2i = ®2 
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The probability distributioQ function, fg{. ), of the focal 
H, after above substitution is the following: 

Fjj(r) = [ r^I^(r)-l2(r)+ {I^(r ) }/a ] /V* 

» ^ 1 . ^ 

5g(r) = [ r‘^I^(r)-l2(r)+ {I^(r )-r^I^(r ) }/a 

+ G( a-y^ )-CJ( a-y ^ ) ] /¥* ;asin v<_r ^ 2a sin(v/2) 

- 6(a-y^ )] /V* ; 2a sin. (v/ 2 ) 

Pj^(r) = [ r^I^(H^)-l2(H^)+ {I^(H^)-r^I^(H^) }/a 

- &(a-y^).]/v* 

where 

= exp {(8.|-i/ 2 me2)n}4[ H^-e^-jaeg )/*^] ; m ^0 
v''= 2 (i-c08 v) [ ) + IgCH^) - 2a I^(H^) ] 

Gr(x) = {r^- 2a^( 1-cos v)} I^Cx) + {4a( 1-cos v) _ i/a r^} 

+ {2 cos V -3 } IgCx) + l/a IjCx) 

^ "■ g ■ 

y^ = a cos V - / r - a sin v 
1 

2 2 2 

y^ = a cos V - /r - a sin v 


distance , 

( A5 .8. 1.60) 

(A3.8.1.6i) 

(A3.8. 1.62) 

( A3. 8. 1.63) 

(A3. 8. 1.64) 
(A3.8.1.65) 
(x) 

(A5.8.1.66) 

(A3.8.1.67) 

(A3.8.1.68) 


=/ a^+ (a-H^)^- 2a(a-I^) cos \> 


(A5.8.1.69) 
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aad 1’(,) is the probability distribution function of N(o,l). 

The corresponding probability density function, fg(. ), is 
gi-ven by 

f^(r) = 2r [l^(r) - i/a I^(r)J /?* ; 0 ^ a sin v 

=[2r {IgCr) - I^(a-y^) + I^Ca-y^) 

- - I^(a-y^) + I^Ca-y^) J/a} 

+ q.(a-y^) + qCa-y^)] /V* j a sin V ^r <_ 2a sin(v/2) 
fj^(r) =[2r - I^{a-y^)- [ I^(r) - I^(a-y^)]/a} 

+ q(a-yp] A* ; 2a sin (vA) ^r 

=[2r - I^Ca-y^) + [ I^Ca-y^ - /a> 

+ q(a-y^)]A* ; 

where 

g(x) = exp [ -(£n z - 9^)V(2e2)] / ( /Sireg x) 

rt ’,2 2 

q(x) = s [ 2a ( 1-cos v)-r + x {r /a - 4a( 1-cos v) } 

+ x^(3-2cos v) - xVe ] g(x) 

^ / / “2 2” ‘.“2~ 
axud s = r/ V r - a sin v 


(A5.8. i.70) 

(A3.8.1.71) 
(a3.8. 1.72) 

(A3.8. 1.73) 

(A3.8.1.74) 

(A3. 8. 1.75) 

(a3-8- 1.76 ) 
(a3-8,1,77) 
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APJMDIX A3. 8. 2 

probability op GISIIE/ilZED PEjffi SEiailC IHTfflSITY, Y 

a) Magnitude of earthquakes is cutoff at lower lewel i.e, 

m e (m^, “>). 

The jrobability of the eweat [Y > y ] is to be derived for the 
defined probability density fmction of magnitude, fj,(- )> Sq. (3.3) 
and for different hypotheses of focal depth distribution Eqs. (3.6-3. 8). 

The attenuation law, Eq. (3.IO), can be v/ritten as 

Y = SW (A3.8.2.i) 

where Z = exp (C^M) (a3.8.2.2) 

and W = exp[ -C^ i!,n(E + C^) ] (A3.8.2,3) 

according to section 3.7. The joint probability density function of 
Z and W is given by Eq. (3.4) as 

= i^iz) f^(w) (A3.8,2,4' 

The probability density function, - )? of 2 is obtained 
by substitution of Eq, (3.3) into Eq, (3.15) and is given by 

f Jz) = esp [ 6m + 6/G_ £a C + to B to G -(6+C )/C iln z] 

; z € (z^> ") (A3.8.2.5 

where e = exp f C„m ] Ca3. 8.2,6 

1 1 i- 2 o -• 

The probability density function, f^(- ), of W is ^vm by (3.14) 
as 
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= 1/ [CjW 


( 0,4-1 )/h 


j fgCW 


■I/C3 _ 


C4); 


w e (w^.Wg) 


viiietre 


^ 



(A5.8.2.9) 
(A3.8. 2,10) 
(A3.8.2.11 ) 


and is the maximuEi focal distance and fjj(. ) is the probability 

density function of focal distance, R. The details of probability 
distribution function, Pg(,), and the probability density function, fg(. ), 
of focal distance, R,are gi-ven in Appendix A3,8.1. 


To find the probability of the event j^Y > y] , there exist three 
different regions in the domain of y, v/hich are shown in iig. (A3.8. 2. 1). 
Erom this figure and Eq. (3.8) the required probability of the event is 
obtained as 


p[Y>y] = SyCy) =1 ; 0 1 y 1 z.,w^= (A3.8.2.12) 

Eor ye[ z^w^] frcsn Rig, (A3.8-2 .i) and Bq. (3*l9)y the 

required probability is obtained as. 


y/ z, " 2 ® 

p[y > y] = Sy(y) ~ j ^ dw / f2^^(z,w)dz + j dw / 


1 


y/w 


yA^ z^ 


Substitution of Eq. (A5.8.2.4) and Bq, (A3.8. 2*5) into this and integrating 
with respect to z results in 

sr/z 

w. 




SyCy) = r ^ S (y/w) f^/w)dw + / f,/w)d' 

y/z^ 


(A3.8.2.13) 


where 


S^Cu) = exp [8{ia^+ Vc^im C^- an n)}] 


(A3.8.2.14} 
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Sibstituting r = 


R 

o 


Sy(y) = / 

X, 


-1/C3 

w into Eq. (A3. 8.2.11) and integrating gires 

C, 

SgECr + C^) ^y]fg(r)dr + 

» ^1 1 y 1 ^1^2"^ ^2 (A3-8.2. 15 ) 


^ere x^= exp [ (^2^^+ Jin (A 3 . 8 . 2. 16) 

5br y from Eig. ( A3. 8. 2.1 ) and Bq,. (3.19), the reqioired 

parobability is 

w 

2 00 

I* ["^ > y] = Sy^y) = J dw / f (z,w)dz 

y/w 

Similaor jrocedure and substitution as above results in 
o C 

Sy(y) = / S^C (x + C^) ^ yj fp(r)dr ; y > yg (A3. 8.2.17) 


Por different hypotheses of focal depth, H, the required probability 
density and probability distribution function of focal distance, R , from 
Appendix A3. 8.1 to be substituted and final form can be obtained, 

(b) Magiitude of earthquakes is cutoff at both ends; i.e. m e 

!I!ae probability density function, fg(. ), of Z is obtained by 
substitution of Dq. (3. 4) into Eq. (3.15) and is 

f„(z) = K exp[ 8 m + S/C„ JtnC + £n S - £n G -(8+G )/c Jin z] 

ii O c 1 c c <- 

; z e (z^jZg] (A3.8.2.18) 

nhere 

2^ = 0^ exp (Cginp (A3.8.2,19) 
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Ihe probability of the event [Y >_y] is obtained frcm Fig.- ( A3, 8. 2. 2) 
and Bi. (3.19) as 


E [Y > yj = S,^(y) = 1 


y < Z W = y 
*^-11 


(3.8.2.20) 


From Fig. (a 3. 8.2.2) and Eq, (3-l9), the probability of the 


event [ Y ^y] \dien J ^ ^min “ ^2^1^ gi-ven by 


y/ 2 z 

j/ 2 


•'12 2 2 
E [Y > y] = S (y) = f dw f f^„, (zjw) + j dw f f^Jz,7ijdz 

y/w y/s, 


Substituting Eq. (A3. 8.2.4) and B:i. (A3.8.2.18) into this and integrating 
vath respect to z results in 


Sj(y) = J S^iy/v) f^^(w)dw + / f„,(w)dw 

y/z^ 


(a3.8.2.2i) 


^ere 


=K {(C^x) exp (-6m^)}exp (em^) 


(A3. 8. 2.22) 


*/ ^*2 

Substituting r = w ^ into Eq. (a 3.8. 2, 2l) and integrating 

yields 


° - % 

Sy(y) = / [(r + C^) ^ yj fg(r)dr + Fg(xp 

^1 


(A3.8,2.23) 


where 


exp [ (an G - an y + C ^ )/c, j - C 


(A3.8.2.24) 


Using the definition of S_(.) Bq. (A3«8,2.22) into Eq, (A3.8, 2.23) 


results in 
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o C 

Sy(y) = 1-K[ l-Pjj(x^)] + / S* [(r + C^) \]fj^(r)dr 

*» ly 1®^'^ (y2»y3) (A3.8.2.25) 

viiere 

S*(x) = K exp [g {m^-(£G x - nn C^/C^} ] (A3.8.2.26) 

and Pjj(. ) is the probability distribution, of focal distance, R. 

Using Bi. (3»19) and frcm lig. (A3.8. 2.2), the probability of 
[ Y ^ y ] ■vdien ±72 ^■'ren. by 

y/z., 

Sy(y) = / dw / f (z,w)dz + / dw / f2.^y(z,w)dz 

y / z ^ y/w ’ y/z^ z^ 

Adopting similar procedure of Bj, (A3,8. 2.25) fcr tMs yields 

*2 C 

5y(y) = (l-«) +K PjjCxp + J S*[ (r+C^) 5y] f^(r)dr 

’ ^3 1 ^ 1. ^2 (A3.8. 2.27 ) 

where 

Xg = exp [(C 2 m^+ £n G^- y)/C^] - (A3.8,2.29) 

Por y 2 f.y < y 3 (A5.8,2.2) and (3.19), the 

probability of [Y ^y] is 

z„ 

2 2 

SyCy) = / dw / (z,w)dz 

w y/w 
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Using the steps as in Bi. (A3.8. 2, 25) reduces to 

R 

0 ' 

Sy(y) = 14C + J [(r+C^) y l < y (A3.8.2.29) 

Itor y ^max (y^,y^)f the suae probability from KLg. (A3.8.2. 2) 
and Eq, ( 3 , 19 ) is 

2 2 

SyCy) = / dw J t^^(zfw)dz 

y/z^ y/vf 

By similar method of Bq. (A3. 8.2.25) this resiilts in 

^2 C 

S^Cy) = (l-K) Rj^Cxg) + J s| [(r+c^) Vl fK(r)dr 

0 

; y ^ max (yg, y^) (A3.8,2.30) 

One of Eq. (A3.8.2.27) and Bq. (A3. 8.2.29) will conform to tte 
numeric^ values of relative parameters. Above reaalts and the restCLts of Appen- 
dix A3B.1 will produce probability of generalized peak seisouc intensity 
for the combination assumed distribution of focal depth, E, and 
magnitude, M, 
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APPMDIX A5.1 

SECTION OP A SPHEBE (EARTH) BY A PLANE (FAULT) 

Bae equation of a faifLt passing throu^ liie points (a,9^,<j>^) 
aMi(a, 02 j^ 2 ^ surface of the sphere and passing throu^ the 

origin of sphere is given by 

A cos 0 + m sin 0 + n tan <{) = 0 

where 


£ 

= (sin 0 ^ 

tan ~ sin 


tan 

<i>^)/A 

(A5.1.1) 

m 

= (cos Sg 

tan - cos 


tan 

4>2 )/ A 

(A5,1 .2) 

n 

= l/A sin (Sg-e^ 




(A5«1‘ * 3 ) 

A 

= / sin^( 62 ' 

p 

- 6 ^) + tan 0 ^+ 

tan' 

? 

■^ 2 - 

2 tan <j>^ 

tan 4'2 cos ( 62 “®i) 







(A5-.1.4) 

Taking OZ* in Hg. (A5 *i), 

the 

new 

axis of 

Z, the normal to 


the fault plane which passes throxigh 0 and nekes an acute angle with 
OZ, the eq-ustion to OZ'jreferred to Ox, OY, OZ, are 

1/a cos ^ cos e = i/m cos ^ sin 8 = i/n sin (A5.1.5) 

(fekiag OY’, ti® new aals of Y, the line in the plane ZOZ’ 
vdiich is at ri^t angle to OZ' and makes an acute angle with OZ, and 
choosing OX*, the new X-axis, at ri^t euagLe to OZ* and OY*, so that 
the system OX*, OY*, OZ' can be brought to coincide with. OX, OY,OZ» 
The gl-ven plane is X*0Y', and since OX* is at right angles to OZ* 
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and OY’, it is at right angles to OZ which lies in the plane Y'OZ’. 
Hence OX' lies in the plane XOY, and therefore is the line of 
intersection of gi-ven plane and the plane XOY. The equation of the 
plane Y'OZ' is 


l/jl cos e = l/m sin 0 


(A5-1.6) 


SLnce OY' makes an acute angle with OZ, the values of 
2.', m', n’ - the direction-cosines of OY’ are obtained as 


2' = - 2n/ / r+ 


m 


(A5.1.T) 


2 2 

m* = - mn/ / 2 + m 


(A5.1.8) 


» / «2 2 
n' = f 2 -i- m 


(A5.1.9) 


are 


The direction-cosines of OX* — at right angles to OY’ and OZ’ 


= -in/ / 2,^+ 


m 


» = 2 / / 2 + 


m 


(A5* I.IO) 
(A5.1.11) 


n" = 0 


(A5.1.12) 


Hence the transformation is given 

X’ = T X (a,5.1.13) 

where 

X* = <cos <1 cos 0 cos 4) sin 8 sin #> (45,1.14) 



4 
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X’ = <cos cos e’ cos <{i' sin e’ sin 4>’> 


and 

r -a/ / £^+ 


-.£n// £^4- 
£ 


a/v'£^4jn^ 


2 2 
-mn/A + m 

m 


0 


2 2 
/ £ + m 


(a5. 1. 15) 


(A5.1.16) 


This transfoimation is used to rotate a feult such that its 
position after station is equator. 



APPEUEEX A5.2 


PROBABILI!FY USPRIBUPION PQNCPIOI MD PROMBILIPY DENSITY 
HMCTIOR OP POCAC, DISTANCE, R, POR AREA SOURCE 

a) Bocal depth, H, has probability density function, fg(. ), and 
probability distribution function, Pj|(-) iu (0,H^j 

iSie part of the feiJlt vshich influences tte seismic intensity at 
a site can be represented as shovsn in Mg. (A5-2. l) without loss of 
generality. The site 0 is at a distance out of plane. The 

transfcsrmation required to represent a fault plane as in Pig. (A5,2, l) 
is given in Appendix A5.1. The focal distance, R, from site to a 
point at source is given by 

E = 

The probability distribution of focal distance, R, may be obtained by 
finding probability distribution of S and transforming it according 
to Bi. (a 5.2. l). Bie location of earthquahes in the fault plane is 
assumed to be uniform in the direction of 0 . Hie focal depth 
posseses a density fg(.) implying that the distance, D, fron the 
centre of earth has a density 

fjj(d) = fjj(a— d) ; d e ^a — H^,a) (a5.2. 2) 

Hie probability distribution of S, assuming statistical independence 
of 0 and H may be given as 
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i-gCs) = 1 [s <. s] = 1/(K6^) // u f^(u) de du (A5.2.3) 


with 


K=E[D] = a- E[H] (A5-2.4) 

g^= arctan (y^/x^) - arctan (A5.2.5) 


and ® [• ] represents expected value of the variable, Bie 
probability distribution of focal distance, R, is obtained using 
Eh, (a 5.2. l) and Bq. (5.2.3) 


Ejj(r) = Eg ( (A5.2.6) 

arcm Eig. (A5. 2. 1 ), the required probability distribution, 






13 




0 


^■5 1^1 = ^ 


(A5.2.7) 


vdiere 

0=/x^ + y^ (A5.2.8) 

= arctan (y^/x^) - arctan (y^/x^) (A5.2.9) 

Oorresponding probability density is 

j^(r) = 0 ; z^5r_<r^ (A5,2. 10) 

Ecr the region £ s <_ ain (s^jS^), the probability distribution 
of S frm Hg. (A5.2. l) and Eh, (A5.2,3) is 

Fg(s) = l/(Kfe,) J u fjj(u) J de du 

h 
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vhere 

g* = arccos [(C^+ u^~ s^/2Cm3 

2 2 2 

u = C cos S- s — C sin 6^ 

1 T ‘ 

_ _ - 

Ug = C cos 8^ + / s - C sin 
82 = 0 sin $2 

= / a^+ C^- 2aG cos 
= arctan (y^/x^) - arctan (y/x^ ) 


(A5.2.11) 


Integration vdth respect to 8 and use of Bi. (A5.2.2) leads to 

""2 

Pis) = 1 /(K 8 _) J u(e*- B ) fjj(a - u) du 
Using the transformation Bq, (A5.2.6) in this yields 


Pg(r) 

2 

= l/(KB,) 1 (6-6 ) u fjj(a - u)du ; 

J + 

(i^.2. 12 ) 

whore 

,2 2 ~i„ 

t^ = C cos 8^ - /r -z^ - C sin 8^ 

(A5.2.15) 


2 2 2 2 

t_ = C cos 8 + /r -z. - C sin 8^ 

(i5.2.l4) 


8 = arccos [(C^+ z^ + u - r )/2Cu] 

(A5.2.15) 


,22 2 
rg = / C sin 02 + 

(A5.2.16) 


22 2 
r„ = / a + c - 2aC cos z^ 

(a5,2. 17) 
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Differ eat iatioa of Eq, (A5.2, 12) leac3sto probability density of R as 
*2 

fj^(r) = f g(u) fjj(a.^)du ; £r _< min (r^jr^) (A5.2.18) 

wiKre 


g(u) = 2ru/[K8^ / {r^-z^-Cu+c)^} {(u-C)‘^+ - r*^} J (A5.2.19) 


2 2 2, 


aLmilar procedure yields probability distribution and probability 
density of focal distance, R, for various regions as mider; 

if ^2 — ^ — ^3 

,*2 

Pg(r) = l/(Ke ) [ / (6-6 V fjj(a-4i)du + J (b- 6 ^)u fg(a-u)du] 

t V 

1 2 


lA / uf„(a-u)du 


(A5.2.20) 


where 


= C cos - /r^- z^ - C^sin^ 62 


(A5.2.21) 


Vg = C cos 6^ + / r^- z^ - C^sin^ 


(A5.2. 22) 


Corresponding probability density is given by 


1 2 

f^(r) = / g(u) fg(a-4x)du + / g(u) fg(a-u)dt 

1 2 


* T < T < T 

* "^2 ^ 3 


(A5.2.23) 


If r„ < r < r., the probability distribution of R is given by 
3 “• 2 


P^Cr) = i/(KB^) j (B- 6 ^)u fjj(a-u)du 


(A5.2.24) 
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Differentiation of thia yields probability density of R as 


a 


g(u) fjj(a-u)du ; 


(A5.2.25) 


/ 0 o o 

Far max < r < r = /a + C + z,, - 2aC cos 6^ , the 

2 ^ 5 — 4 1 2 ^ 

jjrobability distribution of R is gi*^n by 

*^1 S' 

^(r) = i/CkB-) [ J Ce-6^)u f„(a-u)du + / (e-B,)u f„(a-u)(iu] 

*1 ’^2 


XL. 

+ lA / u f (a-u)du 


1 


(A5.2.26) 


Corresponding probability density is 

T 

1 


= / g(u) fg(a-4i)du + J g(u) fjj(a-4i)du 


; max ^r^ 

The probability distribution of R fbr £ r <_ r^ 
/(a - H^)^+ C^+ - 2C(a - H^) cos is 


(A5.2.27) 


1 a 

P (r) = i/(K& ) J (S-Sj uf (a-u)du + t/K J uf (a-u)du 
K ^ t In ' n 

1 1 

(A5.2.28) 


Probability density is given by 

V 

1 

= / g(u) fjj(a^)du j £ ^ < 


(A5.2.29) 


2 2 2 

Probability distribution of R for =/(a-H^) +0 +z^-2(a-H^)c cosBg 


is giyan by 
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= l/CKg^) 



(e-S^V fjj(a-u)du + lA 


cl 

/ uf jj(a-u )du 


(A5.2.30) 


OorrespoMing probability density is given by 

= J g(u) fg(a-4i)du ; ^5 

a-E ^ 

o 


1^6 


(A5.2.31) 


The probability distribution and density of R beyond Tg is given 
by 


Pj^(r) = 1 

= 0 


r > r. 


(15.2.32) 


b) Eocal depth, H, is lognorioal and mixed lognoimal in (0»H^ j. 

The probability distribution and density functions of the focal 
distance, R, fbr any distribution of foc^ depth, H, has been given in 
the previous section. The foimilae given in the previous section have 
to be integrated for the given probability density, fj|(»)* ibwever, 
the focal depth is assumed to be either lognormally distributed or mixed 
lognormally distributed. Itor both these cases analytical integration 
is not possible. Hence numerical integration technique has been emplojred 
to obtain ttie reqxiired probability distribution and density functions. 

c) Bocal depth, H, is uniform in (0,H^ }. 

Bor uniform probability distribution of focal depth, H, Bq. (a5.2, 4) 

yields 

K = a - H /2 

o 


(A5.2.33) 
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Tbs probability distribution of R is obtained using uniform 
distribution of H and integration fbrmulae for different regions 
gi-yen in section 15. 2a. These are 

S_(r) = 0 ; z < r < r = /c^ sin^ 6 ^ + z? (A5.2.34) 


= l/S p^(r,S^) <_r <iQin (A5.2.35) 

Pg(r) = i/E[Ti|(r,B^) - 1^3 (A5.2.36) 

Pjj(r) = i/E [n 2 (r,e^,a) + n^(r,e^) - n^(r,a)] ; (A5.2.37) 

Pjj(r) = i/E[ n^(r,e^,a) + n^Cr,^^) - n^CrjBg) 

- p^(r,a)] ; maz 1^4 (A5.2.38) 

Rg(r) = i/e[ n^(r, 6 ^) - n^Cr^gg)] + G ; r^ <_r < (A5.2.39) 

Pj^(r) = i/E[Ti^(r,a-H^) - n^CryP^) - a-H^ )] + G 

5 ^5 1^6 (A5.2.40) 

Eg(r) = 1 ; 3^1rg (A5.2.4 i) 

where 



(A5.2.44) 
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n^(r,u,v) = [arccos ^ (G^+v^+z^-r^)/(2Cv-)} -6^ ] 


-/{r^- 2 ^ - (v+c)^}{(v-C)^+ z^-irS /2 

(a5,2.45) 

/ 2 2 2 2 

n^(r,u) = C[ C cos u -»r -z^ - C sin u] sin u 


/2 2\r • ( f n ' ^ ^ 

+ (r - z^ ; [arcsin i{C sin u + ^r -z^-C sin u cos 

u)//r^-z^ 


(A5-2.46) 

n^(r,u) = (r^- z^ ) arcsin { (c^+r^- 2 ^-ai^)/( 2 G’^r^-z^ )} 

(A5.2.47) 

G = aV [ (2a - ] 

(A5.2.48) 


Differentiation of Eqs. (A5.2.34 ~ A5-2.4l) lead to probability 


density of E for 

Yarious regions. Biey are the folio-wing 


fgCr) = 0 

; r < and r > 

^ - 1 -*-6 

(A5.2.49) 

I! 

C 5^(r,Bp - 52^r,B^)] ; r^< r <_ minCr^jr^) 

(A5.2.50) 

fg(r) = iA 




f ^2 - ^3 

(A5.2.51) 

f^(r) = 1/3! 

[5^(r,Bp - 5^(r,a)3 i 1 ^2 

(A5.2.52) 

fg(r) = 1/3} 

Cs^(^»B^) - £^{r,a) - ^•j(^>®2^ ^ 



i max < r^ 

(A5.2.53) 

fjj(r) = t/T [ C^(r,g^) - C^CrjSg)] j 1 

(a5.2,54) 


:^(r) = 1/T [ ^ir,a-E^ ) - Bg)] ; 1 ^ 1 ^6 


{A5.2.55) 
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where 


0? = /r^- /[ (a - H^e) 3^ ] (A5.2.56) 

2 2 2 2 2 2 2 

(r,u) = arcsin[ (c sin u +/r - z^- C sin u cos u)/ /r - z^ ] 

(A5.2.57) 

^2 (r,u) = arc sin [ (c sin^-. _ c^sin^u cos u)/ /r^-z^] 

(A5.2.58) 

(r,u) = arcsin[ (r^- z^ - u^)/ (2C/r^- z^)] (A5.2.59) 



APPafDIX 16.4. 1 


THE probability IHSffilBUTIOH HJHCTIOH OF Z 

a) The Ife^iitude, H, is cutoff at lower level i.e. me (O, ») 

The probability distribution function, P„(. ), can be easily 

u 

obtained frcm Fig. (a6.4.i) as 

= 0 ; z 1. exp (x^+ 0 ^ m^) = (a6,4.1.i) 

The probability distribution function of z is given by 
P [ Z <_ z ] = P [ X + C^M £ Jin z ] 

Assuming statistical independence of X = S-n and M and Fig. ( 16,4.1.2) 
3 d.elds 

( 2nz-x )/ Cg 

= J I fJ^:) 4-Cy)<Jy 

" X. m 

1 0 

Infegration of this results in 
jyaz-C^m^ 

P (z) = / f (x) Pjj {(£nz - xyCg} dx 

^1 

Integration aM sln^iLification of this results in 

P^(z) = ilj) + ^ 

; z^ 2 1 ^ 2 ^ ^2^0 ^ 


Ca6» 4« 1 . 2 ) 
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where 


q^(u,v,w,y) = (to u-v-C^w) (to u-v-C^w - 2C^&)y 

(a6.4. 1.3 

= 2y Pj^(v) 

(a6.4. 1.4, 

1/A^ = ix^- ) (xg- x^ ) 

(a6.4. 1 . 5 ) 

= {m z - x^ )/c^ 

(a6,4. 1,6) 

Fjj(x) = 1-exp [e(iB^-x)] ; x^m^ 

(46,4.1.7) 

proceeding alcaig atoilar lines the probability distribution 

of Z is obtained in different regions as xjnder 

= q^(z,x^,in^,A^) + q^iz,e^,h^) - q^(z,X2,m^, A_^ ) 

function 


(a6.4.1,8) 

^ 2 ( 2 ) = l+qgCzje^jA^ - + qg(z,e^,A^)j z 5 

vdiere 

^ (-46, 4, 1 . 9 ) 

= im z - ; i = 1,2,3 

(a6.4,1. 10 ) 

1/42= ^-^1^ 

(A6.4. 1 . 11 ) 

^4,= (*,-x,)(Xj-X2) 

(a6.4. 1 , 12 ) 


b) The ilagnitude, M, is cutoff at both level i.e. Me (m^ ] 

Tbs d(®jaiQ of M and X are shown in Fig, (a 6.4. 1,3). Btom 
Hg, (Afi,4. 1 . 5 ) it is easily seen that probability distribution function, 
is 
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=0 ; z <_ z^= exp (x^+ (a6.4. 1-13) 

Hie probability distribution as well as density fmction of Z in the 
region ^ z ;^min (z^jZ^) is obtained from RLg. (a 6.4. 1,4) and the 
mutual statistical independence of C_^ and M as 

l[Z^z] = p[X4-C2M<.mz3 

Anz-C m ( ilnz-x )/ C_ 

cl O c 

or P (z) = / / f (x,y)dy dx 


gaz~0^m^ 

= K J f^(x) P,^ { (£az-x)/C 2 } dx (A6.4. 1.14) 

Xi 

vshere ^/K = 1 - exp [ B (m - m-) ] (a6,4. 1.15) 

^ O 1 

SabstitutiOE of f,^X* ) f^om S|. (6*22) into (a 6*4*1#14) and 

A 

subsequent integration and simplification results in 


Fg^z) = K [ q^(z^x^r A^) + ^i^^ 

; z^ ^z _<min (z^j z^) (A6.4.1. 16 ) 


where 

Zg = espj (xg + (A6,4.1.17) 

z^ = exp (x^ + G^m^) (A6.4. 1.18) 

She probability distribution of Z for -rarious regions is 
obtained by emplc^ng similar procedure as in Si. (A6.4.1.16), giving 
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•. =4i“l’“2 CA6.4.1.19) 

(46.4.1.20) 


F,(») 


^2 S' 


^ U) = i^A -Kq.,(2,X5»m^,&5) + K^q^(z,X5,m^,^3) 

; max (121=^5) 

F (z) = 1 -«.+ K[ q 2 (z>e.|,A^) + - q 2 Ue 2 »^ 2 ^^ 


(A6.4.I.21) 


(A6,4.1.22) 


5 Z3 1 2 1 ^4 


(A6.4.I.23) 


F (z) = I+OA -K.,q.,(z»x^,m^,A.j) 

Z 1 ' 

- <l^(z,^2,h'> ■" 92 (‘'"' 3 ’‘ 3 ^ ’ - “5 

Fj,(^) =. 1 - ui,+ K,q,(z,*5>i»i.‘i) + K<l2(».'3’‘3^ 

. Z^ 1 2 : 1 2 ig 


5’_(z) =1 f 2i L *^g 

Z 


(A6.4.I.24) 


(A6.4.1*25) 

(A6.4.1.26) 
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?fhere 


=K - 1 

= exp (x^ + ^2 %) 
= exp (xg + ®2 

^6 = ^ °2 


(a6.4. 1.27) 
(a6.4. 1 • 28) 
(a6-4. 1-29) 
(A6. 4-1.30) 
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APPKIDIX A6.4.2 

!EHE PROBABILITY ISa-ISITY AKD DISTRIBUTION FUNCTION OP W 

s 

She domain of W = (e + C^) ^ is shorn in Fig. (A6.4.2. l). 

The probability distribution function of W is gi-ven by 

P[W^w] = P [G^iin(R + C^)£ £nw] {A6.4.2. 1 ) 

and is evaluated using statistical independence of R and G^. Prom 
ELg. (A6.4.2. l) it is obvious that tte probability distribution function, 

), and probability density function, f^(. ) are given by 

a 

P,iy(w) =0 {A6.4-2.2) 

f^(w) =0 (A6.4.2.3) 

For ^ w ^min the probability distribution frc» Bl. (A6,4*2. l) 

and Fig. (A6.4.2.2) is given by 

P^y(w) = J / f^ (y) fjj(r) dy dr 

o 3 

Integration with arespect to y and use of Bq, (6.23) lead to 

P,y(w) = / g(w,r,a^, xp fj^(r)dr ; <_w min (wg,®^) (Afi.4,2.4} 

' o 

Aere 

2 

g(x,y,z,a} = {tn x/£n(y+0^) - z) a 


(Afi,4.2.5 ) 
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l/X 


1 



= (a^- ) (a^- ) 



(A6.4.2.6) 

(A6.4.2-7) 

(A6.4.2.9) 


and fg(.) is the probability density of focal distance, E. Uiffer- 
entiation of Eq,« (A6.4.2.4) with respect to r yields the probability 
density function of W as 

f^(w) = J h(w,r,a^, X^) fj^(r)dr (A6.4.2.10) 


vdiere 


hix,jfZ,a) = 2a [to 3i/to(y +G^ )-z ] / [ x £n(y + C^) ] (A6.4.2. 11 ) 

Hie probability distribution function of W for other regions are 
obtained adopting similar procedure as in Bi. (a 6.4.2.4). Hiese aie 

j, j. 

1 2 

5\y(w) = ^ ^ g(wjr,a^,X^) fj^(r)dr J g(w,r,aj, X^ )fjj(r )dr 


R 


• Wg ^ w £min(wj,w^) (A6.4.2. 12) 


5’,y(w) = J g(w,r,a^,X^) fjj(r)dr ; 1 w £ minCwgjW^ ) (a6.4,2. 13) 


E 


E. 


^(w) = Pg(rg) + J g(w,r,a^,X^) fg(r)dr - J giwjrjS^jX^) fg(r)dr 


; iaax{wg,w^) w ^min (w^,w^) (a6,4.2. 14) 


S 


%(w) = 1- I g(w,r,aj,Xj) fg(r)dr j maxCwgjW^ ) < w <_iain(wj,Wg) 


(a6,4.2. 15) 
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F^(w) = J'jjCrg) - / g(w,r,aj,X^) fjj(r)dr 

^3 

+ J g(w,r,a^,X^) fjj(r)drj w^<_w (a6.4.2. 16} 

^2 
^2 

P/w) = - I gCw,r,aj,Xj) fj^(r)dr 

^3 

H 

o 

+ I g(w,r,a^,X^) f^(r)dr; niax(w^,w^ ) ;^ w <_ (a6.4.2. 1?) 

Pg 

H 

O 

P^(w) = 1- J gCwjPja^jXj) fjj(r)dr ; laaxCwjjW^ ) ( a6,4«2*18) 
^3 

Fyy(w) = 1 J (A6.4.2.19) 


where 


l/Sg 

=. w - 

(a6. 4-2.20) 

''^®3 

r- = w - C. 

3 4 

(a6.4.2,21) 

lA. 

II 

(A6.4.2.22) 

"5 = ‘V 

(a6.4-2,23) 

"6 = 

(A6.4.2-24) 

t/Xj= (a^-a^)(a^-ag) 

(A6.4.2-S) 


and !„(. ) is the probability distribution function of focal distance, 

S' 

Mm and R is the majxmm tocBl distance.* 

.0 
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Different iatioa of the Eq, (A 6 . 4 . 2 . 12 ) to Bq. ( 6 . 4 . 2 . 19 ) leads to the 
probability density function of W and they are given in the followinf^. 

1 2 

%(w) = / h(w,r,a ,X^) fj^(r)dr - J h(w,r,a_,X^) fjj(r)dr 
" 0 


; ^2 (wj,w^) 


(A 6 . 4 . 2 .^) 


R 


%(w) = / h(w,r,a^,X^) fg(r)dr; w^<_ w <_iBin(w 2 »w^ ) (A 6 . 4 . 2 . 27 } 


R 


f^(w) = j h(w,r,a^,X^) fg(r)dr; maxCw^jW^) <, oinCw^jW^ ) 


R 


(A 6 . 4 . 2 . 28 ) 

f^(w) = - J h(w,r,a^,X^) fj^(r)dr; maxCw^jW^ )<_w^min(w^,Wg ) 
o 

(aba, 2,29) 

^1 ^2 

f^(w) = J h(w,r,a^,X^) fjj(r)dr J hCwjrja^jX^) fg(r)dr 


. w, < w < w 

7 3 — — 4 


(A 6 4 . 2 . 50 ) 


R 


* / ii(w,r,a^,xp fj^(r)dr - J h{w,r,a2,Xj)f jj(r)dr 


^2 ' ’ ^3 


; max (wj,w^) <_Wg (A6.4.2.3 i) 

R 

o 

f^(w) = / hCwjrjaj,!^) fjj(r)dr ; max (w^,w^ )<_ w < Wg (a 6 , 4 . 2 . 32 ) 


f^^(w) = 0 


otheiwise 


(a 6 . 4 . 2 . 33 ) 


Only some of the foroulae given above will be valid for the 
numerloal values of the parameters involved. 
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APPENDIX A6.4.3 

PBDBABILITY OP PEAK GfflER/iLI2S)D SEI3IIC INTENSITY, Y 


The probability of [ Y > y ] of peak generalized intensity is 
obtained using the distribution laws for W and Z derived in 
Appendix A6.4.1 and A6.4.2. Since M,R, and are nrutmlly 

independent so will be z and V/, Hiis probability for magnitude 
cutoff at lower ends and magnitude cutoff at both ends are given in 
the following. 

a) lfei@3itude cutoff at lower end i.e, m z (m^, “ ) 

The dcanain of Y is shown in Pig. (A6.4.5. l) 

The required probability 


5CY>y]=1 »yi zy'wg = y^ 

where - exp (x^+ 

K 

and Wg *= (Hq+ ^ 


(A6.4.3.1) 
(A6.4.3.2) 
(A6, 4*3.3) 


This probability is obtained for y^ 
(A6.4.3.2) and is given by 


w. 

6 wy 

!’[ Y> y] = Sy(y) = 1-1 dw / f ^(z,w)dz 

z^/y 

Use of statistical indeperaience of Z and W and subsequent 


integration with respect to z results in 
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w 


SyCy) = 1 - J f (w) f (wy) dw; y. <,y£ yp 


vsliere 


Sl> -2 

y^ = exp (x^+ ^4^ 

■a 


(46.4.3.4) 


(A6.4.3.5 ) 


(a6«4«3«6 ) 


(Hie probability distribution of Z given in Appendix A6.4. 1 
for suitable region has to be ui^d in Bj. (a 6,4.3.4). However fto?ther 
sinplification is not possible analytically. lumerical integration 
has to be used to integrate the integral in Eq, (a 6,4.3.4). Similar 
procedure lead to 

Sy(y) = 1 - / P^Cwy) f^(w)dw ; y^yg Ca6,4.3.7) 

b) Magnitude cuttoff at both ends i.e, m e (m^, a^) 

The domain of Z and W are shown in KLg, (A6.4.3.3). 

Poliowing similar procedure as in previous section the i>equired 
probability is obtained as 


Sy(y) = P[l>y] *1 ; yi “ 2;^Wg 


\(y) « 1 - / ^^(w) f/w)dw ; y^<_y^iain(y2,y^) 


Sy(y) = 1 - / F^(w) f^(w)dwj yg ly 1 y^ 


(A6. 4.3.8 ) 
(a6,4.3.9) 


(A6.4.3. IO) 
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#3© re 


Sy(y) = ^(z/y) / %Cw)<aw5 

zg/y 

Sy(y) = P^(zg/y) - J P^^wy) f/w)dw; max {y2,y3)iy4 


(A6.4.3 .II) 


(a6.4.3» 12 ) 


S^(y) = 0 otherwise 


— 

y-i “ ""As “ ^ °2“o ^ ^ V ^4^ 

(A6.4.3. 13) 


(A6.4.3. 14) 

^3 “ V*6 " f V °2 ®o3 °4^ ^ 

(A6.4.3.15) 

■^1 

(A6.4.3.I6) 
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Appendix A7. 1 
Table A7.1 

Lisft of Earthquakes before 1917 


Bate 

hr 

Time 

- min - Sec 

Latitude 

degree 

Longitude Eocal 
degree depth 

IGas. 

liagnitude 

i RemaxkB 

08930000 


- 


- 

> 7.0 

Biabul (vin- 
kaown) 

15050706 

- 

H 34.45 

S 69. 16 

- 

> 8.0 


16680502 

- 

- 

- 

- 

> 7.0 

delta of Indus 

16680511 

- 

- 

-• 

- 

> 7.0 

delta of Indus 

16690604 

* 

- 



> 6.0 

Ifedrem 

17200715 


I 28.39 

E 77. 13 


6.5 

4-5 sfaocks/day 
for 40 days, 
occasional 4 
to 5 months. 

17620402 

- 

I 22.24 

E 90,42 


> 7.5 

- 

17640604 

- 

- 

- 

• 

> 6.0 

Bank of Ganges, 
Several shocks 

18030901 


H 27.51 

E 77.42 

- 

6.5 


18080000 

- 

- 

- 

- 

6.5 

Eumaon 

181«)5^ 

- 

H 30.00 

E 80.00 

- 

6.5 

2-shocfcs 

18190616 


- 

- 

« 

8,0 

Bann of Kutch 

18210815 

- 

H 23.13 

E 70,08 

- 

> 5.0 


18220129 

- 

H 13.14 

E 79.06 

- 

5,0 

Several shocks 

(s.s.) 

18220403 

- 

- 

- 

I>eep 

6,5 

Manipw-Burma 
Border , S. 8. 
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18230209 

- 




> 6,0 

lodiaii Oceaii 

18^1029 

- 

N 27.50 

E 86.50 


6.0 

8 minor after- 
shocks (a.s.) 

18270900 

— 

— 

- 


-* 

In hills-Bavi 
emerges in 
plain 

18280606 


N 34.04 

E 74.49 


> 6.0 


18280708 


- 

- 

Beep 

> 5.0 

Manipur-Burma 

Border 

18280822 


H 14.04 

E 72.85 


> 5.0 


18301231 

- 

N 22.24 

E 90-42 

- 

o 

« 

A 


18310000 


N 30.00 

E 70,00 

- 

6.0 


18320221 

- 

H 36.50 

B 70.50 


7.0 

S.S. 

18330826 

- 

S 27.50 

B 86.50 

- 

7.5- 

8.0 


18331002 

- 

- 

- 

- 

>6.0 

A.S. of 18330826 

18331018 

- 

- 

- 

- 

>6.0 

A.S. of 18550826 

18390323 

- 

H 22.00 

E 96,08 


7. 5 -8.0 


18420116 


iJ 27.31 

E 77.42 

- 



18420219 


I 34-25 

E 70.26 


7.5 

100 A.S. in a 
atoath 

18420305 


I 30.29 

E 78.06 

- 

5*5 


18421111 


- 

- 

De-ep 

>6,0 

Itanipiir-BiaEa 

Bosrder 

18430401 


I 15.12 

175.57 


5. 5-6.0 


18430411 

• 

H 30.00 

B 00.00 


> 5.0 


18450419 


N 23.48 

E 68.47 


6.0 

To 18450425, 

66 A.S. 
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18450806 

S 24.45 

S 91.55 

- 

6-5 


18461018 

N 24.62 

E 72.70 


6.0 


18520124 

H 55-94 

E 75.39 

- 

6.0 


18520500 

N 27.55 

E 81.10 

- 

6.5 






5.0 

Near Kot^ixr 

18560407 





(Simla hills) 

13580824 

N 19.00 

E ^.50 

- 

7.5 


18651118 

H 22-00 

E 75.00 


5.0 


18640429 

H 24.00 

E 70.00 

- 

6.0 


18651219 



Deep 

6.0 

Manipur-SuTB® 




Border 

18660525 

H 27.5 

E 86.5 

- 

7.0 






> 5.0 

Bampur-^uroaa 

18680650 





Border 






Peshavar, I 

18680800 





Milne 





> 7.5 

Near Cachar 

18690110 




(Assam) 






peshwar, I 

18690400 





mine 

1^90707 

H 27.5 

E ^.5 

- 

6.5 







ifewalpindi, I 

18691220 





mine 






Dacca 1 I mine 

187(X)422 





Sind, I mine 

18701028 





Sind, III Mila 

18721215 





Darjiling, I 

13750426 

- 




mine 






labore, II 

18751212 

mm 




mine 
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187^302 

18811231 

18830400 

18850530 
188507 14 
I888l0ffl 

18881223 

18891017 

18921222 

18931105 

18970612 

18990925 

190)0208 
190^)418 
190206 16 

19020822 

190 » 9 » 

1^30114 

19050404 

1905092S 

19060228 



- 

- 

* Punjab, I 

mine 




>8.0 3®-y of Ben^l 


- 

- 

_ peswan, II 

mine 

H 34.6 

E 74.35 

- 

7.0 

N 23-5 

E 90.25 

- 

>7.0 0^ mcca 

- 

- 

- 

Srina^r, I 
mine 




.. Burma, I Milne 

^ - 

- 

- 

Or iSjSeraj^nj, 
Hiarma- I mine 

I 

B 66.25 

- 

- 

mm •• 

- 

- 

-I Milne 

114)6-00 N 26.0 

E 91.0 

*- 

8.7 A. S. for years , 

Dor j iling, I II 

M ^.7 

E 76.7 

70 

Milne 

6.0 

H 36.5 

1 70.5 

Deep 

7.0 

H 31.0 

1 79.0 


6.0 Pell5 

dun at 07 hr 

H 37.0 

E 77.0 


20 sain. 

7.0 

„ H 36.5 

E 70.5 

'Deep 

7.0 

I 24.0 

E 70.0 

- 

6.0 

S 32.5 

S 76.25 


8.0 ’ 

01-26-09 H 29.0 

E 74.00 


7.0 

S 32.0 

B 77.00 


7.0 
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% 


19060520 

- 

H 36.5 

E 70.5 

Deep 

7-0 

190606 13 

- 

H 31.0 

E 7945 

5. 

.5-6.0 

19060815 

- 

N 24.62 

B 72.7 

- 

5.0 

1906083 1 

14-57-30 

N 28.00 

E 94.0 

- 

6.5 

1906 1024 

14-42-51 

N 40.00 

E 68,0 

Siallow 

7.0 

19070413 

- 

H 36.5 

B 70.5 

260 

7.0 

19071021 

04-23-36 

N 38.0 

E 69.0 

3iallow 

8.0 

19071225 

22-36-00 

N 36.5 

E 70.5 

240 

6.75 

1^a)112 

- 

H 32.0 

B 68.0 

- 

6.0 

19080209 

18-13-00 

H 36.6 

E 90.0 

aiallov? 

7.0 

190^3305 

- 

N 32.0 

E 68.0 

- 

5.0 

19080312 

19-26-24 

I 36.5 

E 70.5 

220 

6.75 

1908 23 

20-14-06 

S 36.5 

E 70.5 

220 

7.0 

19081211 

- 

H 31.0 

E 79.0 

- 

7.5 

19081212 

12-54-54 

I 26.0 

E79.0 

- 

7.5 

'19090217 

_ 

I 27.0 

E 87.0 

- 

5.0 

19090409 

_ 

I 24.62 

E 72.7 

- 

5.0 

1^^707 

21-37-50 

K 36.50 

1 70,5 

230 

7.75 

1»90707 

21-39-05 

N 36.5 

E 70.5 

.Shallow 

8.0 

19090907 

- 

H 36.5 

E 70,5 

- 

6.0 

1909 lOX) 

23-41-12 

N 30.0 

E 68.0 

- 

7.2 

13100712 

07-36-12 

I 38.0 

I 66.0 

Shall 0 ¥/ 

8.7 

19110101 

IO— 1843 O 

I 58.0 

B 66.0 

50 

7.2 

19110704 

, , 

N 36.0 

E70.5 

Deep 

7.6 
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7 moderate 
shocks 
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5?able contd, A7. i 

19111207 

19120425 02-24-06 
19120425 10-27-48 

19120523 02-24-26 
19120523 23-08-18 
19120823 21-14-30 
19121128 20-55--06 

19140206 11-42-18 

19140328 10-44-48 

19141011 16-17-06 

19151114 
191512D5 
19160107 
191608» 


H 30.0 

E 88.0 

- 

5.0 

N 21.0 

E 97-0 

- 

8.0 

N 36.5 

E 70.5 

220 

6.75 

N 21.0 

E 97.0 

- 

8.0 

N 36.5 

E 70.5 

220 

6,25 

N 36.5 

E 70.5 

»0 

7.5 

M 36.5 

E 70-5 

230 

6.5 

il 29.5 

E 65-0 

100 

7.0 

n 25.0 

90.0 

100 

6.7 

N 12.0 

E 94.0 

80 

7.2 

N 63.3 

E 92.7 

- 

5.0 

H 63.3 

E 92.7 

- 

5.0 

H 13.0 

E 77.83 

- 

5.0 

» 30.0 

B 81.00 


7.5 


